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Outline 



This thesis describes theoretical investigations of universahty and finite- 
range corrections in few- and many-boson systems. The major part of this 
work concerns ultra-cold trapped atomic gases, but some of the results and 
methods may be applicable to small molecules, and nuclei too. 

In the introductory chapter [T] we briefly describe some of the basic con- 
cepts and phenomenology used in this dissertation such as universality, Bose- 
Einstein condensation, Efimov physics, and Borromean binding. The typical 
experimental systems of interest are also described. 

Chapter [2] reviews the relevant theoretical background and methods for 
few-body systems as well as mean-field models for condensates. The nu- 
merical procedures we implement are also described. This gives a common 
background for the rest of the dissertation. The main chapters |3j|8] can be 
read almost independently, but they often refer back to chapter [2} 

In chapter [3] we investigate universal finite-range corrections to Efimov 
physics in three-boson systems. Connections to Borromean binding are 
made. We also describe the effects of putting the system in a finite trap. 
Chapter |4] continues this line by investigating the conditions for Efimov 
physics, in particular for large effective range. The content of chapter [3] and 
[4] was published in |TFJ08cl lTFJ^09 , but new material is also presented. 



(See page |xiii| for a list of publications.) 

In chapter |5] we show the existence of a many-body Efimov effect based 
on two-body correlations. The features of the effect and experimental sig- 
natures are discussed. This chapter is based on [TFJOSa] IFJTOSj as well as 
some new results. 

In chapter [6] we numerically consider trapped few-boson systems of or- 
der 10-30 particles including two-body correlations. Condensate-type states 
are identified. We investigate energies and correlations for these condensate 
states and discuss possible experimental signatures. Appendix [A] and [B] con- 
tains some technical details used in this chapter. The content of this chapter 
is composed of selected material from |TI^'J07l ITFJOSbi IFJT08llJKTF07| . 

Chapter [7] approaches the question of effective range corrections in con- 
densates from the mean-field point of view. A modified Gross-Pitaevskii 
equation is combined with a higher-order Feshbach resonance model. Effects 
on stability and decay mechanisms are predicted. In chapter |8] we continue 
the mean-field analysis of condensates. We present the Thomas-Fermi ap- 
proximation in the case of higher-order interactions. Some details are found 
in appendix [C} Chapter [T] and [| are based on |TZJ09| [ZT09bl IZT09aj . 

Conclusions and outlook are given in chapter [9] The back-matter con- 
tains appendices and the bibliography. 
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The results presented in FJT"'"09 do not fit in the hne of this thesis and 
is omitted. 

Online information: This dissertation is available online at the address 



http : //www.martint .dk/thesisl as well as on the arXiv preprint server. 



The bibliography for the thesis is also available on the above address, in- 
cluding abstracts and direct (doi) links to all referenced articles. This will 
hopefully be useful for the reader. 
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Chapter 1 

Introduction 



Universality or model-independence is desirable in theoretical and experi- 
mental physics because the scope becomes broader and the applications more 
transparent and flexible. Universality generally refers to systems determined 
by only a few large-distance parameters, implying that they become inde- 
pendent of the specific short-range interactions and structure. Thus, once 
theoretical formulations are established they can be applied to a wide range 
of fields, e.g. nuclear, atomic, and molecular physics. Much effort has there- 
fore been devoted to extract universal features in many branches of physics. 



1.1 Ultra-Cold Atomic Gases 

One system in which universality can be found is Bose-Einstein condensates 
(BEG) where cold massive bosons occupy the same quantum-mechanical 
state. The BEC concept was proposed almost a century ago |Ein24j . but 
remained an elusive goal for atomic gas experiments for many years. The 
quest for this cold grail was motivated by several factors. The idea of a 
macroscopic number of particles in the same single-particle ground state is 
simple, elegant, and easy to grasp. Also, it is direct evidence of quantum 
mechanics and particle-wave duality on a macroscopic scale. 

Since the experimental realization of the first BEC atomic ®''Rb two 
decades ago |AEM+95] the field of ultra-cold atomic gases has picked up 
speed. Bose-Einstein condensation has now been achieved for most of the 
alkali gases in numerous experimental groups, typically with ^Li, ^^Na, 
39k, -^ ^K, ^^Rb, ^■^Rb, and ^^^Cs [PS02]. In the wake of the BEC mile- 
stone AEM^ 95l a tremendous variety of experiments has been realized 



|BDZ08l IGP S08 . This includes the interference between two condensates 
proving full coherence |ATM"'"97j . structured arrays of quantized vortices in 
rotating condensates |ASRVK01j . and propagating solitons and sound-like 
waves [PS03j . Besides standard (harmonic) traps, numerous experiments 
have been performed with different external confinement, including one- 
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and two-dimensional gases and optical lattices with e.g. Mott-insulating 
phases |BDZ08j . Experiments with ultra-cold fermions have also given rise 
to new physics, e.g. BEC-BCS crossover and super-fluidity have been re- 
alized |GPS08] . Moreover, it is possible to produce cold gases with a wide 
range of mixed bosonic and fermionic species |SZS+n4[llG"0+n4l . Produc- 



tion of bosonic gases with a small definite number of atoms is also in progress 
|DRNn7] . 

Currently, one of the most important experimental tools for ultra-cold 
gases is the use of Feshbach resonances |KGJ06l ICGJT09] , which are found 
in virtually any atomic system. This allows for arbitrary tuning of atomic 
scattering lengths (i.e. low-energy interaction strengths) by varying an ex- 
ternal magnetic field. Tunable scattering lengths have been a key ingredient 
for controlling the production, macroscopic stability, and decay of conden- 
sates [DGPS99| . The Feshbach resonance method also allows for production 
of weakly bound diatomic molecules (shallow dimers) }KGJ06| ICGJT09] . 





Figure 1.1: Experimental realization of Bose-Einstein condensation in trapped al- 
kali gases. Left: The characteristic "spike" of condensed ^^Rb atoms emerges from 
the thermal cloud as the temperature is lowered below a critical value. Right: A 
rotating ^"^Na BEC with vortices arranged in a highly ordered structure. From 
ETH Zurich (2002) and MIT (2001) [ASRVKOlj . 



From the perspective of universality the scattering length is the decisive 
parameter for ultra-cold gases. Most features of the above experiments can 
be understood in terms of this parameter alone. This simplifies the theo- 
retical descriptions immensely and makes the interpretations of the physical 
effects extremely sharp and clear. Most of the theoretical models for con- 
densates are based on some sort of mean-field approximation with a con- 
tact interaction proportional to the scattering length. This means that all 
high-energy effects are integrated out and short-range model-dependent cor- 
relations are neglected. Thus it gives simple models which by definition are 
strictly universal. 

When such universal models are used it is of course essential to know the 
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precise "window" where the results can be apphed to experiments. In many 
reahstic cases this window is in fact very narrow. It is therefore important 
to know the quahtative significance of higher-order effects at the border of 
this window. The higher-order effects may ultimately be expressed in terms 
new model-independent parameters, but can also be strictly non-universal. 



1.2 Few-Body Systems and the Efimov Effect 

Another long-standing prediction of universality is the Efimov effect, which 
was presented in the early 1970s [Efi70, EfiTll IAN 711 IAN 721 lEfi73] and 
has been discussed thoroughly since jFR,S79L lFJ93l IFJR,94[ IAFT97I INFJ991 
IFJOlbl INFJGOll IBH06] . It is an anomalous effect which occurs in neutral 
three-body systems when at least two of the underlying subsystems have 
resonant s-wave interactions, i.e. a diverging scattering length (or at least 
much larger than the range of the physical force). This corresponds to the 
subsystems being close to the threshold for binding. Under these conditions 
an exceptional sequence of infinitely many geometrically spaced three-body 
bound states occurs, with an accumulation point at zero energy. Moreover, 
the ratio between the energies is a universal constant depending only on the 
mass ratios and the particle statistics. For three identical bosons the famous 
ratio is |NFJG01| 

= e^-/^" 515.0 . (1.1) 

Hjrp 

We emphasize that this specific scaling factor is independent of the mass, 
and that the same effect in principle can be found at all physical scales. 
The Efimov states have large sizes which are also geometrically spaced with 
ratios 515.0^/^ = 22.7, and lie outside the classically allowed region. 

The effect was initially proposed by V. Efimov within the field of nuclear 
physics and has been searched for over many years in nuclei - unfortunately 
this search has been futile. The only possible three-body candidate with- 
out Coulomb interactions is two neutrons and a (positively charged) core. 
The possible resonant s-wave interactions should then be between each of 
neutrons and the cor^ Such three-nucleon systems, called two-neutron 
halos, do in fact exist. Common examples are ^He {a + n + n) and ^^Li 
(^Li-|-n -|- n) |JRFG04] where the essential degrees of freedom correspond 
to the core-|-n -|- n structure. Unfortunately the heavy-light-light mass ratio 
gives an extremely large scaling factor (as compared to the already large 



factor in eq. (1.1)) jJF03] . Thus one can only hope for "accidental fine- 
tuning" of the neutron-core scattering length to see the Efimov effect. One 
of the few candidates is ^^C |AFT97I IBH06j . The Efimov effect has also 



^Two neutrons with a spin-singlet function can be in a relative s-wave state. However, 
the singlet n-n scattering length is ~ 20fm, which is not large enough. 
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been searched for in atomic systems. There is a general theoretical con- 
sensus |LDD77l IELG961 IEG061 IMFJ98] that the first excited state of the 
atomic helium trimer ^Hea is an Efimov state, although the experimental 
observation so far proved elusive |BKK"'"05 . 

Another concept closely related to the Efimov effect is Borromean bind- 
ing. The phrase refers to atypical bound few-body systems where none of 
the individual subsystems are bound. This is often symbolized by the inter- 
locked Borromean rings, see fig. |1.2[ Borromean systems are typically weakly 
bound with large spatial probabilities in classically forbidden regions. This 
implies that the structure is mainly determined by the binding energy alone, 
and are thus good universal candidates. The Efimov effect can be seen as 
a critical limit of Borromean binding: Here the two-body subsystems are 
exactly at the threshold but infinitely many bound three-body states exists. 
Examples of Borromean systems are the two-neutron halos ^He and ^^Li 
discussed above as well as ^Be {a + a + n) |JRFG04] . 




Figure 1.2: Left: Schematic two-neutron halo nucleus with the neutrons weakly 
bound to the core. Right: Interlocked Borromean rings: If one ring is removed 
the others fall apart. The name "Borromean" is historical and comes from the 
aristocratic Italian Borromeo family who used the rings in their coat of arms. 



1.3 Efimov Physics in Atomic Gases 

The obstacles for observing the Efimov effect in nuclei are three-fold: i) 
nuclei are electrically charged, ii) the are no methods for tuning the interac- 
tions externally so one must rely on accidental tuning, and iii) the remaining 
candidates have very unfortunate mass ratios. Instead, the approach of ob- 
serving the Efimov effect indirectly via three-body recombination losses in 
ultra-cold atomic gases attract much attention currently. As compared to 
nuclei, atomic gases are very promising for several reasons: i) atoms are elec- 
trically neutral, ii) scattering lengths can be tuned by magnetic Feshbach 
resonances, and iii) by using different species with appropriate mass ratios 
(heavy-heavy- light) the universal scaling factor can be reduced. Concern- 
ing the last part, specific proposals for ultra-cold atomic gas mixtures were 
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given in fDE06], for example the realized boson- fermion mixtures ^^Na-^Li 
SZS+04j or 87Rb-40K ^100+04] . or even better ^^^Cs-^Li or ^''Rb-^Li. 



The basic idea for observing Efimov states using Feshbach resonances is 
to measure the bound state spectrum indirectly. When the scattering length 
is tuned around large values, a series of resonance and interference effects 
occurs in the three-body recombination rate |EG06| IBH06] . These features 
are located at critical scattering lengths with the same scaling factor as 
the Efimov sequence, i.e. 22.7 for identical bosons. The first experimental 
indications of such features has been found recently in cold ^^^Cs gases 
jKMW"'"06l lKFM"'"d9 . Also, a universal trimer state is reported in a three- 



component ^Li Fermi gas |WLO+09|. Recent theoretical |DvSG09l lvSDG09j 



and experimental |FKB"'"09| work indicate the existence of a universal four- 
body effect where an atom is weakly bound to an Efimov trimer. 

The phrase "Efimov physics" was introduced recently in the context 
of atomic gases and Feshbach resonances ^DE05] . but is now used more 
and more frequently. We will adopt this phrase to the extent possible: By 
"Efimov effect" we refer to the infinite accumulation of trimers states at the 
threshold, while by "Efimov physics" we mean the broader concept in cold 
gases using Feshbach resonances and recombinations loss for observing the 
Efimov effect and related features. 

The experiments performed until now strongly indicate that individual 
features of Efimov states have been found. However, the results are yet not 
conclusive |LK J07] . Definite evidence requires two, or even better a whole 
sequence, of the states with the correct relative scaling to be observed. The 
universal window is still quite small and experimental control is not at its 
final level yet. Theoretical physics can help at this point by predicting quan- 
titatively the regions where universality occurs, as well as the corrections at 
the border of these regions. 

The subject of universality in atomic gases will face many new challenges 
during the next decades. In the most recent review on universality in few- 
body systems pBHOGj . three subjects were marked as the frontiers of this 
field: The A'^-body problem for A^ > 4, effective range corrections, and large 
p-wave scattering length. These issues are clearly within reach in ultra-cold 
gases and indicate a promising future. 



1.4 Angle of this Thesis 

In this dissertation we carry out theoretical investigations of universality and 
its limits in few- and many-boson systems. We focus on ultra-cold trapped 
atomic gases, but strive to present the results in universal terms and via 
model-independent parameters. Thus, much of the work may hopefully be 
used or continued in other areas of physics. 

We will investigate effective range corrections and trap effects to Efimov 
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physics, as well as the limits of the universal predictions. Next we investi- 
gate the possibility of an Efimov effect in A^-body systems. We also consider 
trapped BEC-like states from the few-body level, where we will be concerned 
with model-independence and two-body correlations. The approach to an- 
swer these questions will be from the pure few-body level, keeping in mind 
only to include the degrees of freedom needed for physical relevance. Fi- 
nally, we approach Bose-Einstein condensates from the angle of standard 
mean-field theory, but with higher-order interaction terms included. 



Note added. The observation of an Efimov spectrum in an ultra-cold 
gas has just been published (ZDD"'"09j . giving the first definite proof of the 
universal scaling factors. This happened at the very final stage of writing 
this dissertation (I, the author, only became aware of the results a few days 
before handing in the thesis). To keep the chronological order, no contents 
or conclusions have been modified at all. A more elaborate note is found in 
the summary chapter. 
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Theoretical and Numerical 
Background 

In this chapter we consider some of the theoretical background for few-body 
systems. For the two-body systems we discuss basic universal behavior, as 
well as some results on Feshbach resonances, atomic interactions, and model 
potentials. Specifically we discuss the large effective range. For the three- 
body systems we derive the Efimov effect in the hyper-spherical adiabatic ap- 
proximation and explain experimental consequences for recombination rates. 
We also briefly review some mean- field concepts. Finally we describe the 
stochastic variational method which is one of the main numerical tools for 
our work. 

2.1 Two-Body Systems 
2.1.1 Low-Energy Scattering 

Let us briefly review some relevant low-energy concepts for two-body scat- 
tering and bound states. We will be concerned with the limit where where 
only the lowest partial waves, I = 0, contribute |BJ03j . Here the scattering 
can be described by the radial Schrodinger equation 

[-^-Uir)+k^^u{r)=0, (2.1) 

where r = \ri — r2\ is the inter-particle distance, U{r) = 2mV{r)/h? is the 
two-body interaction, m the particle mass, n(r) = R{r)r the s-wave radial 
function, k = \j2mE jU the wave number, and E the relative energy. The 
asymptotic behavior 

u~sin(A;r-h(5(A;)) (2.2) 

defines the s-wave phase shift 5(A;) as function of the incident energy. The 
scattering amplitude /(/c), defined as the ratio of outgoing spherical waves 
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and ingoing plane waves, is in this case spherically symmetric and related 
to the phase shift by 

^^^^ = kcotSik)-tk- ^^-^^ 

The total cross-section is then given by cr = 87r|/(A;)p for identical bosons. 
For small energies one traditionally employ the effective range expansion for 
the phase shift, 

kcot6(k) = -- + —k'^ + ..., (2.4) 
a 2 

where a is the scattering length]^ and Re is the effective range. 

The poles of the scattering amplitude f{k) defines the bound states, 
virtual states, and resonances for the two-body problem. The bound states 
are located on the positive imaginary axis in the complex A:-plane, i.e. k = 
in, K > 0. The binding wave number k is then determined by the equation, 

iKcot 6{iK) + K = 0. (2-5) 

The corresponding energy is given bjj^ 

Ed = = , 2.6 

m m 

where the binding energy is Bd = \E£)\. Contrary to the bound states, the 
virtual states (or anti-bound states) are located on the negative imaginary 
A;-axis, i.e. k = ik, k < 0. The corresponding energy is also real and 
negative, 

Ev = — = (2.7) 

m m 

however, the pole is located on the second Riemann sheet, and can thus only 
be measured indirectly via analytic continuation of the cross section. 

2.1.2 Universal Bound States 

Using only the lowest scattering length term in the effective range expansion 



of eq. (2.4 ) and inserting in eq. ( 2.5 ) , we find a single pole with wave number 



K = 1/a. For a > this gives a bound state with energy 



Ed = - A,, (2.8) 



2 ■ 



which holds when the scattering length is much larger than the inter-atomic 
potential range, tq. In atomic physics such a weakly bound two-body state 
is often referred to as a shallow dimer. This type of state is universal. 



^We use this sign convention for a tliroughout this dissertation, i.e. a weakly bound 
state corresponds to a > 0. 

■^the subscript D referring to "dimer", i.e. a two-body bound state. 
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or model-independent, in the sense that it depends only on the scattering 
length and not the details of the short-range potential: The major part of the 
wave- function is an exponential tail located outside tq, which is an essential 
example of universality. The state can also be described in a zero-range 
model where the usual boundary condition n(0) = is replaced by 



1 du 
u dr 



=0 



/ccot(5(A;)|^_^o 



1 

a 



(2.9) 



1/a and exponentially decreasing wave- 



This leads to the wave number k 
function u ~ exp(— r/a). 

If we also include the effective range in eq. ( 2.4 ) , the binding wave number 
is determined hy R^k^ /2 — k + 1 / a 



1 



1 



= 0. This gives 
Vl - 2Re/a 



(2.10) 



where we have chosen the negative root. The state corresponding to the 
positive root is unphysical, since it is strongly bound with energy E ~ 
h?/{mR'^) comparable to the potential depth |BH06] . Thus for a > we 
have a bound dimer with energy 



Ed = ^ (l- Vl-2i?e/a)' . (2.11) 

The virtual state energy Ey is given by the same expression for a < 0. 



Equation (2.11) holds when the higher-order shape terms in the effective 



range expansion are negligible. Expanding it to lowest order in Re /a gives 



Ed 



1 + 



R. 



+ 0(%) 



(2.12) 



which is then a small correction to eq. (2. 



2.1.3 Feshbach Resonances 

Feshbach resonances is one of the most important experimental tools for 
ultra-cold atomic gases, see |CGJT09] for a recent review. These resonances 
allows for arbitrary tuning of scattering lengths, in particular the s-wave 
scattering length which is relevant for our purpose. This has been a major 
key for controlling macroscopic stability and decay of condensates |DGPS99] , 
and for creating shallow molecular dimers |CGJT09j . 

In the zero-range approximation, the scattering length of a Feshbach res- 
onance can be described by the phenomenological magnetic field dependent 
expression 

a{B) = atg (l - -^^\-) . (2.13) 



9 



2.1. TWO-BODY SYSTEMS 



Here abg is the background scattering length away from resonance, Bq is 
the position of the Feshbach resonance and AB the width. Experimentally, 
these resonance parameters are usually determined indirectly via peaks in 
the loss rate of condensates. 



B„ 



bJO 
1=1 











AB 







abg 



magnetic field B 



Figure 2.1: Scattering length a{B) for the Feshbach resonance model eq. (2.131 as 
function of magnetic field B. The resonance is centered around So with width AS 
and has background scattering length 0^^. 



The phenomenological behavior in eq. (2.13) can be reproduced in vari- 
ous models, a common example is the two-channel model described in |PS02j . 
In this model the resonance width /S.B turns out to be proportional to the 
matrix element connecting the open and closed channels that overlap in 



energy and cause the resonant behavior, see fig. 2.2 



resonant coupling 



closed channel 




open channel 







single channel resonance 














bO 
CD 



















interatomic distance, r 



interatomic distance, r 



Figure 2.2: Left: Schematic representation of a Feshbach resonance. A bound 
state of the closed channel has energy close to the (zero) incoming energy of the 
open channel. The relative position of the channels can be changed by an external 
magnetic field, thereby tuning the coupling. Right: A schematic (shape) resonance 
of a single-channel model. The Feshbach resonance can be modeled by a shape 
resonance if the elastic (open-open channel) scattering is the only contribution. 
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The scattering length, eq. (2.13), is obtained when a two-channel model 
is reduced to an effective single-channel model. The related effective range 
of the effective single-channel zero-range model is determined from the res- 
onance width as (see e.g. jBJKOS] ). 

mA^ABabg ' ^ ^ 

where A^u is the difference between the magnetic moments in the open and 
closed channel. This result holds on- resonance, i.e. near \a\ = oo. In terms 
of common experimental units the effective range reads 

\ mu J \ao J \tJ-B J ^ G ) 

where oq is the Bohr radius, niu is the unified atomic mass unit, the Bohr 
magneton, and /S.B is measured in Gauss. The effective range is always neg- 
ative but can be arbitrarily large for narrow resonances. This is consistent 
with our observations for the zero-range limit of finite-range potentials, see 
below. 



The effective range, eq. (2.14), is constant (independent of magnetic field 
strength, B) and only holds near \a\ = oo. In chapter [7] we derive the more 
general S-dependent effective range, 

2.1.4 Atomic Interactions and Model Potentials 
van der Waals Interaction 

The real atomic potential between neutral atoms can in the Born-Oppenhei- 
mer approximation mainly be described by a short-range and a long-range 
part. The short-range repulsive core arises from the overlapping electron 
clouds, while the attractive long-range van der Waals interaction is due to 
the polarizability of the electron clouds. Asymptotically this van der Waals 
tail goes as 

+ (2.17) 

For ultra-cold gases, where the relative momentum between the atoms is 
small, this is the dominating part of the interaction. The Cg coefficient 
defines the van der Waals length 

I'vdW — I I , (,2.1»j 
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which is the typical interaction length scale between neutral atoms. Ex- 
amples of l^dw are 44.93ao for 23Na and 82.58ao for ^'^Rh |CGJT09] . The 
effective range for the single-channel van der Waals interaction can be esti- 
mated to Re ~ l.Sdlydw for o ^ IvdW 



Model Potentials 

Low energy scattering can be described by model potentials of finite range. 
Typical choices are the finite square well or a van der Waals tail with a 
hard core which can be solves analytically. For example, the square well has 
scattering length and effective range 



1 



tan s 



Re 



1 



a 



(2.19) 



in terms of the finite range ro and dimensionless depth s = ro\/mVo/fi- > 0. 
Each time a bound state is at the threshold, a diverges (s = tt/2 + nvr) and 
the effective range equals the finite range Re = rQ. When a ^ the effective 
range diverges as — 1/a^. The scattering length and effective range for the 
van der Waals tail (— Cq/t^ for r > ro) with hard core (cxd for r < ro) is 
found in |BH06] . 

Another common choice is the Gaussian potential 



V{r) = VoeM-rVrl), 



(2.20) 



which we will use later for numerical calculations. The scattering length 
have features similar to the square well, but the effective range is always 



positive for the attractive case, see fig. 2.3 
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Figure 2.3: Scattering length and effective range for the Gaussian interaction, 



eq. (2.20), in units of the range tq. The first bound state occurs at Vb = 
-2.68lV(mr^). 

Feshbach resonances can also be modeled by single-channel potentials. 
This is of course only possible if the energy is low and elastic scattering is 
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considered. The model potential should have the same low-energy phase 
shift as the real Feshbach resonance. The diverging scattering length can be 
be modeled by any potential, just by tuning a two-body bound state to zero 



energy. However, as we saw in eq. (2.16) the effective range for a Feshbach 



resonance can be very large and negative. Only certain single-channel model 
potentials can reproduce this. 

In general the effective range of a finite-range potential is given by 

/•oo 

Re = 2 [v^oir) - ulir)] dr, (2.21) 
Jo 

where uq is the zero-energy wave-function normalized asymptotically as 
uq — > 1 — r/a, and vq = 1 — r/a is the asymptotic solution extended to 
all r. When considering a potential of finite range rg, we have uq = vq in 
the outer region and the integral only runs to tq. Performing the integration 
over Vq and noting that the integral over Uq is positive we get the bound 
iPC97] 

Re < 2ro f 1 - ^ + (2.22) 

In the resonant limit a = oo, we find Re < 2ro, in particular the effective 
range must be negative if we also approach the zero-range limit We note 
that the large negative effective range occurs when the amplitude of u is 
considerably larger inside the potential as compared to the asymptotic value 
(^^0 = 1). 

This means that to reproduce a large negative effective range with a 
single-channel finite-range potential we need to have a (shape) resonance 
around zero energy, see e.g. fig. |2.2[ The simplest way to do this is by having 
an attractive inner region and a repulsive barrier. This could for example 
be the square well with a square barrier, which can be solved analytically 
[JNW06) . We instead consider the softer potential with barrier 

V{r) = D sech2 {x^^ + B exp (-2(x^ - 2)^^ , (2.23) 



which is shown in fig. 2.4, By tuning the depth of the pocket and height of 
the barrier we can make both a and |i?e| much larger than rg. This leads to 
a large amplitude of the zero-energy wave-function uq inside the potential, 
as compared to the asymptotic region. In the limit \a\ oo,Re — > — oo 
all probability is located inside the potential. This means that universal 
features of the two-body problem can be lost for Re ^ —oo even within the 
normal universal limit of large scattering length. We investigate this point 
for the three-body Efimov effect in chapter |4j 

^The "Wigner bound" [Wig55[ IPC97| says that for zero-range interactions k cot 5 is 



monotonically decreasing with energy, i.e. d{k cot 5 (k)) /d{k^) < 0. Using eq. (2.4 1 this 
also leads to -R. < 0. 
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Figure 2.4: Soft two-body potential V with barrier, eq. (2.23). The parameters 
D = -138.27, B = 128.49 (in units h'^/{mrl)), and x = 4.6667 are tuned to give 
large scattering length a = 556.88 and large negative effective range i?e = —142.86. 
The zero-energy wave-function uo{r) (arb. units) has large amplitude inside the 
barrier. 



In conclusion, when we talk about effective range it means higher-order in 
the two-body scattering dynamics, i.e. the k"^ term of the phase shift expan- 
sion eq. (2.4). It can either be related to the finite range of a single-channel 
model potential (such as the van der Waals interaction) or the higher-order 
term of a coupled-channel Feshbach resonance, eq. ( 2.16| ). However, large 
effective range (compared to the interaction range tq) can only occur near 
a Feshbach or shape resonance, and must be negative in this case. 



2.2 Three-Body Systems 

We now describe the three-body problem with the hyper-spherical adia- 
batic approximation. We consider the zero-range limit with a Faddeev-type 
decomposition and appropriate boundary conditions to derive the effective 
hyper-radial potential responsible for the Efimov effect. Applications and 
experiments are outlined. 



2.2.1 Hamiltonian and Coordinates 

Let us first describe the hyper-spherical adiabatic approximations, details 
can be found in |NF JGOlj IF JOlb] . We consider the three-body system with 
the center-of-mass (cm) Hamiltonian 

H = f + Y,y{rij), (2.24) 

where T is the kinetic energy operator in the center-of-mass frame and V 
is the two-body interaction. We denote the vector from particle j to i by 
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fij = fi — fj 1 and the vector from the center-of-mass of the pair (j, k) to 
particle i by = Ti — {vj + rfc)/2. In the equal mass system the Jacobi 

coordinates {xi,y^) for i = 1,2,3 are defined aa^ 



1 



(2.25) 



where {i,j,k} is a cyclic permutation of {1,2,3}, see fig. 2.5 

2 

1 





Figure 2.5: Relative coordinates r,y and one set {i — 1) of Jacobi-coordinates 
The six hyper-spherical coordinates (p, ai,Qxi,^yi) are defined as 



psmai, Ui = p cos ai 



(2.26) 



where p > is the hyper-radius, ai E [0,7r/2] is the hyper-angle, ^Ixi = 
the angles describing the direction of ccj, and similarly Jl^j for i/j 
The hyper-radius 

1 , 



2 I 2 

Xi +yi 
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+ r]^3 -I- r. 



23) 



(2.27) 



is independent of the chosen Jacobi-set and describes the average size of 
the three-body system. It is the only dimension-full coordinate of the 
hyper-spherical coordinates. All the angles are denoted together by J7j = 
{ai,Qxi,^yi)- For fixed p, the hyper-angle Ui describes (together with 
(Qxij^yi)) the internal configuration. For ~ particles j,k are close 



to each other, and for ai ~ 7r/2 particle i lies between j and k, see fig. 2.6 



The volume elements of the Jacobi and hyper-spherical coordinates are 
related by 

dxidyi = p^dpdn, ^^^^^ 



d^} = sin^ cos^ a,- dJ7^,- dO 



yi- 



The kinetic energy operator in hyper-spherical coordinates is given by 




5 d A2 
p dp p^ 



(2.29) 



^Different proportionality factors are used in literature, we use tliis one to be consistent 
with later definitions for > 3. 
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Figure 2.6: Two configurations with the same hyper-radius p but different hyper- 
angles ai. Left: ai ^ 0. Right: ai ^ tt/2. 



where the square of the grand angular momentum operator A is 



sin 2ai daf 



sin(2aj 



4 + 



+ 



/2. 



sm Q,- 



cos 2 a,- 



(2.30) 



Here l^^ and ly^ are the angular momentum operators corresponding to x 



yi 



and We only consider zero total angular momentum, i.e. Ixi 
Two different sets of Jacobi coordinates i and j are connected through the 
so-caUed kinematic rotation |NFJGOl] 



^3 \ _ 



- COS -fij sin7ij 

- sin 7ij - cos 7ij 



(2.31) 



For three identical particles the "rotation angle" is "jij = a{i, j, k}7r/3, where 
(j{i,j,k} is the sign of the permutation {i,j,k}. By calculating ic? from 
eq. (2.31) and using eq. (2.26), the different hyper-angles are then related 
by 



sin^ aj = cos^ -jij sin^ ai + sin^ jij cos^ ai 
— 2 sin cos 7ij sin Oi cos ai cos 9i 



(2.32) 



where 9i £ [0, 2tt] is the angle between Xi and ^/j. Using | cos 6i\ < 1 one can 
show that for a fixed hyper- angle ai in the coordinate set i, the hyper-angle 
in set j is restricted by 



vr 



— ai 



(2.33) 



2.2.2 Hyper-Spherical Adiabatic Expansion 

In the hyper-spherical adiabatic expansion one first solve the angular part 
(O) of the Schrodinger equation for fixed p, 



P + ^p^Y.V{V2psma^) \ $„(/9, 1^) = A„$„(p, O). 

^ i=l / 



(2.34) 
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This gives a complete set of eigenfunctions $n(/0, 0,) and corresponding eigen- 
values \n{p) as functions of p. We expand the total wave function a^ 



(2.35) 



Inserting this into the full Schrodinger equation gives a coupled set of hyper- 
radial equations for the coefficients fn{p), 

+ |?(^efr,n(p) - E)\ Up) = (2^nn'|- + Qnn')fn'{p), (2.36) 



where the effective hyper-radial potentials are 



^eff,n(p) 



Xnjp) + 15/4 

2m 



Qr 



(2.37) 



and the non-adiabatic couplings are 



Qnn'ip) = ( ^>„(p,S^) 



d 



dp 
52 



$„/(/>, 17) 



<^n'{p,^) 



(2.38) 



The brackets denote integration over 0,. The angular wave- functions are 
normalized to unity for fixed p, i.e. {^{p,0,)\^{p,0,))q = 1. The identity 
Pnn = also holds. In the strict adiabatic limit where all the off-diagonal 
coupling terms Pnn' and Qnn' vanish, the hyper-radial equations, eq. (2.36), 
decouple. 



2.2.3 Hyper-Angular Faddeev Decomposition 



Let us first rewrite the angular equation eq. (2.34) as 

3 \ 



+ Y,U{ai) (sin(2a,)^>) = 0, 



(2.39) 



1=1 



^p2y(\/2psin Oj), the eigenvalue 



where the reduced interaction is U (aj) 

Is v"^ = \ + A, and we have omitted the subscripts n. We now split the 
hyper- angular wave- function into three identical Faddeev components. 



sin(2ai) 



(2.40) 



The phase-factor p ^^'^ is simply conventional to simplify equations below. 
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which fulfill the three Faddeev equations 
52 



da 



(2.41) 



Adding these equations leads to eq. (2.39). The components il)i,il)j,ipk have 



the same functional form, but are expressed in different Jacobi sets. If one 
function, say ■0 = V'ii is given we obtain the "rotated" component Tl[tp] = '4)j 
by expressing tp in coordinates j and projecting onto s-waves. Since we have 
identical particles the two rotated components are identical and we may 
write 



sm(2aj) 

Specifically, the rotation operator is given by 



(2.42) 



sin(2aj) 



(47r)"' sm(2aj) 



1 



2tt 



sin(2aj 



(2.43) 



sin 6i dOi , 



where aj depends on 9i via eq. (2.32) (with fixed a,). From eq. (2.32) we 
also find sin(2aj) sin 0j d6'j = (4/-v/3) sin(2aj) daj. The rotation operator 
can then be written as 



7^ 



TT I TT „ I 



V3 



(2.44) 



The new integration limits were obtained from eq. (2.33). The limits are 



also shown in a useful form in fig. 2.7 



2.2.4 Zero-Range Solution 



We now take the zero-range limit, i.e. [/ = in eq. (2.41). The solutions 
with boundary condition ip^ir/l) = are 



V'(aj) = M{p)sm{iy{ai - -)), 



(2.45) 



with normalization M{p) depending on the eigenvalue i'{p). It is straight- 

(2.46) 



forward to rotate this solution with eq. (2.44) and we obtain 
4 I sin(— i/Oj) sin(z^?) 



, < < f 
\ sin(z^(aj - |)) sin(i/|) , f < < f . 



Note that this function is zero at Oj = and 7r/2. It is continuous but not 
differentiable at = vr/3, which comes from the fact that the zero-range 
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Figure 2.7: Integration limits for the rotation operator eq. (2.44). For a given a 



the shaded region determines the interval of integration for aj . A specific value of 
ai is indicated by the thick blue line. 



model has ip{^) 7^ (in contrast to any finite-range potential). However, 
the total wave-function $ is smooth everywhere. 

We also need to impose a boundary condition at = to determine 
the eigenvalue v. In the two-body case the asymptotic radial wave-function 
is given by u{r) = rR{r) oc sin(A;r -|- 5{k)). The zero-range limit can be 
formulated as the free problem with boundary condition 



d\n{rR{r)) 
dr 



r=0 



1 

a 



(2.47) 



For the three-body problem we apply the same boundary condition in each 
of the three Jacobi sets (for fixed p), 



(91n(rjfc$) 



dr 



'jk 



1 

a 



(2.48) 



In the limit rjk — > we have to lowest order rjk = V^Xi ~ \f2pai. Using 

then becomes 

1 aln(V' + 27^[V']) 



eq. (2.42) the logarithmic derivative then becomes 
91n(rjfe^>) 



dr 



'jk 



rjk=0 



\f2p dUi 

to lowest order in Oj. By inserting eq. ( |2.45 ) and eq. (2.46) we obtain 
1 -i/cos(z^f) ^sin(z^f) 



(2.49) 



V2p 



sin(z^f) 



1 

a 



(2.50) 
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which determines the angular eigenvalues Vnip)- This equation was first 
obtained by V. Efimov [EfiTTj . 



2.2.5 The Efimov Eff"ect, \a\ = oo 



Let us consider the resonant limit \a\ = oo. Then eq. (2.50) has a single 
imaginary solutioij^ fo — 1.00624i determined by 

TT 8 TT 

-iycos{u-) + -^sm{u-)=0. (2.51) 

We will refer to vq or the related effective potential as the Efimov 

solution]^ The related effective hyper-radial potential becomes 



V.^fl(p) = - — ^—^. (2.52) 



where ^ = \ vq\ = 1.60024. Since the angular solutions only depend on p via 
the constant eigenvalues v^, all non-adiabatic couplings Pn,n'^ Qn,n' vanish 
and the adiabatic approximation becomes exact in this limit. The lowest 
hyper-radial equation becomes 

-^-^^J/o(p) = 0, (2.53) 
where k = ^2m{—E')l%. Scaling the variables as /o = fo/^/p, p = up 



eq. (2.53) turns into a Bessel equation. The solution is the modified Bessel 
function of the second kind (of imaginary order) fo{p) = Ki^{p) |AS95| 
sec. 9.6.1]. Thus the hyper-radial solutions with binding wave number 
become 



fPip) = ^K,i:{K(^>p) ^p\T{iO\sm ^ein(K^"V) + 0) ■ (2-54) 

The last approximation is for the low-distance limit Kp <ti 1 where we used 
several identities and expansions for the T function. The phase is given 
by 6* = arg{r(-i^)2-*«} + Tr/2 ~ 0.8747r, and \T{i^)\ ~ 0.5148. 

We see that the solutions have infinitely many nodes at short distances, 
meaning that there are infinitely many low-lying states. This is the well 
known Thomas effect in the three-body problem with zero-range interac- 
tions and fixed two-body binding energy (or scattering length a) |Tho35| 
INFJ(;01] . 

The Thomas collapse is unphysical since the real short-distance details 
must be taken into account. Thus, we need to introduce a regularization. 



^The equation also have real (continuum) solutions i/rn ~ 2m, where m is a positive 
integer, except m = 2 which is spurious. 



7 



Although it should also be attributed to G. S. Danilov [Dan61| . see [EfiTl] . 
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e.g. by setting /o(po) = at some short arbitrary distance pQ. Using 
eq. (2.54) this leads to k*^"^ = exp(— riTr — 0)/po, giving the characteristic 



scaling between the energies, 



= e-""/« ~ 22.7- 



-n 



or :^=e-2™/? 



515.0" 



(2.55) 



This is the essence of the Efimov effect, where infinitely many loosely bound 
states (Efimov trimers) accumulate a zero energy. Three hyper-radial wave- 
functions are shown in fig. |2.8[ The root-mean-square hyper-radii of the 
Efimov states scale accordingly as 



{P 



2\l/2 



/o 



22.7". 



(2.56) 
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Figure 2.8: Hyper-radial wave-functions, eq. (2.541, for the Efimov states in the 
effective 1/p^ potential eq. (2.521. The scale kq is the binding wave number for 



state rt = 0. The solutions have infinitely many nodes as p — > 0, see inset. 



2.2.6 Universal Scaling, |a| < cxd 



Let us finally discuss the case of finite but large scattering length. The 
constant solution uq in eq. (2.50) still holds, but now only for p < |a|, thus 
the effective l/p^ potential will be modified at p > \a\. This means that 
all the Efimov states with hyper-radii larger than \a\ (see e.g fig. 2.8) will 
be affected. However, Efimov states with hyper-radii smaller than \a\ only 
have exponential tails in the affected region and are therefore essentially 
unchanged. Note also, that the regularization can be considered fixed: The 



21 



2.2. THREE-BODY SYSTEMS 




-3-2-101234 



l/a2 

Figure 2.9: Discrete universal scaling of Efimov states. Trimer energy levels, E'^"-' 
are shown as functions of scattering length, a. Points indicate the critical scattering 
lengths where the trimer crosses the trimer and atom-dimer thresholds. Both axes 
were scaled by a power of 1/8 to reduce the large scaling factor 515.0 to 2.18, i.e. the 
axes actually show [fi^ / {171,0^ Eoc)]^^^ and [Et/Eoc]^^^ , with £'00 being an arbitrary 
regularization scale modulo 515.0 (the energy of a trimer at a = 00). 



infinite scattering length is tuned by adjusting the two-body energy infinites- 
imally around zero leaving the underlying finite-range potentials fixed. 

(n) 

The resulting behavior of the Efimov trimer energies, Ej. , are shown 
in fig. 2.9 as function of inverse squared scattering length l/o^. On the 
a > side the two-body system has a shallow bound state (dimer) with 
energy = —h?/ma'^, eq. (2.8) (indicated by the blue diagonal line). For 
the three-body case this defines the atom-dimer fragmentation threshold. 
For a < the two-body system has a shallow virtual state Ey and the 
three-body system has a three-body fragmentation threshold at zero energy 
(indicated by the red horizontal line). At \a\ =00 the infinite sequence of 
Efimov trimer s is shown. 

When a is decreased from infinity to a < the effective hyper-radial 
potential moves up at large distances p > \a\, and the most loosely bound 
Efimov trimers are pushed into the continuum. For a > the effective poten- 
tial decreases, now converging to the dimer binding energy Ves — > E^ < 
for larger hyper-radii. Each Efimov trimer moves down correspondingly, 
until it disappear into the atom-dimer continuum. 

Because of the scale-invariant l//0^ potential all the Efimov trimer en- 
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ergies scale the same way, and only a single universal function needs to be 
calculated. The specific form can be found by solving the Schrodinger or 
Faddeev equations numerically for short-ranged potentials (as done in chap- 
ter Isl), or by using effective field theory (EFT) |BH03l IBHK031 IBH06| . In 

(n) 

the exact universal limit the energy Ej^ of the n'th trimer state can be 
parametrized as 

4") - ^ = e^(^)/«e-2WCi7;^ , (2.57) 

where = 1.00624, A(ip) is a universal function |BH06j of the angle ip E 
[— tt; — 7r/4] given by tan'^ ip = —E^^mo? /ti^ (i.e. Lp is the polar angle in the 



coordinates of fig. 2.9). E^o is the regularization scale (here the energy of 
state n = at a = oo) which cannot be determined within the zero-range 
theory. It has to be supplied from phenomenological information. 

The two critical scattering lengths and where the Efimov trimer 



2.9 



(n) 

Ej. disappears into the three-body and atom-dimer continuum, see fig. 
are given explicitly in the zero-range limit by |BH061 IGME08] 

a^") = -1.50763/k("\ = 0.0707645/k('^), (2.58) 



where k*-"^ = y —mE^\a = oo)/h is the wave number for the n'th trimer 
energy at infinite scattering length. 



Three-Body Recombination in Experiments 

The critical scattering lengths and can be observed via maxima and 
minima in the atom loss rate of cold gases [KMW+06]. In the ultra-cold limit 
the loss due to two-body collisions B + B ^ B2 are suppressed by energy- 
momentum conservation. Therefore three-body collisions B + B + B — > 
B + B2 can be considered as the dominant mechanism for the total loss rate. 
The general behavior for the three-body recombination rate is \a\^ for large 
a |NM99| IEGB99I IBH06] . On top of this background additional resonance 
or interference effects effects from the Efimov trimers should be seen. 

The first case where the Efimov trimer disappear into the three-body 
continuum (a < 0) leads to a resonant peak in the recombination rate. In 
numerical hyper-spherical descriptions this is explained by a tunneling into 
a small /o-region and subsequent decay into strongly bound molecular states 
|EG06j . The second case where the trimer crosses the atom-dimer threshold 
(a > 0) leads to a minima in the recombination ratej^ This can be explained 
as two interfering pathways from the incoming hyper-spherical channel to 
the weakly bound atom-dimer channel |EG06j . 

*These minima are actually located between the positions where the Efimov trimer hit 
the threshold. 
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Until now, only a few experiments measuring such effect have been per- 
formed, namely in ^^^Cs gases [K MW+061 |KFM+09| and also in a three- 



component Li Fermi gas [WLO^ 



The experiments will be conclusive 



when the scaling factor 22.7 is observed (or a correspondingly smaller factor 
for unequal mass systems). 



2.3 Mean-Field Condensates 

The Gross-Pitaevskii (GP) equation is based on a mean-field ansatz with 
the condensate wave-function (or order parameter) ^'(r) together with the 
zero-range contact interaction 

VzR{r) = U^5{T), Uo = . (2.59) 

m 

The coupling constant Uq is chosen such that the correct energy shift is 
reproduced, see chapterjTj This is equivalent to a Born approximation where 
o-born is replaced by the physical scattering length a. This leads to the GP 
energy functional 

E{^) = I dr (I^IV^P + Fext(r)|^P + ^l^r) , (2.60) 

which includes kinetic energy, the external trapping potential Vext, and the 



interaction term from eq. (2.59). A variation of this functional with the 



constraint of fixed particle number N leads to the GP equation 

|^V2 + Vextir) + Uo\^{r)\ - fi) ^(r) = 0, (2.61) 
2m J 

where // is the chemical potential. 

The GP equation (or variations and extensions of it) has been able to 
describe a wide variety of BEG phenomena, e.g. the spatial condensate pro- 
file in different traps, solitons in homogeneous condensates, rotating BEGs 
with vortex profiles, stability of low-energy modes etc. |PS02| IPS03] . A 
specific useful approach is the Thomas-Fermi (TF) approximation where 
the kinetic energy term is neglected and analytical solutions are possible. 
The approximation holds for repulsive gases (a > 0) with many particles 
(e.g. Na/ht ^ 1 in the harmonic trap of size ht = y/h/{muj)). See chapter [s] 
for elaborate discussions. 



From the perspective that the contact interaction eq. (2.59) should pro- 
duce correct energy shifts, it is only the lowest order approximation. In chap- 
ter [7] we derive the higher-order interaction term proportional to d{r)\/^ -\- 
V^5(r), and the related modified GP equation with an extra term. The 
purpose of chapter [7] is to investigate possible effects of such higher-order 
terms via variational and numerical solutions. Ghapter [7] approaches the 
same question analytically in the TF limit. 
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2.4 Stochastic Variational Method 

The stochastic variational method (SVM) |SV98| ISorOS) is a numerical min- 
imization technique to find approximate solutions of few- and many-body 
problems. The method is based on two basic ingredients, i) stochastic ran- 
dom sampling and ii) deterministic linear variation. 

One of the major advantages of the SVM is that one can calculate a large 
number of bound states, obtaining both the spectrum and wave-functions 
directly. Another feature is that one can control the amount of inter-particle 
correlations allowed in the calculated states - depending on the application 
one may e.g. look for highly correlated few-body states or weakly correlated 
many-body states. The method is also numerically efficient and it can be 
parallelized with close to linear scalability. 



2.4.1 The Linear Variational Principle 

Let us consider the stationary many-body Schrodinger equation = 
and denote the eigenfunctions and eigenvalues by and En, respectively. 
The well-known Rayleigh-Ritz variational principle states that the varia- 
tional energy £ evaluated with an arbitrary trial-function ^ gives an upper 
bound to the exact ground state energy of the Hamiltonian H, i.e. 

£ = ™ > (2.62) 

The approach below extends this principle to any number of excited states. 

To find approximate variational solutions to the full Schrodinger equa- 
tion we first restrict the problem to a smaller space spanned by a (possibly 
over-complete) set of K basis functions ipk- Let us first express the wave- 
function in this space, i.e. 

K 

^ = Y.Cki^k, (2.63) 

k=i 

where are expansion coefficients. These coefficients determine ^ com- 
pletely within the space {ipk}, although the expansion may not be unique. 



By inserting eq. (2.63) in the Schrodinger equation and projecting onto the 



j'th basis state ipj we obtain 

K K 



HjkCk = £'^ SjkCk- (2.64) 



fc=i fc=i 



Here the Hamiltonian matrix elements and the overlap matrix in the basis 
{ipk} are given by 

Hjk = {Tpj\H\^pk), Sjk = {il^jlipk), (2.65) 
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where the brackets mean integration over all coordinates in ■0^. In terms of 
matrix multiplication, eq. (2.65 ) reads He = 8Sc, where c = (ci, . . . , ck)'^ 



and H, S have the matrix elements Hjk, Sjk- This is a generalized eigenvalue 
equation of size K with K real eigenvalues £i < • • • < £k and corresponding 
eigenfunctions ci , . . . , cj^ . 

The above procedure shows how to find approximate solutions to the full 
Schrodinger equation within a smaller subspace. Let us now specify what 
we mean by the term "approximate". We take the space {ipk} to be fixed, 



and treat the linear coefficients Ck in eq. (2.63) as variational parameters. 



The variational energy, eq. (2.62 ), expressed in the basis {^fe} then becomes 



2^ He 



(2.66) 



where = (c^, . . . , c|^). The stationary solutions within our subspace cor- 
respond to d£/dci = d£/dc* = for all i. We multiply eq. (2.66) by c^Sc 
and take the partial derivative d/dc* (with fixed and Cj,c*, j 7^ i). This 
leads directly to the generalized eigenvalue equation 



He = £Se. 



(2.67) 



Thus, the energies £ found by projecting the Schrodinger equation onto the 
fixed subspace {^fc} are actually the best possible energies as seen from the 
variational perspective. 

Furthermore, it can be shown rigorously jSV98| that the generalized 



(ordered) eigenvalues £1 < 
(ordered) eigenvalues Ei < 



< £k are strict upper bounds to the exact 
^ Ek < • • • , namely 



Ei<£i, ••• , Ek<£k- (2.68) 

It can also be shown |SV98] that by adding another basis state ipK+i, while 
keeping the first ipi, . . . ^if^K fixed, the upper bounds in eq. (2.68) can only 



get better. This linear variational principle extends the Rayleigh-Ritz vari- 
ational principle to an arbitrary number of bound states. 



2.4.2 Basis States and Minimization Procedure 

To employ the variational principle above we must generate a set of basis 
states V'fc- These states should form a complete set, or at least be able to 
describe the physical degrees of freedoms of interest. For the approach to be 
numerically tractable one should also be able to evaluate the matrix elements 
Hij and Sij analytically. Common choices |SV98| ISorOS] are exponential 
functions or Gaussians (see below). Each basis state ip^ is then parametrized 
by one or more nonlinear parameters, which we denote by the single symbol 
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There are many different minimization strategies for the SVM. The sim- 
plest one is to randomly pick a large number of functions (j)^ (i.e. nonlinear 



parameters {afc}) and solve (2.67). This approach has two minor disad- 
vantages, i) many of the random basis states (p^ may not be necessary to 
describe the real eigenstates, and ii) solving the full generalized eigenvalue 
equation of large dimension is numerically very inefficient. Luckily both 
issues can be solved by adding or modifying only a single state at a time. 

Let us assume that the generalized eigenvalue equation has been solved 
in the X-dimensional space spanned by "01 ; • • • ; with resulting eigenval- 
ues £^1 , . . . , £k and eigenstates (f)i, . . . , (f)K ■ We now expand the space by a 
single basis state tpx+i- Using the Gram-Schmidt process we construct (pK+i 
from ipK+i which is orthogonal to all other (j)^s. In the basis (pi, ... , (pK+i 



the generalized eigenvalue problem (2.67) is almost diagonal and reduces to 



finding the roots of the simple function 

D{£') = J2j^-£'-hK^i^ (2.69) 
k=i ^ 

where hk = {(j)k\H\(l)K+i) ■ D(£') has exactly K + 1 roots <?i, . . . ,£k+i which 
are the new eigenvalues. Thus, by adding only a single basis state at a time, 
the generalized eigenvalue problem reduces to one-dimensional root-finding. 
Let us give a short example on how the minimization procedure can be 



done within the SVM, as shown schematically in fig. 2.10, We start by 
adding one basis state at a time, minimizing the ground state. For the /c'th 
basis state ipk we stochastically generate a set {/?m} of candidates for the 
parameters {a^^^}. For each candidate we solve the generalized eigenvalue 



equation eq. (2.67) from eq. (2.69). We choose the Pm giving the lowest 
ground state £i. All other eigenvalues are also guaranteed to be improved. 
After all K eigenstates are found one may refine the basis: Redundant basis 
states can be replaced by new ones still optimizing £i^ The entire procedure 
can then be repeated for excited states. 

This example is quite simple compared to the actual SVM implementa- 
tion we use |Sor05] . but it contains the essential details. The method can be 
tuned in many ways, e.g. continuous optimization of the random sampling 
algorithm and selection of appropriate basis states |Sor05j . 

2.4.3 Center-of-Mass and Symmetrization 

The specific many-body Hamiltonian used with the SVM is given by 

1=1 * 1=1 i<j 



^This is at the expense of the excited states which almost certainly go up. 
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Generate random set 
{/?„J for basis ^fe({a«}) 



Solve generalized eigenvalue 
problem for each /3m 



Choose the /3m giving lowest 
energy £i . Set a^*^ = /3m 



Output: £i, . . . ,£k 



Possible refining: Cycle through 
{a'*'}, replacing redundant 
basis states with new ones. 



Optimization of excited 
state n: Add new basis 
states k = K + 1, . . . K' , 
minimizing instead of £\. 



Figure 2.10: A simple example of the SVM minimization procedure. The ground 
state is found first, followed by a possible refining and minimization of excited 
states. Since this strategy only changes one basis state at a time, the generalized 
eigenvalue problem eq. (2.67) reduces to a faster root-finding problem. 



with a harmonic oscillator trap with angular frequency u and trap length 
ht = -s/h/ (muj), and two-body interactions V{rij) chosen either as the zero- 
range or Gaussian model potentials, eqs. (2.59) and ( |2.20 ). Since the con- 



fining trap is harmonic, the center-of-mass (cm) R = "^iLifi/N can be 
separated out. We will only consider the center-of-mass ground state mo- 
tion, i.e. 

NR'^ 

M/„n(i?)=exp(-^), (2.71) 

with corresponding energy Ecm = \hoj. 

Since we consider identical bosons the total wave-function must be com- 
pletely symmetric with respect to interchange of every pair. This is accom- 
plished with the symmetrization operator 

^ ' a 

Here Pa is the permutation operator corresponding to the permutation a, 
and the sum is over all possible permutations. The symmetrization does 
not affect the conclusions of subsection 12.4.11 and 12.4.21 since it is a linear 
operation. The specific wave-function is then represented as 

K 

^{n, ...,rN) = ^cm(-R) 5^ Cfc V'fc({a!f }, {Tij}). (2.73) 

k=l 
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The basis functions, tpk, depend on all the internal coordinates {rij}, and 
are parametrized by a, possibly large, number of variables {a^j^}- 



2.4.4 Correlated Gaussians 



As the specific basis functions ■i/'A: 
Gaussians, 



we use the so-called explicitly correlated 



^fuU 



k=l \ i<j J 



(2.74) 



This will also be referred to as the full correlated basis. The basis is complete 
for zero total angular momentum states and it allows all types of clustering in 
the system. This basis choice allows analytical evaluation of matrix elements 
for the Gaussian interaction, see [ SV981 ISorOS] . However, because of the 
symmetrization the computational complexity is of order 0{N\), and the 
method is only possible for relatively small number of particles, say ^ 5. 



The zero-range potential, eq. (2.59), requires an uncorrelated wave-func- 



tion which we choose in the form of the linear combination of the hyper- 
radial basis-functions 



K 

Ck exp 



fe=i 



where p is the hyper-radius for particles. 



N 



N 



N 

E 

i=l 



N 



NW 



(2.75) 



(2.76) 



generalizing the A = 3 definition, eq. (2.27). This function is totally sym 



metric and thus does not require the symmetrization operator S. It is a spe- 



cific choice for the parameters of eq. (2.74) where the different pairs of par- 
ticles have the same parameters a^^'' = oiij \ see eq. (2.76). This is reminis- 



cent of a mean-field approximation, since all particles are treated identically. 
Note that only one basis function K = 1 with a*-"*^^ = corresponds 
to the single-particle product of the non-interacting case. The zero-range 



potential with the hyper-radial variational wave-function eq. (2.75) provides 
results similar to the Gross-Pitaevskii equation |SFJ05] . 

For a typical system of trapped atoms even when the scattering length is 
large the density of the system remains small, nr^ <^ 1, and one can assume 
that only binary collisions play a significant role in the system dynamics. 
In this approximation the variational wave-function can be simplified by 
allowing only two-body correlations in the basis-functions. 



K 



2B = ^cm S'^Ck exp 



k=l 



12 



(2.77) 
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where a^*^-* and P^^'^ are the nonhnear parameters. This form is equiva- 
lent to the Faddeev-Yakubovski expansion \I' = F{p) Yli<j 4'{Pi '''ij) used in 
|SSJF05al [SS JF05b] where each of the two-body amphtudes only depend on 
the distance between two-particles. The symmetrization of this function can 
be done analytically |Sor05| ISFJ 05] by collecting similar terms which greatly 
simplifies the numerical calculations. The two-body correlated approach is 
therefore very useful in dilute many-body systems. 

In the trap of size bt with interactions of range ro, the non- linear param- 
eters a^j^ are typically optimized stochastically in an interval from 1/6^ to 
1/rQ, or even larger. This allows both short and long-range correlations in 
the system. 
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Efimov Physics: Finite 
Range and Trap Effects 



3.1 Introduction 

Much effort lias been devoted to extract universal features in many branches 
of physics, because the applications become more transparent. Several top- 
ics concerned with Efimov physics are of special interest here: i) the Efi- 
mov effect |Efi70| |Efi90j where anomalous three-body properties arise at 
the threshold of binding of two particles, ii) halos and Borromean systems 
|JRFG04] where required scaling properties are equivalent to large probabil- 
ities in non-classical regions, iii) universality for few-body systems derived 
from a zero-range interaction |AFT971 IFTDA991 IBH06j . iv) general prop- 
erties of Bose-Einstein condensates |PS02] . and v) their instability due to 
three-body recombination KMW"'"06] . 



In all these five cases the two-body s-wave scattering length a is the cru- 
cial and only characterizing parameter which is independent of the details of 
the responsible potentials. In fact the same scattering length can be achieved 
by disparate potentials. It is highly desirable to assess the uncertainties in 
the results from the leading order terms and extend to include correction 
terms. An extension inevitably needs more details which again should be 
expressed in terms of model-independent parameters. The obvious choice is 
then to exploit the effective range expansion of the two-body phase shifts 
where the leading term containing a is given by the zero-energy limit, and 
the second term is proportional to the energy and the effective range Rg. 
Inclusion of even higher-order terms is usually not productive because it is 
either inefficient or difficult, conceptually and practically. 

The purpose of this chapter is to go beyond the scattering length ap- 
proximation for Efimov physics. We express the corrections in terms of the 
model-independent, low-energy scattering parameters, scattering length and 
effective range. For negative scattering lengths an attempt in this direction 
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was made by varying the form of the potentials and computing the critical 
strengths for few-body binding |RF94l |MFK+00| . This was aimed at find- 
ing the Borromean window where three particles are bound even though all 
the two-body systems are unbound. The results in |MFK"'"00 are however 



not expressed in terms of model-independent parameters. Other attempts 
to include effective range terms within effective field theory are reviewed in 
|BHn6] . 

Specifically, we calculate the low-lying energy spectrum for three trapped 
identical bosons interacting via finite-range two-body potentials. The spec- 
trum is compared to results from the zero-range model. The thresholds for 
trimer binding and atom-dimer binding are extracted and related to effective 
range corrections. Effective range corrections to Efimov physics and Bor- 
romean binding are two aspects of the same effect, and we connect these two 
regions qualitatively. The Borromean window becomes slightly narrower for 
substantial effective ranges. The structure at the atom-dimer threshold is 
an atom far away from the dimer and the major energy correction is due 
to the change of dimer energy with effective range. This structure becomes 
less pronounced when the effective range increases. Comparisons with recent 
results from effective field theory are carried out. 



3.2 Procedure 



We consider = 3 identical bosons with mass m and coordinates in 
a spherical harmonic trap with frequency uj and corresponding trap length 
bt = \/hJrnu). As the two-body interaction we mainly use an attractive 

We also include results for 
The strength 



Gaussian, eq. (2.20), of fixed range tq, see fig. 2.3 



many other potential shapes based on input from |MFK^00 



Vq is varied within the interval where the potential either cannot support 
bound states (scattering length a < 0) or have only one bound state (a > 0). 
The positive effective range is then also a given function of Vq. However, 
near the resonance, |a| = oo, i?e varies slowly as function of Vq. We choose 
the trap to be much larger than the potential range, 6t/ro = 3965. 

The wave functions and energies are found with the stochastic varia- 
tional method described in chapter |2] using the fully correlated Gaussians, 
eq. (2.74). The non-linear parameters [ajj^] ""^^^ in eq. (2.74) are optimized 
stochastically in a large interval covering values of order tq up to values of 
order ht- This allows both short and long-range correlations. The separated 
center-of-mass motion is given by the lowest oscillator wave function with 
the energy Ecm = 
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3.3 Zero-Range Approximation 

Let us briefly refresh the results of the zero-range model in chapter [2} In the 
zero-range limit of the two-body potential, the scattering length a is the only 
remaining interaction parameter. For a > the weakly bound dimer has an 
energy Ed = —h'^/{ma'^), eq. (2.8), while the virtual state Ey is given by 



the same expression for a < 0. The correction to Ed and Ey in the effective 



range expansion is given by eq. (2.11). In all numerical applications in this 



chapter we have found agreement with eq. (2.11) and higher-order terms on 
the two-body level are not needed. 

For three particles many bound Efimov trimers may exist even when the 
two-body system cannot support any bound states or is weakly bound. The 



energies Ej^ of these states are given in the zero- range limit by eq. (2.57) 



see also fig. 2.9 The energies scale geometrically with a factor of 515.0, 
while the hyper-radii scale with 22.7. The energies are only determined up 
to a three-body regularization scale, E^o = El^\a = oo), which cannot be 
determined within the zero-range theory. It must be fixed by real finite- 
range calculations or experimental data. 

For a < the Efimov trimers cross the threshold for three-body binding, 
while in the opposite direction of a > the thresholds are crossed for binding 
of the atom-dimer system. The two corresponding critical scattering lengths 
a^"^ and oj'^ corresponding to = and E^ = Ed, are given in the 



lowest-order zero-range limit by eq. (2.58), 



a^") = -1.50763/k("), 4"^ = 0.0707645/k('^). (3.1) 

Here k^") = {—•mE^\a = cxd))^/^//i is the wave number for the n'th trimer 

energy at infinite scattering length. The critical scattering lengths {cl'z^) 
can be observed via maxima (minima) in the three-body recombination rate 
of cold atomic gases jKMW+ 06] as described in chapter [l] and |2] Since, in 
this lowest order model, the wave number is given by k = 1/a, eq. ( |3.1[ ) can 
instead be written as 

1/4"^ = -1.50763/k("\ 1/4"^ = 0.0707645/k("\ (3.2) 
where and the respective critical wave numbers. 

3.4 Results 

3.4.1 Three-Body Energies: Overview 

We compute the energies and wave-functions corresponding to the Hamil- 



tonian in eq. (2.70) with a Gaussian interaction by variation of the form 



in eq. (2.74). The resulting three-body energies are shown in fig. 3.1 as 
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function of scattering length. The center-of-mass energy is subtracted and 
the axes have been scaled by a power of 1/8 in order to obtain a well pro- 
portioned figure. The critical scattering lengths and effective ranges for our 



finite-range calculations are denoted, in analogy to (|3.1|), with a subscript 



I.e. a 



a 



W^' and R 



in) 



3 



I 



bt 



ro 



3 



I 




-4-2 2 

sgn(a)|62/a2|i/8 

(n) 

Figure 3.1: Three-body energy levels, E)p measured relative to the center-of-mass 
energy Ecm, £is function of scattering length, a. Finite- and zero-range results 
are shown as points and curves, respectively. As finite-range interaction we use 
a Gaussian with adjustable depth. The upper and lower parts, corresponding to 
positive and negative energies, are shown separately. The central value of \a\ — oo 
corresponds to the threshold for binding two particles. The smallest a (far left and 
far right) corresponds to the interaction range, \a\ = rp. There is a quasi-continuum 
for a > (shaded) which shows the threshold for atom-dimer binding (individual 
states are not shown). 

The three-body energies obtained in the zero-range model are charac- 
terized by the arbitrary regularization scale i?oo- We choose it to equal 

(2) 

our second trimer energy, i.e. Eoo = EL {a = oo). The lowest energy, 



E!^\a = oo), is avoided for regularization because the ground state energy 
could be influenced by the finite-range properties of the potential whereas 
the excited spectrum must have the features of the Efimov sequence. How- 



ever, in fig. 3.1 also the lowest energy coincides with the zero-range result 
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for \a\ = oo. 

The finite-range calculation only shows three bound trimer states of neg- 
ative energy in contrast to the infinite series of Efimov states in the free 
zero-range model. This reflects that the higher lying states are pushed up 
by the walls of the external field. They now appear in the region of positive 
energies, still corresponding to bound states, but now determined by the 
properties of the trap and not the two-body potential. 

For a < the trimer states become less bound and cross the trimer 
threshold. Deviations from the zero-range model is largest for the lowest 
state, n = 1, which is moved to the right in the plot, i.e. towards larger 
\a\. The states n = 2,3 move in the same direction, but the corrections are 
much smaller. 

For a > the straight diagonal line shows the atom-dimer zero-range 
threshold for binding. The corresponding finite-range result is lower and 
quantitatively in agreement with the dimer energy in eq. (2.11 ). Only effec- 
tive range corrections are needed for an accurate description of the dimer 
energy here. Above the atom-dimer threshold a quasi- continuum is present 
arising from the dense atom-dimer spectrum confined by the harmonic oscil- 
lator potential. The small spacing between these levels is then of order htu. 
The three Efimov trimers near the atom-dimer threshold lie below the zero- 
range result. This is mainly explained by corrections on the dimer energy, 
as we discuss later. 



3.4.2 The \a\ = oo Spectrum 

We already commented on the energies at \a\ = oo above, let us now consider 
the quantitative features. The energies e!^^ /fko for n = 1, . . . , 7 are given 



in the first row of tab. 3.1 The normal Efimov scaling sequence without 



a trap, is also shown in tab. 3.1, It can of course only describe the lowest 
bound Efimov states, n = 1,2,3, and not the trap-like state for n > 4. 
Our finite-range results agree well with these zero-range values. The ground 
state energy, n = 1, is remarkably close to the zero-range result. The purely 
attractive Gaussian interaction gives a slightly lower energy, but one could 
have expected more model-dependence for the ground state. The n = 3 
energy is slightly higher, which is due to the trap, see below. The energy 
ratios for the finite-range calculations are {e!^'^ — Ecm) / {e!P~^^^ — Ecm) = 
530,522 for n = 1,2, also in agreement with the zero-range result, e'^'^^^" = 
515.0. 

In [JHP02] the energy of three identical bosons in a harmonic trap was 
considered. Only zero-range interactions and s-waves were included. In the 
limit \a\/bt = oo the three-body energies, e!^\ are given semi-analytically 
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from the equation 



-arg < 



so 
2 



r(i-so) 



(3.3) 



where sq = 1.00624i and is a free regularization parameter or three-body 

(2) 

phase shift (equivalent to -Boo)- We fix = 0.4047r to match E}p of our 
finite-range calculations. The resulting spectrum is given in tab. 3.1 All 
energies agree perfectly with our finite-range results, the deviations being 
due to numerical errors only (except for the ground state). 







Efimov states 




Trap-like states 


n 


1 


2 


3 


4 5 6 7 


FR (trap) 


-3.744 


• 10^ -7067 


-12.04 


2.78 5.10 7.05 7.47 


ZR (free)t 


-3.64 • 


10^ -7067"^ 


-13.72 




ZR (trap)^ 


-3.64 • 


10*5 -7067"^ 


-12.20 


2.76 5.05 6.97 7.20 



t Zero-range model in free space, i.e. universal 515.0 scaling. 

Zero-range model in trap, |JHP02I . with 9 = 0.4047r equivalent to -Eoo- 
* This value was used to fix the three-body scale, Eoo ■ 

(n) 

Table 3.1: Efimov trimer energies E>p in a trap in units of hto for a = oo. The 
finite-range (FR) calculations are compared with the zero-range (ZR) theory. 



3.4.3 Trap-Like States 

We already commented on the trap-like states for \a\ = oo. Let us now 



consider the positive energy spectrum for a < in the upper part of fig. 3.1 
This is a "quasi-continuum" with the spacing of order fvuj. The asymptotic 
limit at small \a\ is simply the oscillator spectrum for three particles which 
remains when the two-body interaction is negligibly small. The lowest pos- 
itive energy level for a = is 3 x 2>hw/2 and the excitation spectrum is 
obtained by adding 2hL0. The almost vertical lines for the Efimov states, 
continuing into the positive energy region, cross the oscillator states. For 
small \a\ the coupling is also small and the avoided crossing appears like 
true crossings. As a increases the couplings increase and smooth avoided 
crossings appear. For a > a similar spectrum exists, but the huge amount 
of avoided crossings is not computed. 



resem- 



More generally, the avoided crossings in the upper part of fig. 3.1 
ble the Zel'dovich level rearrangement for two particles, where an attractive 
long-range potential is perturbed by a strong attractive short-range poten- 
tial |CKR+07| . In our case the long-range interaction is the harmonic trap. 
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and we have three particles instead of two. Also, in FJT^09 we showed 
how to calculate narrow resonances of three-body systems. The approach 
was to discretize the continuum by introducing a large artificial oscillator 
trap. When varying the trap size avoided crossings occurred, and resonance 
energies and widths could be extracted. Such avoided crossing are closely 



related to the crossings in fig. 3.1 The only difference is whether the trap 
or interaction is varied, while the other is fixed. 

Recently, Efimov physics in a finite square box with large positive scatter- 
ing length has been investigated within effective field theory. In the preprint 
|KH09a] it is concluded that "By decreasing the box size, the binding en- 
ergy decreases and eventually the state is shifted into the positive energy 
regime. The finite volume corrections are most important for the shallowest 
states which are largest in size and feel the finite volume first". These two 
conclusions agree with our results above where all but the three lowest Efi- 
mov states were pushed up by the harmonic trap. It is also consistent with 
|JHP02j . The behavior is easily understood from the fact that a confining 
volume induces extra kinetic (zero-point) energy. However, in the published 
version |KH09b| the figures were changed and the conclusion was replaced 
by "If the box size is decreased, the binding increases[ed.]" . This conclusion 
contradicts all physical intuition and is not elaborated in |KH09b . 



3.4.4 Finite-Range Borromean Window Corrections 

The differences between finite- and zero-range results for o < can be 



extracted from fig. 3.1 We focus on the measurable quantities expressing 



that an energy threshold has been crossed. 

(n) 

The finite-range results for the negative critical scattering lengths Up 



are shown in tab. 3.2 for n = 1, 2, 3 together with the corresponding effective 



ranges rP. The related zero-range values were calculated from eq. (|3.l[) 



with the regularization fixed by the finite-range energies e!^^ — Ecm from 



tab. 3.1 The values are shown in fig. 3.2, where we plot the relative shift in 



critical scattering lengths as functions of effective range. The systematics is 
discussed below. 

We note that the first threshold, is also significant by marking the 
interval of scattering lengths where two particles cannot bind while three 
particles can. This Borromean window between the threshold and \a\ = oo, 
has been calculated in [ MFK^OOj for a number of different radial shapes of 
the two-body potentials. Their results are expressed as ratios, S = 53/(721 of 
critical strengths required precisely to bind the two and three body systems, 
respectively. This ratio is limited to 2/3 < S" < 1, but for most ordinary 
potentials S ~ 0.8. We find S = 0.793 for the the Gaussian potential in 
agreement with 0.79 in |MFK"'"00 . With the critical potential strengths 



given m |MFK+00| we calculate Up and Rp ■ The results are compared to 
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Efimov states 




n 


1 


2 


3 


(n) 


-4.376 


-74.09 - 


-1631 


(n) 


-3.089 


-71.10 - 


-1625 


Kp 


1.670 


1.449 


1.436 



Table 3.2: Critical negative scattering lengths, a-p , for binding of the n'th Efimov 
trimer state. The finite-range interaction (F) is a Gaussian with range rp. The 
zero-range (Z) results are fixed by the trimers energies at \a\ = oo. The effective 
range at the position for binding is also shown. All values are in units of tq. 



3.1 



the Gaussian values and the zero-range results in fig. |3.2 As in fig. 
the effect of the positive effective range is a systematic shift of the critica 
scattering length towards larger absolute values. The difference decreases 
with increasing n which implies that most significant effects are related to 
small |a|-values. The available systematics in fig. 3.2 can be described by a 
straight line, i.e. 



a 



W _(n) 

Ui ry 



An) 



1.3 ^ 



(n) 



(3.4) 



The wide range of potential shapes used, indicate that this result most 
probable is model-independent. Although the conclusions are based on cal- 



culations with Re > only, we speculate that eq. (3.4) has the same form 



for Re < 0. This implies that the shift of critical scattering length changes 



sign, i.e. the energies move to the left in fig. 3.1 



The shift in eq. (3.4) should be measurable for cold atomic gases as 



positions of extrema of the three-body recombination rate near a Feshbach 
resonance |KMW+06l . For broad resonances the background channel (the 
van der Waals interaction) will determine the effective range. Thus Re will 
be positive and of the order l^dw^ ^-nd the shift will be toward larger \a\. On 
the other hand, for narrow Feshbach resonances Re is negative and given by 
eq. (2.14), so the shift will be towards smaller \a\. 



3.4.5 Finite-Range Atom-Dimer Corrections 

The effective range corrections are more complicated for a > because the 



atom-dimer threshold is also shifted, see fig. 3.1 This effect is investigated 
in fig. |3.3| (upper part) where the trimer energy is now shown as function of 
the finite-range dimer energy, Ed. Then the thresholds coincide for finite- 
and zero-range models. 

The ground state size is comparable to the range of the attractive finite- 
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0.1 



0.01 




0.001 



0.001 



Gaussian 
Yukawa 
Exponential 
Gaussian 
Morse 
Poschl-Teller 
Linear fit 



0.01 

RF/\aF\ 



0.1 



Figure 3.2: The relative shift between finite- and zero-range critical scattering 
lengths, and a^' as function of the corresponding critical effective range R^p\ 
The Gaussian data from fig. |3.1| are shown with numerical error bars. The other 
data are for the ground states of different finite-range potentials from [MFK+00 . 
The zero-range energy scale E'^^{a = oo) was chosen equal to the individual finite- 
range energies at \a\ = oo. The fit is given in eq. (3.4). 



range potential, and hence decreases more strongly with increasing Eo- As 
expected this state remains in the discrete region below the atom-dimer 
threshold. The energies of the excited trimer states n = 2,3 are in perfect 

(2) 

agreement with the zero-range prediction. Only Erp is slightly below the 
zero-range model near the atom-dimer threshold. We conclude that the 
major shift in trimer energies and threshold scattering lengths are due to 
effective range corrections to the dimer energy alone. 

This is consistent with the explanation in |Efi70| IEfi90| that these Efi- 
mov states disappear into the atom-dimer continuum. As also described 
in |BH06] this implies that their structure, as the threshold is approached, 
converges to the dimer-state with a loosely bound atom at large distance. 

This understanding is tested by computing the root-mea n-squ are hyper- 
radii prms = (p^)^''^! ec[- (2.27), shown in the lower part of fig. 3.3 as function 
of Eu. We first consider the n = 3 state. When E^ increases, both the en- 
ergy and the size of the Efimov state decreases. However, when the state 
approaches the atom-dimer threshold, the structure changes rather abruptly 
towards a dimer and a free atom, and hence the radius increases correspond- 



ingly fast towards the upper limit, bt, defined by the trap. From Fig 3.3 
this occurs at {cip ,Rp ) = (84.59, 1.423)ro. This is consistent with the 

zero-range result = 80.9ro (indicated by the arrow), since the effective 
range is small compared to the scattering length. The energy of the n = 2 
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-2 2 

sgu{a)\E2/nuj\^/^ 



Figure 3.3: Top: Three-body energies, e'^^\ as in fig. 3.1 but plotted as function of 
two-body energy instead of scattering length. On the right part (a > 0) £^2 = Eo 
is the weakly bound dimer, eq. (2.11 ), while on the left part (a < 0) i?2 = Ey is the 
shallow virtual dimer (also eq. (2.11 )). Bottom: The root-mean-square hyper-radii, 



Prms) as function of two-body energy. The zero-range predictions of the crossing of 
the Efimov states with the atom-dimer threshold are indicated by arrows. 
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Efimov state in fig. 3.3 is also in good agreement with the zero-range result. 
However, when the state approaches the atom-dimer threshold no drastic 
increase in size is observed, as for ordinary Efimov states. This can be as- 
cribed to the effective range which is comparable to the scattering length 

(2) -(2) 

in this regime, namely {(ip , Rp ) = (4.1, 1.18)ro. This effect needs a more 
careful analysis in the future. 

3.4.6 Effective Range Corrections in EFT 

After the publication of the above results in |TFJ08c| , many of these results 
have been confirmed qualitatively by effective field theory |PJP09| IPla09 



In this subsection we sum up the common effects and also give quantitative 
comparisons. In |PJP09j the linear effective range corrections to the three- 
body energies were calculated for Rg > 0. The main features are presented 



in their fig. 1, which essentially is identical to fig. 3.3 in this chapter. The 
following conclusions are given in |PJP091 IPla09j : 

i) The three-body spectrum at \a\ = oo is unperturbed when introducing 
a linear effective range correction. 

ii) For a < the critical scattering lengths at the trimer threshold move 
towards higher |o|. The effect is largest for the lowest states. 

iii) For a > the effects are in general small (after the dimer energy 
correction is taken into account). The trimer only becomes slightly 
more bound. 

All these qualitative effects are in perfect agreement with the previous con- 
clusions of this chapter. 

Let us now turn to a quantitative comparison of the shift in critical 
scattering length for a < 0. Since |PJP09| use the two-body bound state 
pole (i.e. energy) instead of scattering length, we must repeat the analysis 



of subsection 3.4.4 using k ^ instead of a. We include the full effective range 



corrections in k as given as in eq. (2.10). The result is shown in fig. 3.4 



The use of compared to fig. |3.2[ mostly affects the n = 1 data. The 



linear relationship becomes marginally better using k ^ instead of a. The 
linear fit gives 



1/kz 

We extract values for the same shift in fig. 1 of |PJP09] . By reading off 
the critical values for the lowest curv^Ewe estimate their lowest order result 
to be klo/'^o = —5.68 and the next-to-lowest order result (including linear 



^Note that they use the notation 7 instead of k, and use a scaled plot. 
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5 

/ n = \ ■ 


0.1 




n = 2 / - 








J 0.01 


n = 3 










0.001 




Linear fit 1 

■■ 




0.001 


0.01 0.1 

Rf\kf\ 


Same data as fig. 


3.2 but 


with scattering length < 



binding wave number, k ^. Legends are the same as in fig. 3^ and the fit is given 
in eq. (13. 5|. 



Re) to be knlo/i'^o = —5.19. Here kq is an arbitrary momentum scale. The 
effective range was cliosen to be constant, ReHQ = 0.01. Tliis gives 

^jVLO -^LO y^^l^^^^l ^ 1.8^0.5. (3.6) 



Tliis value agrees witli tlie coefficient in eq. (3.5). Tlie difference can be 
attributed to numerical and read-off errors. Also, tlie models are based on 
very different backgrounds and formalisms: In jPJP09j only linear effects 
in Rf, are included in an EFT approach, while our finite-range calculations 
in principle can have minor effects from higher-order shape parameters (i.e. 
potential model-dependence) and the trap. Taking all these remarks in 
consideration, the agreement is actually extraordinarily good. 

The discussion above is only based on a single value from [PJPOD], i.e. 
the predictions of the linear relationship in eq. ( |3.5[ ) is not confirmed by EFT 
yet. Such work is currently in progress, and preliminary results confirm the 
linear prediction with the correct coefficienfj^ 



3.5 Conclusions and Outlook 

We have used a basis of correlated Gaussians to calculate the spectrum for 
three identical bosons for a finite-range potential and an external trap. The 
scattering length was tuned to arbitrary values, and the many excited states 
were computed accurately as function of scattering length. The universal 

■^Private communications, David L. Canham, Universitat Bonn, Germany. 



42 



CHAPTER 3. EFIMOV PHYSICS: FINITE RANGE AND TRAP 
EFFECTS 



spectrum at infinite scattering length was in perfect agreement with analyt- 
ical models. 

We then focused on two thresholds at negative and positive scattering 
lengths, i.e. the Borromean window where two particles cannot bind while 
three can form many bound states, and the atom-dimer threshold where 
dimers can bind and an atom may be bound to the dimer. 

We have extracted the universal behavior of both these thresholds, in- 
cluding corrections expressed in terms of effective range divided by scattering 
length. We also conclude that effective range corrections to Efimov physics 
and Borromean binding are two aspects of the same effect and these two 
regions can be connected quantitatively. The linear shift in critical scatter- 
ing lengths at the trimer threshold agrees quantitatively with effective field 
theory. The main effect is that the universal scaling factor becomes smaller 
for the lowest Efimov states (smallest scattering length) on the negative 
scattering length side. 

We also show that the structure as the atom-dimer threshold is ap- 
proached is a weakly bound atom moving away from the dimer. The main 
effective range correction to the trimer energies near the atom-dimer thresh- 
old then comes from the dimer only. This seems to induce a larger scaling 
factor for the lowest Efimov states. 

The results have observable consequences for three-body recombination 
rates in atomic Bose-Einstein condensates and should be taken into account 
when the scattering length is comparable to the effective range 
lKFM+091 . The effects are also applicable for halo nuclei and Borromean 



systems. 

Future investigations could repeat a similar analysis with potentials of 
negative effective range, which could be useful for modeling narrow Feshbach 
resonances. We speculate that the linear relationship eq. (3.4) also holds for 
negative effective range. This conjecture is strongly supported in the next 
chapter, where we analytically derive bounds for the Efimov effect in the 
case of large negative effective range. 
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Chapter 4 

Conditions for Efimov 
Physics 



4.1 Introduction 



We have seen in previous chapters how universal scaUng properties in three- 
body systems arise when the scattering length a is much larger than the 
range ro of the underlying two-body potential. In this regime several three- 
body observables are universal. The universal scaling of Efimov trimers is 
usually said to exist for rms-sizes between vq and a |Efi731 lEfiQH IBH061 
|KMW+06| . The effective range Re from a low-energy phase shift expansion 
is sometimes used instead of ro in this statement |FJ01bl ITFJ08c| IPJP09] . 
This ambiguity occurs because ro and Rg are often of the same order. How- 
ever, for narrow Feshbach resonances in atomic gases Rg can be much larger 
than ro |B JK05] . and the implications for such systems need to be explored. 

Zero-range models, in particular in combination with the hyper-spherical 
approximation |FJ01b| IJon04| IBH06j . have been successful in semi-quanti- 
tative descriptions of three-body systems in the universal regime. Semi- 
rigorous finite-range corrections have been attempted by including higher- 
order terms in the effective range expansion |FJ01b . P JP09] as a step towards 
the full finite-range calculations as in |SEGB02l [TFJO ScJ while maintaining 
the conceptual and technical simplicity of the zero-range approximation. 

The obvious generalization of the zero-range model, eqs. (2.48 )-(2.50 ), 
is to substitute —1/a with —l/a + {Re/2)k'^, where k is the two-body wave 
number, in the relevant expressions for the logarithmic derivative of the total 
wave- function at small separation of the particles. However, in three-body 
systems neither the two-body wave number nor the small separation are 
uniquely defined, and rigorous inclusion of all terms of the given order is non- 
trivial. The lack of rigor in previous works could have serious implications for 
applications where finite-range effects are important, such as the stability 
conditions for condensates in traps, and in particular for Efimov physics. 
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Experimental progress |KMW"'"06 |KFM"'"09] will soon require this increased 
accuracy near the boundaries of the universal regime. 

In this chapter we consider a system of three identical bosons near a Fesh- 
bach resonance in the universal regime with large scattering length usually 
described by model-independent zero-range potentials. We employ the adi- 
abatic hyper-spherical approximation described in chapter |2] and derive the 
rigorous large-distance equation for the adiabatic potential for finite-range 
interactions. The equation is suitable for the analytic studies of the finite- 
range corrections in the three-boson problem. We investigate the finite-range 
corrections to the adiabatic potential and the non-adiabatic term and com- 
pare with the zero-range approximation. The effective range correction to 
the zero-range approximation must be supplemented by a new term of the 
same order. The non-adiabatic term can be decisive. Efimov physics is al- 
ways confined to the range between effective range and scattering length. 
Our analytical results agree with numerical calculations for realistic poten- 
tials. 



4.2 Adiabatic Eigenvalue Equation 



We first generalize the adiabatic eigenvalue equation, eq. (2.50), to arbitrary 



finite-range potentials. The derivation is similar to that of chapter |2] but 
with some modifications. 

We consider three identical bosons of mass m and coordinates rj inter- 
acting via a finite-range two-body potential V, where we assume V{rjk) = 



for r 



'jk 



fk\ > fo- Only relative s-waves are included. We use 



the hyper-radius p 



12 



13 



)/3, eq. 

= k,; 



(|2.27|) and hyper- angles 
In the fol- 



tanaj = {rjk/ri^^jk))V^/2, where r^^Q-fc) = \ri - {rj + rk)/2\. 
lowing we shall use one set of coordinates and omit the index. 

The adiabatic hyper-spherical approximation treats the hyper-radius p 
as a slow adiabatic variable and the hyper-angle a as the fast variable. The 
eigenvalue X{p) = i^'^{p) — 4 of the fast hyper-angular motion for a fixed 
p serves as the adiabatic potential for the slow hyper-radial motion. The 



eigenvalue is found by solving the Faddeev equations, eq. (2.41 ), for fixed p, 



^ = -2un[ip]. 



Here ip{p,a) is the Faddeev hyper-angular component, 

U(p,a) = V(V2psina) — — 5— 
is the rescaled potential, and 



(4.1) 



(4.2) 



nmp,a) 



V3 J\^-a\ 



ip{p, a) da' 



(4.3) 
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is the operator that rotates a Faddeev component into another Jacobi system 
and projects it onto s-waves, eq. (2.44). The hyper-angular wave-function 
of the three-body system is 



$(/9,a) 



V^(p,a)+27^[^](p,a) 



(4.4) 



sin 2a 

In the adiabatic approximation the hyper-radial function /(p) satisfies the 
ordinary hyper-radial equation, eq. (2.36), with the effective potential 

2m 



1/4 



Q 



(4.5) 



where Q is the non-adiabatic term and ^ is normalized to unity for fixed p. 

We first divide the a- interval [0; 7r/2] into two regions: (I) where U ^ 0, 
and (II) where U = 0. The regions are separated at a = ao where sinao = 
ro/{V2p). 

Region (II): Here U = and we have the free solution to eq. (4.1), 



iP'^a) =AA(p)sin(z^(a- -)), 



(4.6) 



with the boundary condition iI)^^{tt/2) = and normalization N{p). 
Region ( I) : We now restrict the problem to be outside the "collapsed" region 
defined hy p > pc = V^rQ. This implies that ao < 7r/6. Since a < ao the 



rotation operator, eq. (4.3), only integrates over the free solution tp \ see 
fig. 4.1, Equation (4.1) then simplifies to 



da^ 



ip^ = -2Un[i;^^]. 



(4.7) 



The solution is = ip^'^ — 2TZ['il>^^] , where ip^^ and —21Z['il)^^] are homoge- 
neous and inhomogeneous solutions, respectively, ip^^ is the regular solution 
to 



m dr"^ 



m 



V''" = 0, 



Vp{r) = V{V2p 



sm 



/2p 



)), 



where kn 



z//(\/2p) and r = \f2ap. When a — > oq, 
'pi^^ oc sin(A;pr -|- (5p), 



(4.8) 
(4.9) 

(4.10) 



where the modified phase shift 6p{kp) arises from the modified two-body 
potential, Vp. 

The solutions <D in region (I) and (II) are now matched smoothly, giving 



da 



Ih 



d_ 
da 



In (V'^^ + 27^[V'■^^]) 



(4.11) 



"0 
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Figure 4-.1: Integration limits for the rotation operator eq. (4.3 1. For p > pc the 



shaded regions does not overlap and the three-body problem reduces to an effective 
two-body problem. When a < ao < 7r/6 the integration limit for a' is strictly 
outside the potential range [0 : ao] giving 7?.['0] = 7?.[V'^^]. 



After inserting eqs. (|4.6| and (4.10), this equation becomes 

= kpCoi5p{kp), 



1 



(4.12) 



V2p sin(z.f) 
which defines u as function of p. This generalizes the — l/o expression for 



the lowest-order zero-range model, eq. (2.50). 



The right-hand-side deviates from the zero-range approximations |Jon04| 
IPJP09) which use the normal two-body phase shift, i.e. the right-hand-side 
is kpCot 5{kp). We use the rigorously defined phase shift 6p for Vp instead of 
6. Our result is exact for p > pc- 



4.3 Effective Range Expansion 

In the limit p » pc, the /o-dependent potential, Vp{r), approaches V{r), and 
consequently 5p approaches S. The p-dependent low-energy effective range 
expansion corresponding to Vp is then to second order 

kp cot 6p{kp) « — ^ + ^^^kl (4.13) 



a{p) 



where a{p) and Re{p) are functions of 1/ p^ that converge to a and Re for 
p oo. Up to l//3^ in eq. (4.13) we get 

1 



a{p) 



1 1 
a zp"^ 



Re{p) ~ Re 



(4.14) 
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The model-dependent expansion parameter Ry, which we shall call the 
"scattering length correction", is found to 



V' {r)r^u{r)'^ dr, 



(4.15) 



where u is the zero-energy two-body radial wave-function, asymptotically 
equal to 1 — r/a. Equation ( 4.12[ ) then becomes 



1 -z^cos(z^f) ^sin(zy| 



V2p 



smu^T 



1 Re 



Rv 
2p2- 



(4.16) 



This equation without the last two finite-range terms has the well-known 
purely imaginary solution uq = 1.00624z, or Aq = —5.0125, for ^/2p <^ 
\a\. This solution gives Ves oc — l/p^ which is the basis of Efimov physics. 
The i?e-term was included in |FJ01b| IPJPOQj , but not the model-dependent 
i?y-term. The latter term makes the finite-range corrections to the zero- 
range adiabatic eigenvalues explicitly non-universal. The last two terms in 
eq. (4.161 restrict the solution Aq to the region |i?o| <^ ^/2p <^ \a\, where 



^0 = 



Rv, 



(4.17) 



as seen in fig. 4.2 where the lowest solution to eq. (4.16) is shown for dif 



ferent parameter choices. Thus, naively one would think that the lower 
limit for Efimov physics is determined by the model-dependent length | .Ro I • 
However, we will show later that Q restores universality and recovers the 
model-independent effective range, i?e- 

First, to illustrate the necessity of both l/p^-terms in eq. (4.16) we con- 
sider a large negative effective range corresponding to a narrow Feshbach 
resonance |BJK05j . To model the large \Re\ we pick the attractive potential 
with barrier, eq. (2.23), 



V{r) = D sech^ x— + ^exp -2(x 



(4.18) 



where D = -138.27, B = 128.49 in units of h^/{mrl), and x = 4.6667. 
The potential is negligible outside the range rg. The low-energy parameters 
are a = 556.88, Re = -142.86, Ry = 73.031, and Rq = -0.71 in units of 
To- In fig. 4.3 we compare \{p) obtained by exact numerical solution of the 
Schrodinger equation [SEGB02] containing the interaction eq. (4.18) with 



^The integral representation of the effective range, eq. (2.211, is usually proved by 



integrating and subtracting the radial equation for two different energies, comparing to 



the free solutions, and then letting the energies go to zero. Equation 14.151 can be proved 



in a similar manner by using different potentials {V and Vp) at the same energy and letting 
p <x. 
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Figure 4-2: Adiabatic eigenvalues A(p) from eq. (4.161 as function of hyper-radius 
p, for large scattering length a and negative Re- Different values of the model- 
dependent length Rq are used, showing that the universal solution Aq exists in the 
region |-Ro| ^ V^P <S^\a\- Lengths are in units of \Re\- 



the solution of eq. (4.16). In the zero-range model (including only 1/a), 
the — /O^ divergence for large p is below the numerical solution. At small 
distances, A approaches Aq, above the numerical solution. Inclusion of the 
i?e-term, as in |FJ01bl IPJP09] . provides a better large-distance behavior 
(since the dimer binding energy is corrected), but overshoots dramatically 
for y/2p < a by approaching A = — 4. Including consistently both Rg- and 
-Ry-terms leads to complete numerical agreement with the exact numerical 
solution except for very small p-values near pc where higher-order terms are 



needed in eq. (4.16) 



4.4 Non- Adiabatic Corrections 



We now show that the non-adiabatic term restores model-independence and 
recovers \Re\ as the limit for the region of Efimov physics. For simplicity we 
only consider the limit |a| = oo and assume |i?o| ^ I - ^e I We first consider 
-v/2p *C l-Rel (region A in fig. 4.2). Expansion of eq. (4.16) to first order in 
— uq) gives a small constant correction 



^0 



Rq 

1^0 Re 



l + 0(f ) 

Up 



(4.19) 



This correction is marked by the arrow in fig. 4.3 (it is out of the range of 



fig. 4.2). This gives 

VeSip) 



2m 



1/4 - 2Ro/Re 

p2 



Q 



(4.20) 



^This assumption holds for the potential eq. (2.231 as well as many other potentials 
with large effective range. 
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Figure 4-3: Exact numerical adiabatic eigenvalues X(p) for a potential with bar- 
rier, eq. (4.181 (solid red line), compared to solutions of the eigenvalue equation, 
eq. (4.16). The zero-range model (dotted line) includes only a. Crosses include 
a, i?e-terms and circles include a. Re, i?y-terms. The inset shows details around Aq. 
The arrow indicates the effect of the correction Aq — 2i?o/^e- The shaded region is 

P < Pc- 



To evaluate Q we note that a large negative effective range (for \a\ = oo) 
implies that the two-body wave-function u is localized mainly inside the 
potential range, as discussed in chapter |2j Then the angular three-body 
wave-function $ can be approximated by u/ sin(2Q!), i.e. u and ne- 

glecting rotated terms The result is Q = c/p^, where c ~ —5/4 as 



confirmed numerically. This term cancels the main 1/p^-part in eq. (4.20) 
and hence prohibits Efimov physics for ^/2p <C \Re\- The intuitive reason 
is that the two-body wave-function is essentially zero outside the potential, 
despite the large scattering length, and hence three particles can not interact 
at large distances. 



When V2p> \Re\ (re gion B in fig. 4.2^ we find 



.„ + .„,„^(l + 0(5^)). (4.21) 
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where 

Co = -A -, ^ ^ ~ -0.671. 4.22 

cos(z^06) - cos(i/o 2) + 2 sin(fo 2) 

This gives the effective hyper-radial potential 

The l//9^ dependence of the correction to the Efimov potential was 
expected |Efi91j . The model-independent term proportional to Re/ was 
recently calculated in |PJP09| . However, we also get a model-dependent 
term Ry / which is of the same order. Q generally receives contributions 
both from distances inside and outside the finite-range potential. Zero-range 
models only have the external part of the wave-function, which depends on 
p only though the eigenvalue v{p). The zero-range result for Q is then 

QzR = Mo(|)^ = Mocg.o^^, (4.24) 
where Mq = {<^\d'^<^ / dv'^)\u=uo- This fourth order correction can be ne- 



glected in eq. (4.23), as was done in |PJP09j . However, the internal part 
of the wave- function contributes to order 1/p^. To estimate this 1/p^-term 
we take the analytically solvable finite square well potential of range tq and 
\a\ = oo. This fixes Re = and Ry = n^7r^ro/24 where n is the number of 
bound states (including the zero-energy state). We find 

Q^o. = coyl{^-2Ry)-^^, (4.25) 



neglecting 1/p -terms. The model-dependent i^y-terms in eqs. (4.23) and 



(4.25) cancel exactly, giving 



So the effective potential receives a Re/ p^ correction where the model- 
dependent coefficient is different from zero-range models |PJP09| because 
of the inclusion of Q. We also expect the i?i/-terms to cancel for general 
potentials. In conclusion, the Efimov effect persists for \/2p ^ |-Re|- 



4.5 Atom-Dimer Potential 

We have seen that model-dependent corrections to Aq are cancelled by equiv- 
alent terms in Q. A similar effect occurs for the atom-dimer channel poten- 
tial. Suppose the binding energy is B£) = h^kl^^/m with corresponding wave 
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number kf) > 0. Then v = ik^V^p is an asymptotic solution to eq. (4.12) 



and A diverges as —p corresponding to a bound dimer and a free particle. 



For this solution, the effective range expansion eq. (4.13 ) does not hold, since 
asymptotically kp — > ik^ is finite. Instead eq. (4.8) reduces to the radial 



two-body equation, with a normalized bound state s-wave function ti/)(r). 
Treating Vp — V oc 1/p^ as a perturbation gives the correction 

Since TZ[ip] is exponentially small for the atom-dimer solution, Q can be 
computed using the unperturbed wave-function ip = ^U£){\/2ap), giving 

1 1 1 

Q = -—+ UD{ru'D+r\'[))dr^ + 0{^). (4.28) 

The term — l/(4p^) for the atom-dimer solution is well-known and sometimes 
referred to as the Langer correction |NM99] . 

By using the two-body radial equation and partial integration the two 
integrals in eqs. (4.27 ) and (4.28 ) cancel. Thus the 1 / /?^-terms in the effective 
potential eq. (4.5) cancel exactly, giving Fofr(/o) = —Bd up to order l/p^. 



Thus Veff only depends on through B^. 

4.6 Applications near Feshbach Resonances 

The effective range near a narrow Feshbach resonance can be estimated as 



in eq. (2.14). As an example with a noticeable effect we take the alkali 
atoms K with the very narrow Feshbach resonance at i? = 825G having 
parameters /S.B = — 32mG, A^u = — 3.92/^^, and a^g = — 36ao |DZF"'"07 



This gives the large effective range Re = —2.93 x lO'^aQ. For ^^K, rg is of 
the order of the van der Waals length /^dw = 1-29 x IO^oq [BHOG'. Since 
I -Re I ^ ^'O) I -Re I determines the lower limit for Efimov physics and corrections 
to the universal regime are of order Re/ a (not ^vdw/^)- Thus, the window 
for universal physics is reduced. 

Another example is the Feshbach resonance at Bq = — 11.7G (Ai? = 
28.7G, a^g = 2.30ao, Ap = 22pB [CGJT09]) in ^^^Cs where jKMW+06| 
apparently have observed a single Efimov state. In this case the effective 



range estimated from the Feshbach model, eq. (2.14), is Re = — 0.34ao. This 
is much smaller than l^dw = 202. Ooq |CGJT09| . thus l^dw is the dominating 
scale for effective range corrections. The linear effective range corrections of 
chapter [3] could be applied here. 

4.7 Conclusions and Outlook 

We have considered a three-body system of identical bosons with large scat- 
tering length modeling a Feshbach resonance. The Efimov physics occurring 
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in this universal regime is customarily accounted for by zero-range models. 
We used the adiabatic hyper-spherical approximation and derived rigorously 
a transcendental equation to determine the adiabatic potential for a general 
finite-range potential. We have solved this equation for large scattering 
length, investigated finite-range effects, and compared with exact numerical 
results. 

Inclusion of the effective range correction to the adiabatic potential is 
insufficient in general. Crucial corrections of the same order must also 
be included from both a "scattering length correction" and the non- 

adiabatic term. These two contributions may separately be large but they 
tend to cancel each other. Accurate results in zero-range models must ac- 
count for these new corrections. In conclusion, for large negative effective 
range the window for Efimov physics is precisely open between the effective 
range (not the potential range) and the scattering length. 



54 



Chapter 5 

A^-Body Efimov Effect 



5.1 Introduction 

The Efimov effect appears in quantum tliree-body systems wlien attractive 
interactions between at least two pairs of particles are such that the scatter- 
ing length is much larger than the range of the interaction; in other words 
two of the three two-body subsystems are close to the threshold of binding. 
Under these conditions a characteristic series of weakly bound and spatially 
extended Efimov states appears in the system. These states appear due to 
specific long-range two-body correlations between particles caused by the 
large scattering length. The effect is easiest to see in the hyper-spherical 
adiabatic approximation, chapter |2] where the slow adiabatic variable is the 
hyper-radius p. It has been shown |FJ93j that close to the two-body thresh- 
old the effective adiabatic potential V^ff(p) is attractive and asymptotically 
proportional to the inverse square of the hyper-radius, 

.»(.)^-|,^, 

where m is the atomic mass, and ^ is a constant depending on masses of 
the particles |NFJGOl| INFJ98] . For three identical bosons this constant is 
given by ^3 = 1.00624 leading to the scaling factors 

53 = exp(7r/^3) ~ 22.7 and 5| ~ 515.0, (5.2) 

for radii and energies, respectively. This is equivalent to a geometric series 
of bound states with exceedingly small energies. En oc e"^'^"/^^, and exceed- 
ingly large root-mean-square hyper-radii, {p^)]!"^ oc e™/''^, where n is the 
state number. 

It was shown by Amado and Greenwood that in an A^-body system with 
N > 3 the Efimov effect does not exist at the — 1 threshold |AG73j . In 
particular, for A = 4 there is no infinite sequence of four-body states at 
the trimer threshold. This does not rule out an A-body Efimov effect at 
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the two-body threshold where the scattering length a is infinite. Of course, 
weakly bound A^-body states at a = oo can not be true bound states, since 
at this point the clusters with three and higher number of particles are 
generally deeply bound. The states would then be unstable due to lower 
lying thresholds and would decay into deeply bound cluster states. 

However it has been sugg ested in |SFJn2] that a sequence of meta-stable 
A-body states with the characteristic exponential energy dependence can 
yet show up at the two-body threshold. Using the A-body hyper-spherical 
method it has been shown that an A^-body system at the two-body threshold 
has a hyper-spherical adiabatic potential with inverse-square dependence. 
This peculiar adiabatic potential appears due to the same mechanism as for 
three particles and thus gives rise to A^-body Efimov states with a structure 
similar to that of three-body Efimov states: An (otherwise) uncorrelated 
system with very specific two-body correlations caused by the large scatter- 
ing lengths. 

This specific hyper-spherical adiabatic potential is not the lowest one as 
different bound clusters with lower thresholds create lower lying adiabatic 
potentials. However, although not truly bound, these A-body Efimov states 
might still exist as meta-stable states slowly decaying into clusters, much 
like the Bose-Einstein condensate states. The structure of the Efimov states 
is determined by the long-range two-body correlations and is thus quite 
dissimilar to the structure of clusterized states with short-range many-body 
correlations. Therefore the overlap between Efimov states and clusterized 
states should be small and consequently the life-time should be large. 

The conclusions about meta-stable A-body Efimov states were obtained 
in |SFJ02] in an hyper-spherical adiabatic approximation where the cou- 
plings to all other channels were neglected. In this chapter we report 
on a more realistic calculation of A-body Efimov states with a different 
method, namely the two-body correlated Gaussian stochastic variational 
method where no adiabatic approximation is assumed. We calculate spec- 
tra of trapped systems of A = 3, 4, 5, 6, and 7 bosons: For each system the 
calculations reveal a series of Efimov states where the energies and the radii 
exhibit the characteristic exponential dependence upon the state number. 
We obtain the new scaling factors for A > 3 and compare with analyti- 
cal results from the hyper-spherical adiabatic approach. Finally, we discuss 
observable consequences. 



5.2 Numerical Results 
5.2.1 The System and Procedure 

We consider a system of A identical bosons with mass m and coordinates 
rj in a spherical harmonic trap with frequency uj and trap length bt = 
y^h/{'muj). The two-body potential is an attractive Gaussian, eq. (2.20), 



56 



CHAPTER 5. iV-BODY EFIMOV EFFECT 



with range rg and depth Vq. As in chapter [3] we choose the trap to be much 
larger than the potential range, bt/ro = 3965. The depth is tuned around 
Vq ~ — 2.68^^/(mrQ) where the first bound state is at the threshold, giving 
\a\ > bt. 

The Schrodinger equation for the Hamiltonian in eq. ( |2.70 ) is solved 
numerically with the stochastic variational method of chapter [2 using the 
two-body correlated basis, eq. (2.77). The non-linear parameters in the 
wave-function are optimized stochastically by random sampling in a region 
that covers the spatial distances from ro to bt, which allows us to describe 
states over many length scales and energies. 

We look for A^-body Efimov states with the structure analogous to that 
of the three-body Efimov states. The spatial extension of such state should 
be much larger than the range of the potential ro and smaller than the trap 
length bt. 



5.2.2 Universal Energies and Radii 

The energies of Efimov states in a trap should be on one hand much larger 
than the typical energy scale of a cluster state, h?/{mrQ), and on the other 
hand smaller than the oscillator energy Hui = h? /(mb^). We have calculated 
the spectrum for N = 3,4,5,6, and 7 bosons. The calculated energies are 



shown in fig. 5.1 Indeed, in the indicated energy region each of the A-body 
systems have a series of states with exponential dependence upon the state 
number. En oc e~^'^"/'>^. The exponential fits give the scaling constants 
Sn = expvr/^jv, namely ^3 = 23.6(4), S4 = 5.47(40), ^5 = 2.45(4), Sq = 
1.93(2), and = 1.66(1). The indicated statistical errors are from the fits, 
while the uncertainty related to numerical convergence is not included. This 
is expected to be largest for A = 3, 4 which use only a few states for the fit. 
The A = 3 result agrees with the known analytical result of 22.694 |BH06) . 

The spatial extension of Efimov states in trapped systems must also be 
much larger than the interaction range ro and much smaller than the trap 
length 6f The average size of the state can be described with the (A^-body) 
hyper-radius p, 

N N N 

i<j i=l i=l 

where rij is the pair distance and R the center-of-mass. On fig. 5.2 are 
shown the calculated r.m.s. hyper-radii {p'^)}/'^ as function of state number 



n. The radii of the Efimov states identified on fig. 5.1 are reproduced well 
with the exponentials {p^)]!"^ oc e™/^^ where the parameters ^tv are taken 
from the fits on fig. |5.1[ Apparently all these states fall within the correct 
boundaries and the values of the radii follow the correct exponential trend. 
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Figure. 5.1: Energy En as function of the state number n for negative energy states 
of a system of N bosons in a harmonic trap with length bt interacting via an 
attractive potential of range tq with scattering length much larger than ht. The 
horizontal lines demarcate the region where the Efimov states can exist in trapped 
gases. The lines show the exponential fits of the form £"„ oc e~2'™/{]v drawn through 
the points in the indicated region. 



There are a few states in the = 6 and 7 systems with radii much 
smaller than those of the typical states in the series with similar energies, 
see fig. 5.2 Clearly these states are not Efimov states but rather relatively 
compact states with a different structure. We believe that these states are 
caused by missing convergence or numerical inaccuracies. 



5.3 Comparison with Analytic Results 

The geometric scaling found in the previous subsection indicates that an 
underlying 1/ effective potential exists. In this section we briefly outline 
the A-body results of |SSJF05al ISSJF05b] and find analytical scaling factors 
Sn = exp(7rn/^Ar) for arbitrary N . 



5.3.1 iV-Body Hyper-Radial Potential 

In |SSJF05a| ISSJF05b] the A-body zero-range problem was investigated 
within the hyper-spherical adiabatic approximation. The method is equiv- 
alent to the A = 3 discussion in chapter [2] still using the hyper-radius 
as the adiabatic variable and denoting all other coordinates by J7. The 
basic approach is now to use a Faddeev-Yakubovski expansion <I'(p, Jl) = 
Si<j "^(z^' '''ij) angular wave-function where each of the identical two- 

body amplitudes only depend on the distance between two-particles. This 
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Figure 5.2: Root-mean-square hyper-radius {p^)l/'^ as function of the state number 



n of a system of N bosons from fig. 5.1 The hues show the fitting curves of the 
form const xe'^"/^" drawn through the same points and with the same parameters 



as on fig. 5.1 



approximation is equivalent to our two-body correlated basis, eq. (2.77), 
used in the previous section. 

The adiabatic eigenvalue is denoted N) which is similar to the eigen- 
value u(p) for the three-body case. If the effective hyper-radial poten- 
tial is written in the form eq. ( |5.1[ ), the relation between ^(p) and ^{p) 
is 4z/ = 5 — 3A^ -|- 2 Y^— ^2 (note that .^^ is real). By imposing the zero-range 



boundary condition similar to eq. (2.48) and collecting similar terms from 



the two-body amplitudes (j), the adiabatic eigenvalue equation becomes 

1 B{D) + R{D) 1 



(5.4) 



This form is similar to eq. ( |2.50 ) in the three-body case. The terms A and 
B are found to be |SSJF05al ISSJF05b] 



A 



3Af-6^ 



vr 



3N-5 
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2cos{TTv)T[i> + 1) 



vr 



while the "rotated" terms R = 2{N - 2)R^i + UN 



2){N -3)R34. are 
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with the integrals /13 and I34 given by 



/■v^ ^ 1 

/i3 = y ^ dx vi + x (- 

/34 = y dx \/r+~X (1 - 



(3Af-ll)/2p_(3Af/2-4,l/2) 



(3Af-ll)/2p(3Af/2--4,l/2) 



(5. 



(x). 



Here P^"'*^ are the Jacobi functions, see e.g. [NFJGOll App. A]. Equa- 



tion (5.4) defines ^(/o) and hence the effective hyper-radial potential. 



5.3.2 Scaling Factors 



The solutions of eq. (5.4) was discussed in |SSJF05al ISSJF05b] . but only 
in the asymptotic limit — > 00. We now solve the equation numerically 
for finite and infinite scattering length. The solutions are defined by the 
roots of the real function 



2, _ _B{e)+R{ei 
Me) 



{5.9) 



which is calculated numerically and shown in fig. 5.3 for N = 3, . . . , 10. 



Each root defines a constant solution ^^r, leading to the corresponding scale 



factor Sn = exp(7r/^Ar) (see section 2.2.5). 
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Figure 5.3: Roots £,n of the function Z]s[{^) for boson number = 3, . . . , 10. Only 
the lowest adiabatic potential (largest ^) is shown. 



The values are summarized in tab. 5.1 and in fig. 5.4 where we also 
compare with the numerical results. There is generally very good agreement 
between analytic and numerical results. We note that the scaling factors 
decrease with increasing N. Already for four particles it has been reduced 
to 4.29 as compared to 22.7. We also show the asymptotic expression for 
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large A'", 




(3A-4)(3A-6) 



1 
4 



(5.10) 



which was established in |SSJF05al ISFJ03j . Our results are consistent with 
the asymptotic estimate. 



N 




Sn 




Sn{SVM) 


3 


1.00624 


22.6942 


515.028 


23.6(4) 


4 


2.15584 


4.29413 


18.4396 


5.47(40) 


5 


3.40602 


2.51523 


6.32638 


2.45(4) 


6 


4.74528 


1.93875 


3.75877 


1.93(2) 


7 


6.15939 


1.66537 


2.77347 


1.66(1) 


8 


7.63821 


1.50878 


2.27641 




9 


9.17426 


1.40837 


1.98352 




10 


10.7618 


1.33900 


1.79291 





Table 5.1: Scaling factors Sn — exp(7r/^Ar) for the A^-body Efimov effect. Both 
analytical and numerical (SVM) results are shown. Indicated errors are from the 
fits in fig. |5.1| The deviation from the analytical results is largest for = 3,4 
which use only a few states for the fit. 



5.4 A^-Body Recombination 

The A-body Efimov states could in principle be identified by their contribu- 
tion to the recombination rates in a cold gas as function of scattering length. 
As the scattering length is increased (using the Feshbach technique), the N- 
body Efimov states should produce peaks in the recombination rate when 
crossing the threshold similar to the three-body case. Such features would 
occur at critical scattering lengths spaced by the scaling factors Sn- 

The non-resonant background rate for A-body recombination can be 
estimated as follows. Let us first consider a three-body reaction. The rate 
of the loss of particles from a cold gas due to the three-body recombination 
reaction into a shallow dimer with the energy jima^) is given by Fermi's 
golden rule, 

dAo _ A3 27r 2 dv^ 

where the factor 3 is there since each recombination reaction removes three 
particles from the cold gas, Aq/6 is the number of triples in the gas of Aq 
particles, dvj is the number of final states (dimer plus atom) with relative 
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Figure 5.4: The scaling factors Sn as function of the boson number N. Both 
analytical roots and numerical (SVM) results are shown. The = 3 value is the 



well-known factor 22.7. The asymptotic form is given by eq. (5.101 



momentum qj = 2/(\/3a) and the relative kinetic energy 



E 



f 



2Um) 



ma" 



(5.12) 



and Tfi is the transition matrix element from the initial three-body state to 
the final atom-dimer state. The number of final states is given by 



dv 



Vm 



Vd^qf _ 



dE 



f 



(5.13) 



In the typical experimental regime, where the scattering length is still much 
smaller than the size of the trap, a <^ bt, the transition matrix element 
for the non-resonant three-body recombination rate from a cold gas state 
into a shallow dimer state can be estimated perturbatively, substituting the 
asymptotic expressions for the initial and final wave- functions. 



dr dR 



V 



C/(R-^) + C/(R+^ 



gikr giqR 



V vv 



(5.14) 



where r is the distance between two particles, R is the distance between 
their center-of-mass and the third particle, V (x is the normalization 
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volume, k and g (A; ~ g oc 6jr^ <ti a~^) are the initial momenta of the 
cold gas particles, and ipdif) is the s-wave function of the shallow dimer 
with binding energy ti^ /{mo?'). Using the zero-range approximation for the 
transition interaction, U{v) = Uo6{r), Uq = Anh'^a/m and for the dimer 
wave- function, ipdi^) oc exp{—r/a)/{ra^^'^), the matrix element eq. (5.14) in 
the limit k ^ q <^ is estimated as (cf. |FRS96] ) 



Tfi oc 



Y3/2 



m 



(5.15) 



Finally, the non-resonant factor of the three-body recombination rate be- 
comes (cf. |NFJG0ll|FRS96]) 



dn 
~dt 



, ha 



oc n 



3— body 



m 



(5.16) 



where n = Nq/V \s the density. 

For an A^-body recombination reaction into a shallow N — 1 body Efi- 
mov state with the binding energy of the order of h?/{ma'^) the modifica- 
tions to the expression for the rate in eq. (|5.11 ) include the extra 3{N — 3) 
spatial dimensions in the integral in eq. ( 5.14[ ), which gives an extra factor 
(^3/2y-i/2)2(7V-3)^ and also the factor ^"3/6 is substituted by N^/Nl. Thus 
for the A-body recombination rate (A > 3) we have 



dn 
"dt 



N-body 



oc n3 — (na3)^-3 
m 



(5.17) 



We emphasize that these simple estimates only refer to the non-resonant 
contributions and they also do not include other types of decays where the 
final state structure might be substantially different from the shallow A — 1 
cluster plus one particle. The result eq. (5.17) was recently confirmed in 



"09 



Although in the regime na^ <^ 1 the A-body recombination has an 
additional small factor {na^)^~^ it might still be possible to observe the 
A-body Efimov states as a sequence of resonant peaks in the recombination 
rate as function of scattering length, scaling as Sn- 



5.5 Connection to Universal Four-Body States 

Recent theoretical work |vSDG09[ IDvSG09[ lMRD+09| has predicted a uni- 
versal four-body effect where a single atom is loosely bound to an Efimov 
trimer. For each trimer energy two four-body states should exist below. 
This leads to two extra recombination peaks at the critical scattering lengths 
|vSD(;n9j 

045,1 ^ 0.43a3B, aiB,2 ^ 0.90a3B, (5.18) 
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relative to the Efimov trimer resonance a^B (denoted in chapter 2). The 
second state 2 is very close to the Efimov trimer. Two such features were 



reported experimentally in |FKB^09 



The universal scaling factor for this atom- Efimov effect is identical to the 
three-body case, i.e. o,4B~j^/'^4Bj ~ '^'^•'^ ~ since these positions 
simply scale with the three-body states. This is different from our four-body 
Efimov effect which have the scaling factor §4 = 4.29413. However, we note 
that twice the ratio of the two special four body states is 2045, 2/0^45, 1 — 
4.2. This is very close S/^ and probably within the numerical accuracy of 
|vSDG09] . It could indicate an underlying fundamental relationship or be a 
mere coincidence. 

We point out that the two four-body effects are distinct, since different 
degrees of freedom are considered. For both of the effects the four-body 
states are embedded in the complete Efimov spectrum and in this respect 
all four-body states are meta-stable states. To treat these states as bound 
states one must pick out the relevant degrees of freedom, and this is where 
the models differ. The degrees of freedom in the atom-Efimov model is 
obvious, while the four-body Efimov effect presented in this chapter includes 
full two-body correlations. 

Both effects could be relevant for Borromean four-body systems, as also 



pointed out in |MRD"'"09 . The connections and applications need to be 
explored. 



5.6 Conclusions and Outlook 

We have calculated the spectrum of trapped A^-boson systems with = 
3,4, 5, 6, 7 using the stochastic variation method with a correlated Gaussian 
basis. Only two-body correlations were allowed in the variational space. 
Thus the cluster states were a priori mostly excluded from the variation 
space, which made the calculations technically possible. Inclusion of the 
cluster states would turn the Efimov states into meta-stable states. However 
the life-time of these states should be comparable with that of Bose-Einstein 
condensates which have similar structure and decay modes. 

For each system a series of states is found with specific exponential de- 
pendence of the energies and r.m.s hyper-radii on the state number, which is 
a characteristic feature of Efimov states. For the A = 3 system the obtained 
scaling factor agrees well with the known analytical value. We also calcu- 
lated the scaling factors analytically for finite A^ by employing a specific 
Faddeev-Yakubovski decomposition in two-body amplitudes. The results 
are in perfect agreement with the numerical results. 

It might be possible to observe the four-body Efimov states as peaks in 
the recombination rate of a condensate as function of the scattering length, 
similar to the three-body case. For a dilute gas the four-body recombination 
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rate is a factor na^ smaller than the three-body rate, therefore the accuracy 
constraints on the experiment would be higher. On the other hand, the 
large scaling factor 22.7 for three-body recombination is reduced to 4.29 for 
four-body recombination. This would make the universal scaling easier to 
observe. The four-body case may be related to other recently investigated 
atom-Efimov states and could be applicable to Borromean systems. 

We finally note that if the A^-body Efimov states are found, the higher- 
order corrections can be derived by combining the Faddeev-Yakubovski de- 
composition of this chapter with the analysis of chapter |4j This would lead 
to shifts in the critical scattering lengths as discussed in chapter |3| 

In conclusion we have lent theoretical support to the possibility of exis- 
tence of long lived meta-stable A^-body Efimov states in trapped Bose gases. 
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Chapter 6 



Two-Body Correlations in 
Condensates 



6.1 Introduction 

The density of trapped cold gases is generally low under typical experimental 
conditions |PS02] . such that the parameter nrj] is small, where n is the 
particle density and rg is the range of the inter-particle potential. In other 
words the typical distance between particles is much larger than the range of 
the potential, and the typical relative momentum between particles is much 
smaller than the inverse range of the potential. 

In this regime the system of particles exhibits universality (also called 
model-independence or shape independence) |BH06j . The system is not 
sensitive to the details of the potential and the properties of the system are 
essentially determined by only a few low-energy parameters of the potential. 
Very different interaction models then provide quantitatively similar results 
as soon as the low-energy parameters are the same. 

In a two-body system the universality is manifested in the well known 
effective range expansion, where the low-energy s-wave phase shift 5 is deter- 
mined by only two parameters, the scattering length a and the effective range 



Re, see eq. (2.4). The effective range is typically of the order of the range 
of the potential while the scattering length can vary greatly. In three-body 
systems the universality manifests itself in the Thomas collapse |Tho35| , the 
Efimov effect |Efi70] . the Philhps hue |CPI^'G91l iBHvKOOl |FJ02| . and other 
low-energy phenomena |NM99] . also in two dimensions jNF J97] . 

Cold gases also exhibit universality: In the dilute limit their properties, 
in particular the energy per particle, are independent of the shape of the 
inter-particle potential and are determined by the scattering length alone. 
This universality is customarily employed by using the zero-range (pseudo) 
potential for theoretical descriptions of Bose- Einstein condensates. The 
zero-range potential has only one parameter, the scattering length. Com- 
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bined with the Hartree-Fock product wave-function the zero-range potential 
leads to the Gross-Pitaevskii equation |DGPS99l IPS03j . 

The limits of the zero-range model have been tested by numerical calcu- 
lations with finite-range potentials. In particular, repulsive potentials have 
been employed with Monte-Carlo methods |BG01| IGBC99| IDGOl] in the 
case of large positive scattering length. These investigations showed that as 
the scattering length is increased the energy of the Bose gas with a repul- 
sive potential exceeds the zero-range predictions and the condensate fraction 
becomes considerably depleted. 

However, repulsive potentials have a problem when modeling an increas- 
ingly large positive scattering length: The potential range, tq, has to be 
increased essentially linearly with the scattering length. This does not seem 
to match the experimental conditions where the scattering length is adjusted 
by tuning the atomic resonances in an external magnetic field. The range 
of the inter-atomic interaction is then left essentially unchanged. 

Instead, an attractive finite-range potential might be a more realistic 
interaction model for descriptions of trapped Bose gases with Feshbach res- 
onances. Indeed, with an attractive potential an arbitrary large positive 
scattering length can be achieved by fine-tuning the energy of the bound 
two-body state around the threshold, while maintaining the given realistic 
potential range. However, attractive potentials with bound states bring in 
a major complication for numerical calculations: A large number of many- 
body self-bound negative-energy states appears in the system and the con- 
densate state in the trap becomes a highly excited state. 

For a homogeneous Bose gas an approximate Jastrow-type wave- function 
was employed where the pair-correlation function was essentially a solu- 
tion of the two-body equation ^CHM"'"02( IGBC99] . In contrast we propose 
a direct numerical diagonalization of the many-body Hamiltonian where 
the condensate state of trapped bosons appears as a many-body excited 
state which is automatically orthogonal to all the self-bound negative-energy 
states. This proposal opens up the possibility to investigate the correlations 
which are completely absent in standard zero-range mean-field models. 

When correlations are allowed in the BEG wave-function these may in- 
duce a depletion of the fully condensed state. The complete qualitative 
connection between the correlations and the condensate fraction still needs 
to be understood in detail. 

The purpose of this chapter is to investigate the energy, correlations, 
and condensate fraction of a system of trapped bosons. The interaction is 
attractive and supports a single two-body state, i.e. the scattering length is 
positive. We will show that when the scattering length is small compared 
to the trap length the system is model-independent: All potential models 
- attractive, repulsive and zero-range - provide similar results. When the 
scattering length is large the attractive model differs qualitatively from the 
repulsive and zero-range models. In this regime the system with attractive 
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potential becomes independent of the scattering length, with both the energy 
and the condensate fraction converging towards finite constants. The pair- 
correlation function is shown to reflect the two-body interaction at short 
distances and the mean-field solution at larger distances. 

6.2 Condensation and Correlations 

Let us start by introducing the rigorous definitions of condensation and 
correlations. This is done via the one-body density matrix and the two- 
body correlation function. The results apply to the T = formalism used 
in this dissertation. We also comment on some subtleties for condensation 
in attractive systems. 

6.2.1 One-Body Density Matrix and Condensate Fraction 

Assume we have a pure state described by the symmetric many-body wave- 
function ^{ri, . . . ,rjv), which is normalized to A^. The one-body density 
matrix (OBDM) is then defined as 

n'^^\ri,r[) = j ^*{ri,r2, . . . ,rN)'if{r[,r2, . . . ,rN)dr2 . . . dr^- (6.1) 

The diagonal part of the OBDM is the density, n(r) = n^^\r,r), which en- 
sures that the trace equals N. The OBDM is hermitian so we can diagonalize 
it, 

J n(r, r')x^ir') dr' = Na^{r), (6.2) 

obtaining a complete set of orthonormal eigenfunctions Xiif) ^^d corre- 
sponding eigenvalues Aj. The OBDM can then be expanded uniquely as 

n{r,r') = J2NiXt{r)Xi{r'). (6.3) 

i 

The eigenvalues Ni are non-negative, and taking the trace of the OBDM in 
the x-basis gives 

Y^Ni = j n{r,r)dr = N. (6.4) 

i 

This allows us to interpret Ni as the occupation numbers of the single par- 
ticle states or natural orbitals Xi- We can now define a system to be a 
Bose-Einstein condensate if the largest eigenvalue, say Nq, is of order A^, 
while all other eigenvalues are of order 1 |P0561 IYan62j . This means that 
just one single-particle state, xoi is macroscopically occupied. The ratio 
Ao = Nq/N is referred to as the condensate fraction. 
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If more than one single-particle state is macroscopically occupied, the 
system is called a fragmented condensate [Bay 01 IMHUB06] . The factor- 
ization of the OBDM into a single product X*('")x(^') is also referred to as 
perfect first order coherence |NG99j . 

From the OBDM we can get the expectation values of all one-body ob- 
servables, such as position or momentum density. Equation (6.1) can be 
interpreted as the amplitude to stay in the same state, if we destroy a parti- 
cle at r']^ and create one at ri. It can be related to the contrast of interference 
fringes measured in spatial coherence experiments [BHEOOj . 

In practice we calculate the condensate fraction by expressing the OBDM 
in the harmonic oscillator basis and diagonalize numerically. Details can be 
found in appendix [A} 



6.2.2 Condensation in Attractive and Bound Systems 

The definition of condensation introduced above is the conventional one. 
However, [WGS98j considered a gas of weakly attractive bosons in a har- 
monic trap with fixed angular momentum. They found that the angular 
momentum is absorbed by the center-of-mass motion, and that the system 
is fragmented, i.e. non-condensed, according to the definition used above. 
This is quite surprising, because the system is just a Bose-Einstein con- 
densed ground state performing a center-of-mass motion. Ref. |PPOO| then 
proposed an alternative definition in terms of an internal density matrix, 
where the center-of-mass is integrated out. According to this definition the 
system will be a pure condensate. However, ||Gaj06j argued that the internal 
density matrix is only useful when describing observables that are indepen- 
dent of the center-of-mass. It is not useful when measuring global properties 
of a gas such as condensation. 

These subtleties are easily explained in terms of conditional measure- 
ments. In a perfect BEC two subsequent measurements of individual parti- 
cles are completely uncorrelated. The spatial probability distribution for the 
first particle is determined by the mean-field density, and so is the second. 
However if all the particles move around in a correlated way, e.g. because 
the system is bound, a measurement of the first particle will determine the 
systems center-of-mass with high accuracy. A measurement of a second 
particle will be close to the center-of-mass, and condensation is lost. 

A concrete example on this is given in appendix |B] where we show that 
strong hyper- radial correlation can deplete the condensate completely if the 
internal length scale is very different from the center-of-mass scale. Also see 
appendix |B] for a discussion on how other types of correlation affects the 
condensate fraction. 

In the light of all these subtleties we have chosen to use the conventional 
criterion for Bose-Einstein condensation. 
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6.2.3 Two-body Correlation Function 

The two-body correlation function is defined as 

n(^)(ri, r2) = j ^* {ri,r2, . . . , r7v)^'(ri, r2, . . . , tn) dr^ . . . dr^, (6.5) 

and is interpreted as the joint probabihty of detecting two particles at ri 
and r2 simultaneously. If the particles are totally uncorrelated the two-body 
correlation function factorizes as n^'^\ri,r2) = n(ri)n(r2). The normalized 
two-body correlation functioi|^ is defined as 

Ci^>(rur.) = "^p^. (6.6) 

n{ri)n{r2) 

It describes the tendency of particles to cluster (C*^^-* > 1) or to stay sepa- 
rated (0 < C^^) < 1). If no two-body correlations are present, C*^^) equals 1 
identically. 



In a homogeneous gas the correlation function (6.6) is translational in- 
variant and depends only on ri2 = — T2I. In a spherically symmetric trap 
it also depends on the lengths ri and r2. However, for relative distances ri2 
much smaller than the density variations and the trap length, the correlation 
function will only depend on ri2. Thus we will typically take r2 = 0. 

The correlation function in quantum gases can be measured via noise 
correlations |ADL041 IBDZOSj . The experimental procedure is to release 
the gas from the trap and let it expand. The expanded density profile is 
then measured by light absorption imaging, and correlations in the noise 
are extracted. It is preferable that the gas is weakly interacting during the 
expansion. For Feshbach resonance experiments this is usually obtained by 
also ramping the scattering length down to zero when the trap is turned off. 

The technique of noise correlations have recently been used to measure 
the correlations for fermionic gases and for bosons in optical lattices |ADL041 
IBDZ08] . but it also works for normal trapped BECs. 

6.3 System and Numerical Techniques 

We consider a system of N identical bosons with mass m and coordinates 
Ti, i = 1, . . . , N , in a spherical harmonic trap with frequency uj. The Hamil- 



tonian of the system is given by eq. (2.70), where the system parameters are 
taken from [BGOlj : m = 86.909amiipRb) and to = 2it x 77.87Hz, giving 
the trap length bt = (muj) = 23095au=1.22/im. 



The zero-range model potential, eq. (2.59), has the scattering length a 



as the only length parameter. For dilute bosonic systems this parameter is 



^Simply referred to as the (two-body) correlation function below. Note that several 
different (but equivalent) definitions and notations are used in the literature. 
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customarily chosen to be equal to the inter-atomic scattering length. The 
zero-range potential provides the correct low-energy scattering amplitude in 
the first order Born approximation. The zero-range potential can only be 
used with an appropriate non- correlated functional space, like the Hartree- 
Fock product wave- functions. Using the zero-range potential in a correlated 
space will lead to a complete Thomas collapse of the system |FJ01a] . 



For the finite-range potential model we use a Gaussian, eq. (2.20), with 
the range ro = 11.65au and a varied negative strength Vq. The variation of 
the strength is limited to the region where the potential provides exactly one 
two-body bound state and has a positive scattering length. The range of the 
interaction is much smaller than the trap length, i.e. ro/bf = 5.04 x 10~^. 

Our numerical procedure is the stochastic variational method described 
in chapter [2] The wave-function of the system is represented as a linear 
combination of basis-functions taken in the form of symmetrized correlated 



(k) A I 

Gaussians, eq. (2.74). The non-linear parameters • in eq. (2.74) are op- 



timized stochastically by random sampling from a region that covers the 
distances from ro to bt. The center-of-mass motion is assumed to be in the 
oscillator's ground state. 



The zero-range potential eq. (2.59) requires an uncorrelated variational 
wave- function, and we use the hyper-radial wave- function in eq. (2.75). 



The zero-range potential with this variational space provides results similar 
to the Gross-Pitaevskii equation ^SFJ05| . 

The calculation of a highly excited state with the fully correlated basis 
^fuU) eq. (2.74), is a difficult numerical task and is only possible for rela- 
tively small number of particles. However, for a typical system of trapped 
atoms even when the scattering length is large the density of the system 
remains small, nrj] <C 1, and one can assume that only binary collisions 
play a significant role in the system's dynamics. In this approximation the 
variational wave-function can be simplified by only allowing two-body cor- 



relations in the basis-functions as in "^23^ eq. (2.77). The symmetrization of 



this function can be done analytically |SFJ05j which greatly simplifies the 
numerical calculations. 



During the calculation of a given system the number of Gaussians in the 
basis is increased and the stochastic optimization is carried out until the 
number of negative energy states and the energy of the lowest state with 
positive energy is converged. The convergence within four digits typically 
requires about 5 x 10^ Gaussian and about 10^ random trials per nonlinear 
parameter. When the energies are converged, the relevant properties of the 
wave- function can be extracted. We choose to consider densities, condensate 
fractions, and two-body correlation functions. 
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6.4 Results 

6.4.1 Obtaining the BEC State 

For repulsive potential models it is simply the ground state of the trapped 
many-boson system that is identified as the Bose-Einstein condensate (BEC) 
state. With deep attractive two-body potentials, however, the many-body 
system in a trap has a large number of self-bound negative-energy states 
and identification of the BEC-state is not obvious. 

The typical spectrum of a trapped many-body system with attractive 
potentials and including two-body correlations is shown in fig. 6.1 The 



system has a number of deeply bound states with negative energies and then 
a positive quasi-continuum spectrum which starts at about {2>N/2)huj and 
has the characteristic distance between levels of the order huj <^ / (2mro)^. 
Apparently the BEC-state should then be the lowest state of the quasi- 
continuum spectrum or, equivalently, the lowest state with positive energy. 



3 




Figure 6.1: A typical spectrum (in the vicinity of zero energy) of a system of 
A'^ = 20 bosons in an oscillator trap, eq. (2.701, interacting via an attractive two- 



body potential, eq. (2.201, with one bound state and a positive scattering length. 



The inset shows the beginning of the so-called quasi-continuum spectrum. 



To verify this conjecture we calculate the central density, no, of the 
system for the negative- and positive-energy states around zero energy. The 



results are shown in fig. 6.2 in the form of the inverse central density (the 



volume per particle) rig . In the BEC-state the atoms should occupy the 
whole volume of the trap and thus the volume per particle should be close 
to one (in the correspondingly scaled oscillator units). The states with 
negative energy are self-bound states with much higher density and thus 



much smaller volume per particle. And indeed that is what fig. 6.2 shows - 
a sharp increase in the volume per particle from a small value to about unity 
exactly at the lowest state with positive energy where the quasi-continuum 
starts. 
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Figure 6.2: The inverse central density ^ (in oscillator units nose = 7r~^^'^Nb. 
for a system of N bosons in an oscillator trap, eq. (2.70 1, interacting via an attrac- 
tive two-body potential, eq. (2.20), with scattering length a = 119au as function of 
the state number. The lowest state with positive energy is numbered zero. 




74 



CHAPTER 6. TWO-BODY CORRELATIONS IN CONDENSATES 



Another test is the condensate fraction, shown for several states around 



zero energy in fig. 6.3 The self-bound states with negative energy must have 
smaller condensate fraction compared to the BEC-state, and the excitations 
from the BEC-state must gradually deplete the condensate fraction. Ap- 



parently this is what is seen in fig. 6.3 - a sharp increase of the condensate 



fraction to about 100% at the lowest state with positive energy with the 
subsequent gradual depletion. 

We have thus verified that in the case of attractive potentials the BEC- 
state of a system of trapped bosons is the lowest state with positive energy. 

6.4.2 Accuracy of the Two-Body Correlated Basis 

In |BG01) the energies of several low-density systems of trapped bosons were 
calculated using the Gross-Pitaevskii equation as well as repulsive hard- 
sphere models, both with the same scattering length of 100 au. In this 
regime the systems exhibit universality and the energies calculated in both 
models were very close. 



To test the accuracy of our two-body correlated basis, eq. (2.77), (which 
is expected to be a good approximation in the low-density regime) as well 
as the identification of the BEC-state for attractive potentials we consider 
the same systems with the same scattering length but with the attractive 



Gaussian potential, eq. (2.20), and calculate the energy of the BEC-state 
according to our prescription: The BEC-state is now an excited state and 
is identified in the calculations as the lowest state with positive energy. We 



also calculate the energies for the zero-range potential model, eq. (2.59) 



with the hyper-radial trial wave- function, eq. (2.75). 



The results are given in tab. 6.1 As expected, these low-density systems 
with relatively short scattering length seem to exhibit universality as all the 
potential models give essentially the same results. We conclude that our 
identification of the BEC-state is correct and that the two-body correlated 
basis has an adequate accuracy. 

To check the accuracy of the two-body correlated basis also for large 
scattering lengths we perform a test calculation for = 4 particles with fully 
correlated and with two-body correlated basis for vastly different scattering 



lengths. The results are given in tab. 6.2 Although the accuracy of the two- 



body correlated basis decreases somewhat with the increase of the scattering 
length, the relative accuracy is better than 1% even for exceedingly large 
scattering lengths. 

6.4.3 BEC Energies 

We first show the zero-range results for = 10, 100, 1000 in fig. |6.4[ The 
energy depends on N only through the combination {N — l)a/bt and is in 
perfect agreement with the Gross-Pitaevskii result. In the large (A — l)a/bt 
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N 


GPt 


HSt 


ZR 


A 


3 


4.51032 


4.51036(2) 


4.5103 


4.510 


5 


7.53432 


7.53443(4) 


7.5342 


7.534 


10 


15.1534 


15.1537(2) 


15.1533 


15.154 


20 


30.638 


30.640(1) 


30.6394 


30.640 



t Data taken from fBGOT| . 

Table 6.1: The energies in units of Hlxj for the BEC-state of a system of A'' bosons in 



a harmonic trap, eq. ( 2.70 1 , for different interaction models with the same scattering 
length a — lOOau. For the Gross-Pitaevskii (GP), hard-spheres (HS), and zero- 
range (ZR) models the BEC-state is the ground state. For the attractive model 
(A) the BEC-state is the lowest state with positive energy. The zero-range model 



the two-body correlated basis, eq. (2.771 



employed the hyper-radial basis, eq. (2.75), while the attractive model employed 



Vo X 10^, au o, au 


E{2B) E(full) 


-1.400 119.4 


6.025 6.025 


-1.300 327.0 


6.067 6.067 


-1.290 402.4 


6.083 6.083 


-1.280 525.4 


6.108 6.108 


-1.270 761.0 


6.155 6.156 


-1.260 1400 


6.282 6.283 


-1.255 2430 


6.478 6.481 


-1.252 4370 


6.818 6.848 


-1.251 5962 


7.059 7.112 



Table 6.2: The energies, in units of hu, of the lowest state with positive energy for 
a system of iV = 4 bosons in a harmonic trap, eq. (2.70). The bosons interact via 
an attractive Gaussian potential, eq. (2.20), with the strength Vq and the scattering 
length a. The results from the fully correlated basis, eq. ( |2.74 1, and from the two- 
body correlated basis, eq. (2.77), are designated i?(full) and i?(2B), respectively. 
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limit we reproduce the energy of the Thomas- Fermi hmit |PS02j ^ 



( 



15Na 



14 V h 



2/5 



HllJ. 



(6.7) 



For smah {N —l)a/ht the interaction can be treated perturbatively: With the 
zero-range perturbation H' = '^i^j{^T^fi^a/m)5{rij) on the non-interacting 
single-particle product ground state we find the energy shift 



AS 

Nhuj 



1 {N-l)a 



'2-K 



(6.8) 



Our results also reproduces this regime. The agreement of our results with 
the GP equation was expected since both models are based on the same 
zero-range interaction. Also, the hyper-radial variational space for the zero- 
range interaction is reminiscent of the Hartree-Fock product space used in 
the GP equation. 
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Figure 6.4: The energy per particle E/N as function of the scattering length a 



for a BEC-state of a system of N identical bosons in a harmonic trap eq. (2.701 



for the zero-range model eq. (2.591 (points). The results are TV-independent in this 
parametrization and in perfect agreement with the Gross-Pitaevskii (GP) equation. 
Also shown is the perturbative and Thomas-Fermi (TF) limits, see text. 



We now consider the energy per particle for the attractive potential 



model in fig. 6.5 and compare with the zero-range model. For small scatter- 
ing lengths the different models give the same universal results - the system 
is model-independent. For larger scattering lengths the energies from the 
attractive model are systematically below the zero-range model, quite unlike 



■^See also chapter [s] for a complete discussion on the Thomas- Fermi approximation. 



77 



6.4. RESULTS 



the repulsive model which goes above the zero-range model [BGOIJ IDGOl] . 
For very large scattering lengths the attractive model, unlike the zero-range 
and repulsive models, becomes insensitive to the scattering length and the 
energies converge to a constant. This is consistent with the Jastrow-type 
approximation used in |CHM^02 . 



in 




0.01 



{N - l)a/bt 

Figure 6. 5: The energy per particle E /N as function of the scattering length a for a 
BEC-state of a system of N trapped bosons. The results for the attractive Gaussian 
interaction (points) is compared with the zero-range/Gross-Pitaevskii result from 



fig. 6.4 



In the regime close to the two-body threshold, where the scattering 
length is large, an arbitrarily large change in the scattering length needs 
only an infinitesimally small change of the depth of the attractive poten- 



tial (see tab. 6.2). Therefore, when the scattering length is larger than the 
trap length it ceases to be a physical length scale for the system. Since the 
external oscillator potential turns all continuum states into discrete states 
all singularities due to various thresholds are removed. Clearly an infinitesi- 
mally small change in the potential, despite the large change in the scattering 
length, only leads to an infinitesimal linear change in the energy which then 
becomes independent of the scattering length. 

On the contrary, for finite-range repulsive potential models an increase of 
the scattering length needs an almost proportional increase in the potential 
range. Thus the system never ceases to depend on the scattering length. 

If we plot the data using a different parametrization, namely E as func- 

1 /2 

tion of {N -l){a/bt)'', the energy data points seem to follow a universal 
curve as shown in fig. |6.6[ The specific value to which the energy converges 
probably depends on an independent energy scale, as for the Efimov effect 
in chapter [3] and [5} We speculate that other non-Gaussian interactions will 
reproduce the same universal curve, by tuning only a single parameter. 
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6.4.4 Condensate Fractions 



Our results for the condensate fraction of a system of trapped identical 
bosons in a BEC-state are shown in fig. |6.7[ For small scattering lengths 
the system is 100% condensate. When (A^ — 1) {a/htf'^'^ is about 3 the 
condensate fraction rapidly drops a few percent before stabilizing again. 
Although we could not reach further due to numerical convergence problems, 
we expect that the condensate fraction will not appreciably change with 
further increase of the scattering length. This seems consistent with the 
energies being stagnated when approaching the two-body threshold. 




1 2 3 4 5 6 

(7V-l)(a/&0* 



Figure 6.7: The condensate fraction of a BEC-state of a system of N identi- 
cal bosons in a harmonic trap eq. (2.701 with attractive inter-particle potential 
eq. (p:20| as function of {N - \){a/bt)^. 
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This behavior is quahtatively different from what happens within the 
repulsive models. In the hard-sphere Monte-Carlo simulations |DG01] . the 
condensate fraction starts to deviate from 100% much later when a/bt is 
of the order of 0.1, and with further increase of the scattering length the 
condensate is completely quenched. Again our results are different from the 
Jastrow-type approximation |CHM"'"02 , where the condensate also becomes 



fully quenched. However, the authors of |CHM^02 note that their estimates 



for the condensate fraction may be rather crude. 

We have to note also that our wave-function includes only two-body 
correlations which may lead to an over-estimate of the condensate fraction. 

6.4.5 Two-Body Correlation Functions 

We now discuss the results for the two-body correlation function C^^^(ri, r2). 



eq. ( 6.6 ) . For simplicity one of the coordinates is set at the center of the trap, 
r2 = 0. This measures the probability to find a second particle at radius r 
when the first one is at the center. In this case the two-body correlations 
depend on one coordinate only, and we can easily investigate the properties 
over many length scales, including vq and ht- However, the main conclusions 
will of course hold for arbitrary r2- 

Figure 6.8 shows the two-body correlation function C^^^(r, 0) for N = 
20 bosons in the case of an attractive Gaussian interaction with one two- 
body bound state (a > 0) and two-body correlated functional space. Three 
different states are shown: The BEC-state determined by the lowest state 
in the quasi-continuum, the first state (BEG— 1) right below the BEG-state, 
and the ground state. 

Let us first consider the asymptotics: The correlation function for the 
BEG-state approaches unity for r — > co, thus the particles are completely 
uncorrelated at very large distances. This is what we expect, since we have 
shown earlier that the BEG-state has condensate fraction close to unity. The 
BEG— 1 state is self-bound (bound without the trap) and hence the particles 
are strongly correlated. The correlation function goes to zero asymptotically, 
but is very large within the potential range rg. If one particle is measured 
at a particular point in space it is very likely that all other particles are 
close to that point. Thus the state is clearly non-condensed as discussed in 
subsection I6.2.2[ 

Let us only consider the BEG-state in the rest of this chapter. We 



observe that the two-body correlation function in fig. 6.8 is large and positive 
within the range ro of the potential. Outside ro we observe a large dip in 
the correlation function. The value at the minimum of this dip is very small 
but finite. The position of the dip occurs exactly at the scattering length, 
r = a. 

To understand this we consider the two-particle problem without a trap. 
If the two-body potential supports a single two-body bound state the relative 
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o 



ground state 




0.0001 



0.001 



0.01 

r/bt 



0.1 



Figure 6.8: The two-body correlation function C^^-'(r, 0) for some of the different 
states found with the stochastic variational method. The correlation for the BEC- 
state (lowest state in the quasi-continuum) reflects the two-body interaction at 
short distances and goes to unity at large distances. The first state below the 
BEC-state (BEC— 1) is self-bound and the correlation function approaches zero at 
large distances. This is even more pronounced for the ground state (GS). 



wave- function R{r) has a single node. If a is larger than rg, the position 
of this node is approximately at a. So the joint probability, and hence the 
two-body correlation function, is zero at the scattering length. 

To investigate this in more detail we take a new two-body potential with 
two bound states and a large a > 0. We calculate both the radial function 
R{r) for the two-body problem and the correlation function for = 20 



particles. The result is shown in fig. 6.9 



At small distances the correlation function closely follows the two-body 
wave-function while at large distances it converges to unity. This behavior 
agree well with the Jastrow product approximation used in [CHM^02]. We 
conclude that the correlation function is completely model-dependent for 
distances ri — r2 comparable to the potential range. However, for larger 
distances the scattering length determines the properties. 

Let us now investigate the large scattering length limit, where the details 
are independent of the short-range potential. In fig. 6.10 the correlation 
function is plotted for two different scattering lengths. The ai data are the 
same as the BEC-state in fig. |6.8[ When the scattering length is increased the 
dip moves correspondingly outwards but becomes shallower. The decrease 
of the dip happens since the correlated pairs starts to feel the mean-field 
of other particles between them. Thus the probability to find two particles 
within a distance a is increased considerably. This is easily understood in 
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ro 



o 



0.0001 



N — 2 radial func. 
= 20 corr. func 




Figure 6.9: Correlation function (points) for the BEC-state of = 20 trapped 
bosons with scattering length a and two bound two-body states, compared with 
the squared two-body zero-energy radial function R{r)^ (line). R is normalized to 
to match the correlation function in r = 



o 




r/bt 

Figure 6.10: Correlation function for = 20 bosons in a trap. An attractive 
Gaussian interaction (range tq — 11.65au) was used in combination with the two- 
body correlated functional space. The bosons are in the BEC-state, identified as 
the lowest state of the quasi-continuum. Two calculations with different scattering 
lengths (fli and 02) are shown. When the scattering length is small compared to 
the trap length ht a dip occurs at a. When a increases the dip moves outwards and 
becomes shallower. 
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terms of the universal shallow (Feshbach) dimers with sizes of order a. We 
conclude that in the large scattering length limit the two-body correlation 
function has a universal behavior. We propose to measure this effect via 
noise correlations as discussed earlier in this chapter. 

6.5 Conclusions and Outlook 

We have calculated the energy and the condensate fraction of a system 
of N bosons in a harmonic trap as function of the number of bosons and 
the scattering length a. Specifically, we considered the regime where the 
scattering length is positive and comparable to the trap length. The positive 
scattering length is modeled using an attractive two-body potential with a 
bound two-body state. The many-body system then has a large number of 
negative-energy self-bound states and the condensate in the trap is identified 
as the lowest excited state with positive energy. 

When the scattering length is small compared to the trap length the 
system shows model-independence (universality) - the results from the at- 
tractive potential model are very close to those from the zero-range and the 
repulsive potential models. 

In the limit of large scattering length the system properties become inde- 
pendent of the scattering length, contrary to the zero-range and the repulsive 
models. For the attractive potentials the energy per particle of the system 
of trapped bosons tends to follow a universal curve. The condensate fraction 
decreases with the scattering length and reaches a finite constant at large 
scattering lengths contrary to the repulsive models where it reaches zero in 
this limit. 

The two-body correlation function reflects the finite-range interaction 
a short distances and the mean-field solution at larger distances. When 
the scattering length becomes large the two-body correlations become long- 
ranged and universal. 

In general, the stochastic variational method allows for calculations of 
correlations in BECs with attractive potentials. Such calculations have tra- 
ditionally not been possible because of large amount of low-lying states. The 
obtained BEG states contain the short-range correlations while the normal 
mean-field features are retained. Thus, this is an important step towards 
understanding the non- universal corrections in condensates. 
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Chapter 7 



Mean-Field BEC with 
Higher-Order Interactions 



7.1 Introduction 



The stability of Bose-Einstein condensates (BECs) in ultra-cold alkali gases 
is determined by the sign of the scattering length a |DGPS99] . For a < 0, 
one has effectively attractive interactions and the condensate will collapse 
to a dense state when the number of condensed particles, A^, is larger than 
a critical number |RHBE95l IBP961 IDGPS991 IGFTOl] . This has been 
beautifull y demonstrated in experiments with ^Li |BSH97j, ^ '''Rb |RCC+Oli 
iDCC+ni) . and recently with a dipolar ^^Cr BEC |LMF+n8) . The findings 



indicate that the theory based on the surprisingly simple Gross-Pitaevskii 
(GP) equation can reproduce and describe most features of the experiments, 
thus the scattering length is a universal parameter for these many-particle 
systems. 

The GP equation includes two-body terms through a contact interaction 
which is parametrized by a. This is equivalent to a Born approximation 
but with an effective coupling that is obtained by replacing afeorn by the 
physical scattering length a. However, higher-order terms in the expansion 
of the phase shifts at low momenta, determined by the effective range Re, 
the shape parameter etc., give corrections to the simple GP equation. In this 
chapter we explore the influence of the effective range term on the quantum 
properties of a BEC. In particular, we show that the critical number of 
condensed atoms depends strongly on the higher-order scattering term when 
the scattering length approaches zero (zero-crossing). We also show how 
the macroscopic quantum tunneling (MQT) rate, in which the entire BEC 
tunnels as a coherent entity, can be modified for small condensate samples. 

The considered effects depend on a combination of a and Re which yield 
different behavior for wide and narrow Feshbach resonances. Recent mea- 
surements on ^^K found many both wide and narrow resonances which allow 
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for tuning of a over many orders of mag nitude [ RZCD+O?! [DZF+n7| . We 
therefore consider a selected example from '^^K in order to elucidate the 
general behavior for realistic experimental parameters. 

In this chapter we first introduce the modified GP equation with effective 
range dependence. Using both a variational and numerical approach we find 
a phase diagram describing the stability of the condensate. We then consider 
an extended Feshbach resonance model including effective range variations. 
The behavior of the critical particle number near a scattering length zero- 
crossing is derived. We discuss MQT and show numerically how the rate is 
modified. We finally discuss other possibilities for probing the higher-order 
interactions. 



7.2 Modified Gross-Pitaevskii Equation 



We assume that the condensate can be described by the GP equation and 
we focus on the a < attractively interacting case. Since we are interested 
in the ultra-cold regime, where the temperature is much smaller than the 
critical temperature for condensation, we adopt the T = formalism. In 
order to include higher-order effects in the two-body scattering dynamics, we 
use the modified GP equation derived in jCMPOT] for which the equivalent 
energy functional is 



dr 



It 
2m 



V^\^ + Vecct{r)\^\^ + 



Uo 



(7.1) 

where m is the atomic mass, Vext is the external trap, Uq = ATrh^a/m, 
and (72 = — aRe/2 with a and Re being the s-wave scattering length 
and effective range, respectively. The corresponding stationary GP equation 
found by variation of ^ is 



2m 



V2 + 14xt(r) + C/o(|^P + 92V2|^| 



(7.2) 



where the chemical potential, was introduced to fix the particle number 
N. 



7.2.1 Effective Zero-Range Interaction 

To derive the modified GP equation we first need an effective two-body 
zero-range interaction, which by construction will give the same interaction 
energy as the real two-body potential. Partly following |RFOH ICMP07j , we 
first consider two non-interacting particles and restrict ourselves to s- waves 
only. In order to count energy levels and evaluate the corresponding energy 
shifts we choose the radial wave-function to be zero some large but arbitrary 
radius L. This gives the discrete energy levels En = h?k'^/m, kn = nn/L 
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and normalized eigenstates ipnii^) = Y^2/(47rL) sin(/c„r) /r. Inclusion of the 
true two-body interaction gives the new energy levels = fi?K^/m and 
corresponding small energy shift A£'„ = En — En- The criterion 



AE„ 



{lpn\yZR\lpn) 



(7.3) 



then defines the coupling constant (s) of the effective zero-range interaction 

VZR. 

We first evaluate AEn- In case of interaction the asymptotic wave- 
function is proportional to sin(knr + 6(kn)), where 6 is the phase shift. Since 
the wave-function must still vanish at L, the wave numbers are given bj|^ 
knL = nvr — 5(kn)- Since Akn = kn — is small we have 5{kn) — S{kn). 
Using AEn — 2h^knAkn/m we find 



5{kn) 2En 



L 



(7.4) 



Since we are only interested in terms up to order /c^, we use the low-energy 
effective range expansion 



k cot 6 



-- + ^e + o{k') 

a 2 



defining the scattering length a and effective range Re- Using 5 
[tan 6)^/3 we find 



m 
k 



1 



^ 3 



aRp 



)k^ + 0{k^ 



(7.5) 

tan 6 — 

(7.6) 



(7.7) 



We now construct the zero-range interaction. The simplest hermitian 
ansatz up to order k'^ is Vzr = Vq + V2 with 

Voir) = UoS{r), 

Mr) = Uog20l-S{r) + 6{r)Vl.], 

where Uq and §2 are unknown coupling strengths to be determined from 
eq. (7.3). Using the non-interacting wave-functions ipn we get 

Uom 



{i'n\VzR\ll'n 



(1 - g2k 



2En 

L ■ 



(7.8) 



Comparison of eqs. (7.4), (7.6), and (7.8) gives Uq = AnhPa/m and 52 = 
a^/3 — aRe/2. Thus we now have an effective zero-range two-body inter- 
action to order k^. By constructing the many-body mean- field Hamilto- 
nian from this interaction (integrating out relative pair-distances), the last 
higher-order term of eq. (7.1) can be obtained |.FWG031 ICMP07) . 



^If the two-body interaction supports nt bound states the phase shift must be reduced 
by ntn according to Levinsons theorem. 
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7.3 Variational Calculations 

We are interested in the stability properties of the ground-state and we 



therefore perform a variational calculation on eq. (7.1 ) using the mean- field 
trial wave-function 



exp(- 



VT' 



(7.9) 



where q is the dimensionless variational parameter and bt = y^h/moj is the 
trap length. The normalization is N = f dr|^'(r)p. For simplicity we only 
consider isotropic traps with Vextir) = ^mcu^r^. However, the effects found 
should hold for deformed traps as well (along the lines of the analysis in 
pTOSj). The variational energy is 



Nhw 



2 3 1 

' +7^ + 
4 



1 iV|a| / 1 



2tt bt 



_|_ 3£f 

j3 ^2 g5 



(7.10) 



In fig. 7.1 we plot E{q) for different parameters. As shown, there are many 
possibilities for 52 / 0, including stable, unstable, and metastable systems. 
We see that the 52 term modifies the barrier for N\a\/bt = 0.5, implying 
that tunneling rates will be altered. For N\a\/bt = 0.7, the 52 = case has 
no barrier at all, and here the §2 term can in fact produce a small barrier 
on its own. The q~^ dependence of the term means that the effect is small. 
However, the plot clearly shows that a new stability analysis is needed. In 
addition to the variational approach we have numerically solved the full 



time- independent GP equation corresponding to eq. (7.1). 
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Figure 7.1: Energy of a BEC with fixed N\a\/bt as function of the variational 
parameter q, i.e. the size of the condensate. The higher-order interaction term 92 
modifies the height and shape of the barrier. 



88 



CHAPTER 7. MEAN-FIELD EEC WITH HIGHER-ORDER 
INTERACTIONS 



7.4 Phase Stability Diagram 



To determine the ground-state stability one looks for the vanishing of the 
barrier towards q = 0. For 52 = eq. (7.10) leads to Nc\a\/bt ~ 0.671 
|DGPS99] . The full integration of the CP equation gives Nc\a\/bt « ko, 
ko = 0.5746 |RHBE95llGFT01j . These values are indicated by filled points in 



fig. 7.2 The stability coefficient has been determined precisely to 0.547(58) 
using bound-state spectroscopy near a Rb Feshbach resonance jCKT+03 



The deviation is well understood in terms of the asymmetric trap used 

|(;ftoi] . 
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Figure 7.2: Phase stability diagram of a BEG with higher-order interactions. The 



solid line obtained from the variational ansatz, eq. (7.9 1, divides the stable (S), 



metastable (M), and unstable (U) regions. The white region indicates where sta- 
tionary solutions exists for the full GP equation. Filled points show the well-known 
32 = results. Arrows correspond to the values in fig. |7.1[ The dashed line is the 
fit in eq. (7.171. Two specific "^^K values are chosen for zero-crossing (cross) and 



MQT rates (open square) calculations (see text). 



In the general 32 7^ case we first take the variational energy eq. (7.10) 



and solve for multiple roots of dE{q)/ dq. The resulting "phase diagram" is 
plotted in fig. |7.2| (solid line). In the upper left region (S) we have complete 
stability of the condensate: only one minimum exists at large g ~ 0.8 (see 



fig. 7.1 ) and the potential goes to plus infinity at small q. In the metastable 



region (M) a barrier and a minimum exists for large q: for (72 > another 
minimum exists at small q < 0.2, while for g'2 < the potential goes to minus 
infinity. In the unstable region (U) the barrier has vanished: the potential 
either has a minimum at small q (§2 > 0) or no minimum at all ((72 < 0). In 
the variational approach the stable and unstable regions are connected via 
the upper right part of fig. 7.2, Going from (S) to (U) corresponds to an 
adiabatic change, where the macroscopic {q ^ 1) minimum is transferred to 
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a microscopic (g ^ 1) high-density minimum. 

Next, we numerically solve the full the GP equation. The stationary 



(white) and non-stationary (shaded) regions are shown in fig. 7.2 Since 



both stable and metastable solutions are considered stationary, the white 
region covers both (S) and (M). Both our variational and numerical results 
agrees with the known g'2 = results, and have similar behavior for small 
and negative 92 /bl- 



7.5 Feshbach Resonance Model 

In order to predict effects of the higher-order term, we need a realistic model 
for g2. Since 52/^4 is the relevant parameter, and the trap length, bt, is 
usually orders of magnitude larger than atomic scales, it is necessary to look 
for divergences of (72- Since §2 depends on a and Re-, a Feshbach resonance 
is the obvious mean. The standard single-channel models are inadequate 
since only a is considered. 

We therefore use a multi-channel model |BJK05] , which describes both a 
and Re as a function of resonance position Bq, width AS, magnetic moment 
difference between the channels A//, and the background scattering length 
aiyg. In |BJK05j the on-shell T-matrix in the open-open channel is given by 



T 

00 



abg 



y- + fi^k'^/m-Afi{B-Bo) ) + '^'^^9'^ 
Airh'^ 1 



(7.11) 



where k is the relative momentum of the atoms and we introduced R, 



eO 



2h'^/mabgA^AB < from eq. (2.14) and r] = {B - Bq)/AB. Expanding 



the denominator of eq. (7.11 ) to second order in k we get 



By comparing this to the effective range expansion of the vacuum T-matrix, 
on-shell T-matrix jPSn2] 

A-Kh'^ 1 47r^2 I 

X = . ~ . (7.13) 

m k cot 6{k)—ik m 1 _ -|- iA: ' 

the scattering length and effective range become^ 

a{B) = at,g (1 - , (7.14) 



^We note that eq. ( |7.15| l can also be obtained from [Jon04l eqs. (16) and (44)]. 
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In |BJK05] only the regime \B — Bq\ <^ AB was considered where Re ~ Reo 
is constant. The modification above includes the field dependence of Re and 
since we later want to consider a broad region around Bq we retain this 
explicit field dependence. 

We note that Re can become large (and negative) in two different ways. 
The first way is to consider narrow resonances where AB is small, and hence 
Reo is large. The other option is to tune the scattering length near a zero- 
crossing to small values, a ^ oq. The effective range goes to zero far away 
from the resonance where a ~ abg. 

We now have a field-dependent g2{B) for given values of i?o, AS, A//, 
and Ofeg. Notice that with this Feshbach model 

aReo ( ai,g\2 
52(a) = y-^(l--j , (7.16) 



hence 52 diverges when a — > (referred to as zero-crossing) or a — > 00 
(on resonance). These are well-known features for model potentials like the 
square- well and van der Waals interactions. 

As a concrete example we use the extremely narrow Feshbach resonance 
in found at B p = 825G, with AB = -0.032G, Afi = 3.92^iB, and atg = 
— 36ao |DZF"'"07| . For this resonance we find a large Reo = —2.93 x IO^oq. 



The variation of a. Re and §2 as function of B is shown in fig. 7.3 We use 
trap length bt = 1.84^m = 3.48 x lO^ao in all calculations unless indicated 
otherwise. 
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Figure 7.3: Scattering length a, effective range Rf,, and coupling constant 32 as 
function of B-field for the narrow ^^K Feshbach resonance at Bq = 825 G. The 
inset shows the MQT rate T2 (and Tq for 52 = 0) for N = 242 and bt = 1.84^m. 



The cross and open square are as in fig. 7.2 
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7.6 Critical Particle Number Near Zero-Crossings 



As noted above, large values of §2 are possible at zero-crossings (a = 0) 
since Re diverges here, see e.g. fig. |7.3| around B — Bq — 32mG. The 
scattering length is currently being tuned with extreme accuracy near such 
zero-crossings m ^Li |PD.T+n9l . 

The advantage is that many particles Nc oc l/\a\ can be accommodated 



in the condensate. However, for a — > we have ag2 



-Reoio-bg) /2, 



I.e. a 



finite limit. Remembering that Reo and V |^'| are negative, the last term in 



the energy functional eq. (7.1) also becomes negative. Thus larger densities 



or density fluctuations gives lower total energy. This implies less stability 
and smaller Nc near a = 0. 

To calculate quantitative effects on Nc we focus on the critical line be- 



tween the unstable and metastable regions in fig. 7.2 for §2 < 0. We fit the 
dependency as 

-1 



Nc = k 



01 



bt 



1 - ki 



92[a) 



(7.17) 



where k\ = 14.5, as shown in fig. 7.2 With 52(0) given by eq. (7.16), Nc 



7.4 



becomes a function of |a|/6t for fixed values of a^g/bt and Rco/bt- In fig. 
we plot 1/A^c as function of \a\/bt for the ^^K resonance with various trap 
lengths. The curve deviates from the linear g'2 = result both at \a\ ~ 
and \a\ ~ bt- Squeezing the trap makes the effect even more pronounced. In 
the limit \a\/bt the critical particle number approaches a finite value. 



2kob'l /{ki\Reo\alg). Large effects also occur near \a\ ~ bt, however. 



here the allowed particle number is very small. 




10-^ 10-^ 10-3 10-2 10-1 
W\/bt 

Figure 7.4: Critical particle number Nc as function of scattering length a. The 
linear dashed (red) line is the well-known g2 — result N\a\/bt = ko. The solid 
(black) lines are for a narrow ■^^K resonance and trap length bt ~ 1.84/im, (bottom), 
0.5bt (middle), and 0.2bt (top). The cross and open square are as in fig. 7.2 
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To exemplify with the ^^K resonance we use B — Bq = — 32.3mG (cross 
in fig. 7.3) giving a = — 0.334ao and 52/^? = —0.0461. The critical particle 
number is then = 3.58 x lO'*. The values are indicated by crosses in 
fig. 7.2 and 7.4 The (72 = prediction is Nc = 5.975 x 10'*. In conclusion, 
we predict a 40% reduction in Nc for macroscopic particle number Nc ~ 10*. 



7.7 Macroscopic Quantum Tunneling 

Close to a = 00, g2/b^ is still a small term in typical traps. As the MQT 
rate |MMT07] is exponentially dependent on the integral under the barrier, 
we expect changes due to non-zero g2 to be amplified. A similar point was 
raised in |UL98j where the MQT rate was shown to increase dramatically 
close to Nc- The MQT rate, can be obtained from field theory ^Sto97j by 
considering the bounce solution of the effective action in variational g-space 
with potential E{q) /Nhuj. For comparison, we denote by Fq the rate with 
g2 = 0. All MQT calculations are done within the variational approach. 

We are interested in metastable states with small barriers and large MQT 
rates, thus we have to work close to Nc |UL98j . Below we find that g2/bf is 
only important for the rate when is of order 10^ or less. Since the three- 
body loss depends on the density to the third power |Adh02) . we expect it 
to be small in the outer metastable minimum and to be large in the inner 
one. Thus, the physical picture is that of a metastable BEC in the outer 
minimum that knocks on the barrier as a coherent state with a common 
quantum tunneling probability. When it tunnels to the inner minimum it 
swiftly decays as the large increase in density amplifies the three-body loss. 
Thermal fluctuations can of course also induce MQT and we therefore need 
to operate at very low temperature. An estimate of the thermal tunneling 
rate is given in |Sto97j . and leads to the following criterion for thermal 
fluctuations to be suppressed: E{qm) — -E'(q'o) ^ ksT, where qm and qo are 
the positions of the barrier maximum and the outer minimum, respectively. 



For the ^^K resonance we now consider B — Bq = 20mG (see fig. |7.3[ ) 
where 52/^? = —4.3 x 10~* and a = — 93.6ao. The critical particle number 



becomes Nc ~ 249 (see fig. 7.4). We now pick N = 242 < Nc to obtain 



metastability (see fig. 7.2). The rates are found to be F2 = 0.49 s ^ and 



Fq = 0.05 s^ (shown by the open squares in the inset of fig. 7.3), i.e. a 
tenfold enhancement of the MQT rate. The temperature must be below 
8 nK for thermal fiuctuations to be small. It is possible to obtain larger 

with the cost of a smaller effect on the MQT rate. However, going to 
N ~ 10^ does not seem possible. 

We propose to start from the a > side where the condensate is stable. 
Preparing a sample with low enough N and temperature of a few nK for the 
effect to be observable is the very difficult challenge for the MQT scenario. 
A sudden ramp of the magnetic field to the appropriate a < is then 
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performed. The system has to retain a small barrier so that the MQT rate 
is considerable and to avoid collapse. Then one monitors A'^ as a function 
of time, as was done in e.g. |DCC"'"01 . The low N is also challenging 



for optical imaging of the atomic cloud, but should be possible with current 
techniques. We expect not only to see modified MQT rates, but also changes 
in the decay after the condensate tunnels to a high density state: The three- 
body recombination is very sensitive to the density, and the inner barrier 
caused by the 52 term should affect this process as well. 



7.8 Other Signatures 

From the structure of the 52 term it is clear that density variations are needed 
to see effects. As we have demonstrated already, tighter traps amplify the 
contribution, and we expect optical lattices to do the same. Solitons and 
vortices are other features with density variation. For a simple one vortex 
state, the |^'pV^|^'p term will leave the core and the asymptotic regions 
unchanged but change the profile in between. Work is in progress |ZT09a] . 



7.9 Conclusion and Outlook 

We have explored the effect of higher-order terms in the Gross-Pitaevskii 
description of a Bose-Einstein condensate, particularly the effective range 
correction near a Feshbach resonance. Using both a variational and numer- 
ical approach we find an interesting new phase diagram for the stability of 
a condensate with negative scattering length. 

The critical particle number is strongly affected near zero-crossings of the 
scattering length for narrow resonances. Effects of 40% are found for particle 
number of order 10^. These deviations increase for more narrow resonances 
or tighter traps. We find that the critical number will be reduced as the 
zero-crossing is approached from either side since the higher-order term is 
attractive here. Macroscopic quantum tunneling is also modified by these 
higher-order corrections and we have discussed some experimental conditions 
for exploring the physics. Narrow Feshbach resonances are the best way to 
isolate the effect of tunneling. However, one needs small samples of order 
10^ particles or low temperatures. Other possible experimental signatures 
are tight traps, optical lattices, solitons, and vortices in rotating BECs. 

In a more general sense, the fact that higher-order interactions can be- 
come large and dominant in determining the stability properties means that 
it is no longer a correction and that still higher terms might become im- 
portant. This is the subject of future work. For the moment we have 
demonstrated that experiments targeting the regions discussed above could 
find interesting new stability properties beyond those that were already un- 
derstood about a decade ago. 
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Chapter 8 

Higher-Order Thomas-Fermi 
Approximation 



8.1 Introduction 



The Gross-Pitaevskii (GP) equation |DGPS99l lP502l IP503] has been ex- 
tremely successful in describing a wide range of mean-field features for exper- 
iments with Bose- Einstein condensates (BECs). In particular, the Thomas- 
Fermi (TF) approximation ^BP96| IDPS961 ILPS97| . where the kinetic energy 
is neglected, has been very rewarding |VHBL"'"98 . This approximation holds 
for repulsive condensates with positive scattering length a and large parti- 
cle numbers. In the regime of validity of the TF approximation, the total 
energy is distributed between interaction energy and potential energy from 
the confining trap, while the kinetic energy becomes negligible. 

Because of the non-linear nature of the GP equation, it is only solved 
analytically in a few cases, e.g. vortices and solitons in homogeneous con- 
densates |PS021 IPS03j . The TF solution is also analytical, although it only 
holds in the bulk of the condensate. At the surface the approximation breaks 
down and is usually patched by including the kinetic energy at the surface 
|DPS96|lLPS97l . 

The interactions of the ordinary GP equation is based on the lowest order 
zero-range potential, which is governed by the scattering length alone. Al- 
though this approximation is usually very good, the higher-order corrections 
to the scattering dynamics |RFOH |FWG03| lGMP07j can be crucial in certain 
cases, e.g. for Rydberg molecules embedded in BECs fCMP07] and for nar- 
row Feshbach resonances |ZT09b| . Inclusion of higher-order terms are well 
known and applied in Skyrme-Hartree-Fock calculations in nuclear physics 
|BGH85j . Here they often play a crucial role in order to get bulk nuclear 
properties right |Sky56[ ISJ87] . However, the effects of similar higher-order 
terms in the GP equation have been less investigated. 

In this chapter we solve the modified GP equation with higher-order in- 



95 



8.2. MODIFIED GROSS-PITAEVSKII EQUATION 



teractions analytically in the TF approximation. The chapter is organized 
as follows: We introduce the modified GP equation and its parameters and 
show how it is derived from an appropriate energy density functional with 
careful treatment of boundary terms. We present the analytical solution 
in the TF approximation and discuss the condensate size and chemical po- 
tential as function of the interaction parameters. The density profiles and 
energies are discussed and we address the consistency of the TF approxima- 
tion by considering the kinetic energy of the solutions. We compare to some 
relevant atomic systems and finally present our conclusions. 



8.2 Modified Gross-Pitaevskii Equation 



We assume that the condensate can be described by the GP equation. Since 
we are interested in the ultra-cold regime, where the temperature is much 
smaller than the critical temperature for condensation, we adopt the T = 
formalism. In order to include higher-order effects in the two-body scattering 
dynamics we use the modified GP equation derived in |iCMP07j . which in 
the stationary form reads 



U) 



where m is the atomic mass, Vext is the external trap, Uq = Anh'^a/'m, and 
52 = — ai?e/2, with a and Re being respectively the s-wave scattering 
length and effective range |CMP07] . We assume an isotropic trap, Vextir) = 
mw^r^/2, and introduce the trap length bt = ^/h/muj. The single-particle 
density, p(r) = |*I'(r)|^, is normalized to the particle number, N = J drp(r), 
and /i is the chemical potential. 

As the boundary conditions are important for the TF approximation 
applied below we now discuss the procedure for obtaining the modified GP 
equation from the corresponding energy functional which is 



E{^) = j dr(ei^ + ey + e/ + e/2), 
with kinetic, potential and interaction energy densities 



(8.2) 



^K = i^\y^\, ev = Vext{r)\^\ (8.3) 

e/ = ^^ol^l', e/2 = ^t/o52|^|'v2|^|2. (8.4) 

The corresponding integrated energy contributions are denoted Ek, Ey, Ef, 



and Ej2 respectively. To obtain eq. (8.1), we vary eq. (8.2) with respect to 
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^* for fixed ^. To first order in 5'^* we have 



5E = E[^* + 5^*] - E[^*] 



dr 



2171 



Here S is the outward-pointing surface normal. In the usual analysis one 
assumes that ^ and V\I' vanishes at infinity, drops the boundary terms, and 
eq. (8.1) is obtained by varying E — fiN. However, the existence of these 



surface terms are essential for the inclusion of higher-order interactions as 
discussed below. 

In the rest of this chapter we use trap units, huj = bt = 1, i.e. energies 
{E, Vext, etc.) are measured in units of Huj and lengths (a. Re, r, etc.) in 
units of bf. Note that 52 has dimension of length squared. 



8.3 Thomas-Fermi Approximation 

Let us briefly review the standard Thomas-Fermi approximation jBP96l 
IDGPS991 IPS02| IPS03] . Neglecting the kinetic energy term, as compared 
to the trap and interaction energies, the GP equation has the solution 

PTF = ^i^^TF - ^r^), (8.6) 

with chemical potential ^tf- This solution is used out to the surface, 
Rtf, while outside ptf = 0. The normalization and surface condition 
Ptf{Rtf) = gives 



PTF 2- 



Ir^p, Rtf = (15A^a)^/^ (8.7) 
The total energy becomes 

^TF 5 R'^p 

The trap and interaction energies are Ey = 3E/5 and Ej = 2E/5, respec- 



tively. Since Rtf > in eq. (8.7), these results only hold for a > 0. The 
TF approximation is good for Na ^ 1, except at the surface region where 
the kinetic energy density diverges. Here the solution can be corrected as in 
|DPS96l ILPS971 lP502ll^03] . essentially giving a smah exponential tail. 
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8.3.1 Inclusion of Higher-Order Interactions 



We now consider the TF approximation with the higher-order interaction 
term, e/2. Ignoring the boundary terms in eq. (|8.5p, the modified GP equa- 

as 



tion can then be written in terms of the density p{r 

1 



|^(r)i2 



A* 



+ Aira {p + 52 V^) • 



(8.9) 



With scaled coordinate x = r j ^Jlfi (assuming 52 > for the moment) and 
density /(r) = Airaxpir) / §2-, this becomes 



dx^ 



+ / 



92 



1 



.10) 



The inhomogeneous and homogeneous solutions with boundary condition 
/(O) = are 



m = ( 



^92 



1 



+ 3)x, fh{x) 



A 

92 



■ sm X, 



Ul) 



where A is a constant (with dimensions of length squared) to be determined 
later. The full solution is 



p{x) 



92 
Ana 



92 



1 



-x^ + 3 + 



A sin X 

92 X 



.12) 



For a given A, the chemical potential p and the condensate radius R are 
determined by the normalization and the surface condition, 



f-xo 

/ 4:'Kx'^p{x) dr = N and p{xo) = 0, 
Jo 



i.l3) 



where xq = Rj ^J~g2- The solution p should be positive for x < xq which 
must be explicitly checked. Outside xq we use p = 0. 

We now consider the boundary terms in eq. (8.5). Above we assumed 
that /o(xo) = at some finite radius xg which we identify as the condensate 
size. However, only the first two boundary terms in eq. (8.5) vanish on 



account of this condition. For the last term in eq. (8.5) to vanish we need 
Vx^'(xo) = 0, which implies that 



dp 
dx 



(xo) = 0. 



^.14) 



Notice that this latter derivative is in fact non-zero in the 52 = case, 
which is the root of the divergence of the kinetic energy at the condensate 
surface as we discuss later. Equation (8.14) gives a closed expression for the 
remaining free parameter A, 



A 

92 



X 







Xq cos Xq — Sm Xq 



.15) 
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This additional requirement on the derivative at the edge of the conden- 
sate implies that higher-order terms require a smoothing at the surface of 
the cloud. In addition, the discussion of which kinetic operator structure to 
use (|V^p or \I'*V^^' |LPS97j ) is obsolete in our treatment since the bound- 
ary term 6^*V^ vanishes. In this sense the inclusion of a higher-order term 
neatly removes some of the difficulties of the traditional TF treatment. 

The solutions with a finite boundary R of the modified GP equation only 
minimize the energy functional if eq. (8.14) holds. We note that extremal 
states of the energy functional always satisfy the virial theorem. Thus, 
enforcing the virial theorem on the GP solutions is equivalent to eq. (8.14). 
For completeness, we show in appendix [C] that the virial theorem approach 
also leads to eq. 



^.15). 



8.4 Size and Chemical Potential 

We now determine the condensate size R and chemical potential /i. The 
normalization condition is 



Na 

5/2 
92 



^3 

Xq 



3^2 



10 



while the surface condition reads 



ii 

92 



xl/2 + 3 + 



A sin Xq 

92 Xq 



0. 



Combining eq. (|8.15|), eq. (|8.16|), and eq. (|8.17|) gives 

1 + 



Na 



5/2 
92 



xl/3 



15 



1 — Xq cot Xq 



i.l6) 



.17) 



.18) 



which determines R for given Na and §2, and upon back-substitution also 
fi. 

The 52 < case can be worked out analogously by replacing trigonomet- 
ric functions with hyperbolics and keeping track of signs. The two cases can 
in fact be combined into one equation 



Na 



\92 



5/2 



|Xo|Xo 



15 



1 + 



xl/3 



1 — \xo cot Xo\ 



(8.19) 



This equation determines Xq 
The result is shown in fig. 



= R? /g2 implicitly as function of Na/\g2\^^'^ ■ 
1 ) . We notice that in principle R becomes a 



multi-valued function. However, all the higher solutions for 52 > (dotted 



in fig. (8.1 )) are spurious, since the density becomes negative on one or more 
intervals inside R. The non-spurious solutions (solid line in fig. (8.1 )) defines 



i? as a single- valued function of a and g2, which was not guaranteed a priori. 
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The four quadrants in fig. (8.1) correspond to the different sign combinations 



of a and 52- The sign of the extra interaction energy, Ej2, is determined by 
052 V^p. For a typical concave density the Laplacian term will be negative. 
We therefore see that for ag2 > the higher-order interaction is attractive, 
whereas for ag2 < it is repulsive. The TF solution only exists for ag2 < 0. 
We discuss both cases separately below. 




-10000 -5000 5000 

Na/\g2\'/^ 



10000 



Figure 8.1: Condensate size (i?) as function of Na and 52 as found in the modified 
TF approximation, eq. ( |8.19 1. The solutions (a) and (b) correspond to the sign 
combinations (a > 0, (72 < 0) and (a < 0, (72 > 0), respectively. No solutions exist 
for ag2 > 0. The spurious solutions (dotted) have negative densities for one or more 
intervals inside R. The branch (a) approaches the normal TF result eq. (8.7 1 when 
Na -foo or 32 — > —0. Note that the convergence is only relative, see eq. (8.191, 



and the TF limit is better represented in the logarithmic inset. Points indicate the 
data from tab. 



ll. All values are in trap units. 



8.4.1 The Attractive Regime: ag2 > 

For a < 0, 52 < (third quadrant in fig. ( |8.1| )) there are no solutions, which 
is expected since the normal TF approximation has no solutions for o < 
as the interaction energy Ej is negative and the kinetic energy that could 
prevent collapse is neglected. 

The 52 > 0, a > case in the first quadrant has only spurious solutions. 
Here the 52 term is attractive for the typical concave density and a collapse 
towards a high-density state is possible in complete analogy to the usual 
discussion of attractively interacting condensates within the standard CP 
theory. Whereas there can be metastable states at large values of Naj ^ 
these are stabilized by kinetic energy and thus are not present in our TF 
approach. Thus, even when the the total kinetic energy is small, it is still 
needed to prevent the attractive higher-order term from amplifying local 
density variations. 



100 



CHAPTER 8. HIGHER-ORDER THOMAS-FERMI APPROXIMATION 



This important point can also be established by considering the stability 
of the homogeneous condensate through linearization of the GP equation. 
By repeating the analysis of [PS02^ with the higher-order term, we find that 
for 52 > and a > the kinetic energy term is crucial for the stability of 
the excitation modes. In fact, exponentially growing modes will always be 
present if the kinetic energy is neglected. This will be discussed elsewhere 
in relation to the numerical solution of the full GP equation [ZT09aj . 



8.4.2 The Repulsive Regime: ag2 < 



For 52 < 0, a > a single solution (a) exists. This was expected since 
Ei2 > gives extra stability. The solution approaches the normal TF result 
^.7) when Na/\g2\^^'^ — > +oo, as can also be seen from eq. (8.19). 



m eq. 



Of course in this limit Ej2 *C Ej. However, the co nverg ence in terms of 
Na/\g2\^/^ is only on a relative scale, see inset in fig. 



1) and eq. (8.19). 



For 52 > 0, a < there is a single solution (b) which connects smoothly 
to the (a) solution. In the limit Na/\g2\^^'^ — oo, which is determined by 
xocotxo = 1, we find B? / g2 = 20.1907. This solution is possible when the 
52 term provides just enough repulsion to cancel the usual a < collapse 
behavior. 



8.4.3 Chemical Potential 



In fig. (8.2) we show the chemical potential for the smoothly connecting 
solutions (a) and (b). Again we see that (a) approaches the normal TF 
limit for large A^a/|52|^^^. Here it is interesting to note how /x turns around 
near the origin (amplified in the inset in Fig ( |8.2[ )) and maintains a positive 
value. This occurs in the region where the lowest-order interaction gives 
a large negative energy contribution which the 52 term is still able to bal- 
ance yielding a well-defined TF solution. This behavior is analogues to the 
balancing of attraction by the kinetic term in the usual a < 0, 52 = case 
|BP96| IDGPS99] . As a becomes increasingly negative so too does jj, and 
collapse is inevitable (and likewise when 52 — > O"*"). 



8.5 Densities and Energies 



With R and determined we can find the density profile, energy densities 



and integrated energy contributions. With eq. (8.12) the energy densities 
are given by 



e/ = 27rap , en 



1 , A sin X , 
■-p3 + . 

2 52 X 



.20) 
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CO 



g2 > 



g2 < 




-10000 -5000 5000 

iVa/|g2|5/2 



10000 



Figure 8. 2: Chemical potential /i as function of Na and 52 as found in the modified 
TF approximation, using the solutions (a) and (b) from fig. (8.1 1. For branch (a) 
and the upper part of branch (b) (see inset) we have /i > 0. 
has /i < 0. Points indicate the data from tab. (8.1 1 



The lower part of (b) 
All values are in trap units. 



Using eq. (8.9) the total energy density (without ex) becomes 



e = ev + ei + en 



\p{x){V,,t{x) + ^). 
2 92 



.21) 



In fig. 



^.3 ) we show the density profile of the (a) solutions for Na = 10^ 
We clearly see that the higher-order term tends to 



and selected 52 < 0. 
expand the condensate through its repulsion. Importantly, at the boundary 
there is a smoothing caused by the condition in eq. (8.14), see inset in 
fig. (8.3). We will discuss how this affects the estimated kinetic energy in 



the next section. As \g2\ grows we see the condensate flatten and in the 
limit of very large \g2\ it becomes a constant density. 

Figure (8.4) displays the density profile for the (b) solutions with a < 
for selected 52 > 0. Here we see the profile collapse towards the expected 
delta- function with decreasing §2- It is interesting to follow the (a) solution 
through the origin in fig. (8.1 ) and onto branch (b), passing from §2 = —00 
to 92 = CO- On the (a) branch the solution flattens as 52 decreases and 
eventually becomes effectively constant in space. This is also true for the 
(b) branch at (^2 = 00, and as §2 is decreased the solution proceed to shrink 
as the (72 term becomes unable to provide the repulsion needed to prevent 
the a < collapse induced by the lowest-order term. 

From the figures we see that large \g2\ induces large changes in cloud size. 



As the condensate can be imaged with very good resolution |VHBL"'"98| , this 
should be measurable if the regime of large \g2\ can be accessed. 

We now discuss the energy contributions which are interesting since the 
release energy are in fact measurable quantities |PS03j . Since we neglect the 
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Figure 8.4: Same as fig. (8.1 1 but for solutions (b), i.e. opposite signs 92 < and 
Na = -lO''. 
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kinetic term in the TF approximation, the release energy is simply Eji = 
E1+E12 = E — Ey. In tab. (8.1 ) we give the integrated energy contributions 
for some relevant values of §2 calculated for A^|a| = 10^, whereas fig. (8.5) 
gives the energies as function of Na/\g2\^^'^ ■ We note that for smaller values 
of A'^lal the same overall behavior is found, however, the kinetic term is more 
important and the TF approximation becomes worse. 



52 > 



52 < 



40 
30 
20 
10 


-10 
-20 
-30 

-300 -200 -100 100 200 300 400 500 

Na/\g2\^/^ 

Figure 8.5: Different total energy contributions. N\a\ — lO'*. Values are in trap 
units. 
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We observe that E/N grows towards the I52I =00 point. This is due to 
the trap energy increasing as the density flattens {Ey diverges around the 
origin in fig. (8.5)). Furthermore, as §2 — > 0"*" the energy diverges towards 
—00 as the collapse sets in {Ej diverges on the 52 > side in fig. (8.5)). The 
boundary where the energy vanishes is around 52 = 5.14 for A^|a| = 10^, but 
this depends on the choice of A^|a|. With respect to the release energy, we 
find that somewhere in the region 10 < g(2 < 50, En becomes negative. This 
is a result of the unavoidable collapse, and also indicates that kinetic energy 
cannot be ignored at this point. Notice, however, that the release energy 
changes considerably and could provide a way to measure the influence of 
the 52 term. 



8.6 Thomas-Fermi Approximation: Consistency 

We now address the validity of the TF approximation with the §2 term 
included. In order to do so we must consider the contribution of the kinetic 
energy. The kinetic energy density can be written 

02 ( Ax sin X — cos x \ ^ , , 

= « U ^2 0^ + • 8.22 

8/>(47ra)^ V 52 x J 
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Strictly speaking, this is not the true kinetic energy, since the kinetic term 
were neglected from the start. However eq. (8.21 ) and eq. (8.22 ) can be used 
to test whether the TF approximation holds locally, i.e. e;^ <C e should hold 
for the solution p to be consistent. In tab. (8.1) we calculate the integrated 
contribution of the kinetic energy relative to the total TF energy, and we 
find that the contribution is small everywhere except the point where E = Q 
on the g2> ^ side of fig. (8.5 ). Here the kinetic energy is of course the most 



important term and the TF approximation is poor. 

In the standard TF, the kinetic energy causes trouble at the boundary 
of the cloud. Here V^' oc Vp/y^ and since the density vanishes and the 
derivative is finite (see eq. ( |8.6[ )) this diverges at Rtf- When including 
the higher-order term we need to use the additional boundary condition 
V^' = at i?, so the kinetic energy will be strictly zero at R. However, 
as one approaches the boundary the kinetic energy density grows rapidly 
before it descends towards zero within a very small interval at R. The total 
energy density in eq. (8.21 ) goes to zero at this point and we find that e^^/e 
is very large near the boundary as in the usual (72 = case. 

We conclude that the inclusion of the higher-order term does not alleviate 
the difficulties with kinetic energy at the boundary. The techniques for 
addressing this problem described in |DPS96| ILPS97j should therefore be 
generalized to include the higher-order interaction term in order to improve 
the description at the boundary of the cloud. 



8.7 Comparison to Atomic Systems 

The considerations above show that deviations from the usual TF approx- 
imation can be strong when 52 is large. In the following we reintroduce 
explicit units for comparison with real systems. We have to consider g2/bf. 
Of course the factor means that this quantity is generally very small since 
52 is of order Oq and bt is of order lO^ao. 

We first consider some typical background values for bosonic alkali atoms 
away from Feshbach resonances. We estimate the effective range to be of 
order the potential range, and assuming a van der Waals interaction we 
have i?e ~ 50 — 200oo. For typical one-component gases we have — 450ao < 
a < 2500ao |CGJT09| . Since c/2 = - aRe/2, we see that the term 
will dominate and in all cases < §2 ^ lO^ao. In trap units this becomes 
92/bt ^ 5 X 10~^(l/xm/6t)^. In typical traps of 6f ~ 1 — 10/im the higher- 
order term is therefore very small. These values also predominantly lie in 



the first quadrant of fig. (8.1 ) and thus no TF solution exists. 

We now consider Feshbach resonances in order to increase the influence 
of the g2 term. We use a multi-channel Feshbach model |BJK05| IZT09bj , see 
chapter [7] which describes both a and Re as a function of resonance position 
Bq, width AB, magnetic moment difference between the channels A/i, and 
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the background scattering length abg. We have a = abg{l — AB/{B — Bq)) 
and Re = i?eo/(l - {B - Bo)/AB)'^, where Reo = -2h'^ /mabgAfiAB < 0. 
Notice that Re = Reo{^ — cibg/a)'^ and 

92{a) = -- —(1 - — ) . (8.23) 

Hence g2 diverges when a ^ (referred to as zero-crossing) or a — > oo 
(on resonance). Near a zero-crossing g2 behaves as ag2 — \Reo\alg/2, where 

ReO < 0. 

As a concrete example, we consider the alkali isotope ^^K where several 
Feshbach resonances of vastly different widths were found recently jDZF"'"07 



First we focus on zero-crossing and consider the very narrow resonance at 
Bo = 28.85G with AB = -0.47G, A// = 1.5hb and abg = -33ao. We obtain 
Reo = — 5687ao and ag2 93.8 x IO^Oq for a ^ 0. It is important to notice 
that ag2 > around a = 0. This means that we are looking for solutions 



in the first and third quadrant of fig. (8.1), and again we have to conclude 
that no TF solutions can be found when higher-order terms are taken into 
account. 

Another case of interest is around resonance where \a\ = oo. Here we 
have Re ^ ReO and g'2 oc > on both sides of the resonance. Thus the 
a > side will be in the first and the a < in the second quadrant of 



fig. (8.1). This makes it difficult to imagine sweeping the resonance from 



either side to probe the solutions on branch (b) in fig. (8.1). One could 
imagine starting on the a > side with small 52 > 0. The full GP equation 
will have perfectly sensible solution here, however, when one approaches 
the resonance the 52 term will diverge and induce collapse already on the 
a > side. If we approach from the a < side then we face the problem 
that the critical number of particles decreases dramatically before 52 grows 
sufficiently, and one therefore needs a very small condensate since Na/bt ~ 
0.5 jZT09b] . At this point the TF approximation is no longer valid. 

From the examples above we see considerable problems in accessing the 
TF solutions presented above in current experiments with ultra-cold alkali 
gases. In particular, we notice that realistic systems which have been used 
for creation of BECs in alkali gases for the last decades have parameters that 



predominantly lie in the first quadrant of fig. (8.1). As we have discussed 
there are no well-defined TF solutions in that region. Therefore we see 
that the kinetic energy plays a decisive role and we are forced to consider 
it in principle, even if it is small for all practical purposes. The physical 
reason is that for a > and 52 > the higher-order interaction is effectively 
attractive and induces collapse which will have to be balanced by a barrier 
from the kinetic term, similar to the a < 0, 52 = case |DGPS99] . Since we 
neglect the kinetic term in the TF approximation we should not expect to 
find solutions in the 0*72 > case. 
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8.8 Conclusion and Outlook 

We have considered the effect of higher-order interactions in Bose-Einstein 
condensates within the Gross-Pitaevskii theory. We derived the GP equa- 
tion with effective range corrections included and solved it analytically in 
the Thomas-Fermi approximation. Higher-order interaction terms act as 
derivatives on the condensate wave function which means that the bound- 
ary conditions on the solutions of the GP equation must be carefully con- 
sidered. We then discussed the solutions for various parameters, presented 
the chemical potential, density profiles, and the energy contributions. 

We find that no TF solutions are possible when the higher-order term 
is attractive. This conclusion holds both in the trapped system and in the 
homogeneous case. An estimate of the relevant parameters for alkali atoms 
showed, however, that they typically lie in the region where the effective 
range correction is effectively attractive. In order for this to be consistent 
with current experiments that have found the standard TF approximation 
to be quite good we conclude that the kinetic energy, even if very small, is 
crucial in order to stabilize collapse due to higher-order interaction terms. 
A full numerical analysis is in progress |ZT09a| . 
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Chapter 9 

Summary and Outlook 



In this dissertation we presented theoretical investigations of universality 
and finite-range corrections in few- and many-boson systems. The main 
focus was on ultra-cold atomic gases seen both from the three-body, many- 
body, and mean-field perspective. 

In the introductory chapter [T] we gave a soft introduction to the concept 
of universality in few- and many-body systems, in particular the wide range 
of effects in atomic gases and Bose-Einstein condensates. We also focused on 
the search for the elusive Efimov effect and its relation to nuclei and atomic 
gases. In chapter |2] we presented the relevant theoretical and numerical 
background for universal few- and many-body systems. 

The purpose of chapter [3] was to go beyond the scattering length ap- 
proximation for three-boson Efimov physics and express the corrections in 
terms of model-independent parameters. We found universal scaling with 
corrections determined by the effective range. We showed that for nega- 
tive scattering lengths the effective range corrections to Efimov physics and 
Borromean binding are two aspects of the same quantitative effect. This 
leads to a linear shift in critical scattering lengths at the trimer threshold. 
For positive scattering lengths near the atom-dimer threshold the effective 
range corrections to the trimer energies are mainly determined by the dimer 
corrections. Results agree quantitatively with the newest effective field the- 
ory predictions. This has consequences for three-body recombination rates 
in atomic gases when the scattering length is comparable to the effective 
range. The main effect is that the universal scaling factor becomes smaller 
for the lowest Efimov states (smallest scattering length) on the negative 
scattering length side, corresponding to a smaller Borromean window. We 
also described the effects of putting the system in a finite trap. 

Chapter |4] continued this line by analytically investigating the conditions 
for Efimov physics, in particular for large effective range. We used the adi- 
abatic hyper-spherical approximation to derive rigorously a transcendental 
equation for the adiabatic potentials for general finite-range potentials. So- 
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lutions agreed with exact numerical results. We concluded that in the hyper- 
spherical adiabatic approximation it is insufficient to include effective range 
only, as another term of the same order appears. However, non- adiabatic 
corrections restores model-independence. For large negative effective range 
the window for Efimov physics is precisely open between the effective range 
(not the potential range) and the scattering length. 

In chapter |5] we provided theoretical support for the existence of long 
lived meta-stable A^-body Efimov states in atomic Bose-gases. The inclusion 
of two-body correlations (or equivalently a Faddeev-Yakubovski decompo- 
sition in two-body amplitudes) gave an effective hyper-radial potential as 
for the three-body Efimov effect. This in turn led to universal scaling, but 
with new scaling factors. Results were confirmed both numerically and ana- 
lytically. The experimental signatures in atomic recombination experiments 
were discussed including possible connections to other universal four-body 
effects. The four-body Efimov states would be harder to detect due to the 
low density of experiments, but also easier to probe because of the smaller 
universal scaling factor. The effect could be applicable to Borromean sys- 
tems too. 

In chapter [6] we considered trapped few-boson systems with large positive 
scattering length. This was modeled using attractive two-body potentials, 
giving the many-body system a large amount of low-lying bound states. 
We presented a novel numerical technique to identify the BEC state as a 
highly excited many-body state. The obtained state contained short-range 
correlations determined by the interaction, while the normal mean- field fea- 
tures were retained at long distances. We found that for scattering lengths 
smaller than the trap the system shows universality. For larger scattering 
lengths the system properties becomes independent of the scattering length, 
contrary to the zero-range and the repulsive models. The correlations be- 
come long- ranged and universal in the large scattering length limit. This is 
an important step towards understanding the non-universal corrections in 
condensates. 

Chapter [7] approached the question of effective range corrections in con- 
densates from the mean-field point of view. A modified GP equation with 
effective range dependence was introduced. Using variational and numerical 
approaches we found a phase diagram describing the condensate stability. 
We then considered an extended Feshbach resonance model including effec- 
tive range variations. Effects on macroscopic quantum tunneling were small 
for realistic systems. The behavior of the critical particle number was mod- 
ified near a scattering length zero-crossing with observable consequences. 

In chapter [8] we continued the mean-field analysis of the condensates 
using the modified Gross-Pitaevskii equation. We solved this equation an- 
alytically in Thomas-Fermi approximation where the kinetic energy can be 
neglected, but keeping higher-order interaction terms. Boundary conditions 
of the solutions were carefully considered. We presented the chemical po- 
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tential, density profiles, and the energy contributions. No Thomas-Fermi 
solutions were possible when the higher-order term was attractive. This 
holds both in the trapped and homogeneous system. However, realistic 
atomic systems typically lie in the region where the term is attractive. In 
order for this to be consistent with current experiments, the kinetic energy, 
even if very small, is crucial in order to stabilize collapse due to higher-order 
interactions. 

In conclusion, we have carried out new theoretical investigations of uni- 
versality and its limits in few- and many-boson systems. The focus was 
on ultra-cold trapped atomic gases, but results were presented in universal 
model-independent terms. Thus, much of the work may hopefully be used 
or continued in other areas of physics. The subject of universality in atomic 
gases will face many new challenges during the next decades. Interesting new 
physics is clearly within reach in ultra-cold gases and indicate a promising 
future. 



Note added. The observation of an Efimov spectrum in an ultra-cold gas 
has recently been pubhshed |ZDD+n9| giving the first definite proof of the 
universal scaling factors. We note that this experimental group is indepen- 
dent of the group measuring the first isolated Efimov state [KMW"'"06 



In |ZDD^09 five experimental features were found. For positive scat- 
tering length two interference minima were observed in the recombination 
loss rate, corresponding to two Efimov states. Also on this side two small 
peaks were found, corresponding to a secondary loss mechanism. For neg- 
ative scattering length a single peak was found in the loss rate (the second 
Efimov state being out of the experimental range). 

All five features agree well with the theoretical universal scaling factors 
except for systematic shifts. They conclude that the scaling factors are 
reduced on the negative scattering length side and increased on the positive 
side. They attribute this to finite range corrections and compare with the 
effects in |TFJ08c] (i.e. chapter |3]) and the identical conclusions of effective 
field theory |PJP09j . The theory and experiment agrees. A more precise 
analysis is needed at this point, but the need for effective range corrections 
to Efimov physics is more evident now than ever before. 



Ill 
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Appendix A 

Calculation of Condensate 
Fractions 



To calculate the condensate fraction of the states found with the stochastic 
variational method in chapter [6] we first need to calculate the one-body den- 
sity matrix (OBDM), eq. (|6.1[), for the two-body correlated wave-functions, 



eq. (2.77). We choose to change basis to the harmonic oscillator eigen- 
functions before diagonalization. This basis turns out to be good, since the 
OBDM is almost diagonal and numerically only a few basis states are needed 
to ensure that the occupation numbers sum up to A^. 

Here we give an outline on how the matrix elements are calculated and 
the condensate fraction is obtained. First we expand the eigenfunctions of 



eq. (6.2) in the harmonic oscillator basis, 

X^{r) = Y,4LRni{r)YU0,^), (A.l) 



I'lri 



Rnl{r) = ,1 ^ „ e'i^''-'(gr)'Li+^(gV), (A.2) 




where Yim are the spherical harmonics, L„ ^ are the associated (or gener- 
alized) Laguerre polynomials [AS95j . and g is an arbitrary oscillator scale 
(typically chosen as b^^). In this basis the eigenvalue equation eq. (6.2) 
takes the form 

n'/'m' 

where the matrix elements are given by 

Anim,n'Vm' = J n{r , v') R^^n' {r') Ki{r)Yi,^, {n')YC^{n) dr dr' , (A.4) 

with dr = dr dil. Most work lies in the evaluation of these matrix ele- 
ments. It can be shown that the OBDM (in the r-basis) for the two-body 
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correlated wave- function eq. (2.77) has the form 



9 

n(r,r') = X] 



,kk' 



exp 



(A.5) 



kk' t=l 



kk 



, dfl'^' } are functions of N and {Cfc , a^'^^ /J^^) } 
k' corresponds to the linear-combination in 
eq. (|2.77l). The 9 different terms labeled by t comes from the symmetriza- 



The parameters {a^^' , b^^' , c\ 
The summation over k and 



tion, which was done analytically. By expanding the factor exp((ij*^V • r') in 



eq. (A.5) in spherical harmonics and inserting it into eq. (A.4), it is straight- 



forward, although tedious, to calculate the matrix elements. The result is 



^nlm,n'l'm' —OU'Omm' I ) I / ) / ^ 

\Ti/ \n J 



1 1 

* 3/2+Z+n+n' 



a. 



kk't y++ (A.6) 

X 2Fi{—n, —n; —k; 1 



where 



y++ ={bf +g^/2){cf +9^2) - {df /2)\ 
={bf + 9V'^){cf - - {dfl2)\ 

y-+ ={bf - a'micf + g'm - {df /2)\ (^-^^ 

x=g^d\^'/2, k = n + n' + l + ^. 

The hypergeometric function 2-^1 is here reducing to a Jacobi polynomium 
in the last argument of degree min(n, n') |AS95j . 

We see that A is diagonal in / and m, and also independent of m. This 
is because eq. (2.77) is rotationally invariant. We can therefore diagonalize 
one block at a time. The eigenvalues and eigenfunctions are re-labelled as 
g^J^^-^ T^^'™^(r), with j labeling the eigenvalues in decreasing order. The 

condensate fraction will usually be the first eigenvalue Aq = Nq'^^^ /N. 

The matrix A can now be diagonalized numerically and the condensate 
fraction obtained as the largest eigenvalue. 



^The exact forms of the functions are analytical but very complicated. It is more 
convenient to do the transformation partly numerically. 
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Appendix B 

Examples of Strongly 
Correlated Wave- Functions 



The correlated wave-functions used in the stochastic variational method 
opens up the possibility for the condensate to be depleted considerably. It 
is therefore useful to know how the different types of correlations are able 
to lower the condensate fraction. In this appendix we give some examples 
of strong hyper-radial and two-body correlations. 



B.l Strong Hyper-Radial Correlation 



If only one term {k = 1) is included in the hyper-radial variational function, 



eq. (2.75), it reads 



ci exp 



2h\ 



exp 



(B.l) 



In this case the occupation numbers Ni defined in eq. (6.2) can be obtained 



analytically. Since we are considering a spherically symmetric trap, the total 
angular momentum / and its projection m are good quantum numbers, and 
we can re-label the occupation numbers as _/vj''"\ The subscript j labels the 
eigenvalues in decreasing order for given /, m. By expanding in the harmonic 
oscillator basis as in appendix [A| the result is 



ilm) 



8NC 



l+2j 



[1 + VTTI]2'+4j+3' 



with 



7 = Na^^h^. 



(B.2) 



(B.3) 



7 iV2 ' 

Here a^^-* is the non-linear parameter in eq. (B.l). Note that 7 = 1 corre- 
sponds to the fully condensed gas where A''- = 5jo5iQ5moN . When 7 is 
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far away from unity the condensate fraction can be made arbitrarily small. 
This means that the condensate can be destroyed just because the internal 
length scale (measured e.g. by y/ < >) is different from the center-of-mass 
length scale. This conclusion still holds for more /c-terms in eq. (2.75). The 
same conclusion is also found in |Gaj06| . 



B.2 Strong Two-Body Correlation 

To investigate the impact of two-body correlations on the condensate frac- 
tion we choose one term {k = 1) in eq. ^rri\ with iVa(^)6j2 ^ i_ ^j^g 
wave-function then has the form 

^2B = IJexp (^-^) X Sexp (-^/J^'V?^^ . (B.4) 

This is a special case of the more general form 

N 

^ = \{<l^{r,)xSf{r^2), (B.5) 



i=l 



which we will investigate instead. It is a mean-field product state modified 
by a correlation function. The real correlation function / only includes one 
pair ri2, and the function is subsequently symmetrized. We can assume 
without loss of generality that and / are normalized to unity, 

j |</.(r)|2dr = j |/(r)|2dr = l. (B.6) 

We now consider the strong two-body correlated limit, where / goes to zero 
over a range much smaller than the range over which cj) varies, say 

fir) ~ for Irl > e, 
II 

(j){r + 6r) ~ (j){r) for \6r\ < e. 
Here e is some small length going to zero. In this limit we are able to calculate 



the eigenvalues and eigenfunctions of the OBDM, eq. (6.1). Since eq. (B.5) 



is not normalized we must first calculate the overlap integral. Because of 



the symmetrization in eq. (B.5) we get \N{N — l)/2]^ terms. Many of these 



terms are identical, however, and only 3 different types occur, namely 
($1$) = Af(Af-l)[(12|12)-F2(iV-2)(12|13)-F(iV-2)(7V-3)(12|34)], (B.8) 
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where 

(12|12) = J |0(ri)|2|(/>(r2)|V(r-i2)'dridr2~ J |0(r)|^dr, 
(12|13) = I |</.(n)|2|</,(r2)|2|</.(r3)P/(n2)/(r-i3)dridr2dr3 

^ J \Hr)fdriJ /(r)dr)2, 
(12|34) = J |</.(ri)|2|</>(r2)|2|</.(r3)|2(/,(r4)P/(''i2)/(r-34)dridr2dr3dr4 

c^ij l^irtdrfij /(r)dr)2. 

(B.9) 

The terms J f{r) dr are of order e. To lowest order in e we find 

($1$) ~ Ar(Ar _ 1)(12|12). (B.IO) 

The symmetrization of the non-normalized OBDM (denoted n(r,r')) can 
be done in a similar manner. It consists of 9 different types of terms, but to 
lowest order in e only two terms remain, 

n(r, r') = (A^ - 1)(A - 2)n2323(r', r') + (A^ - I)ni2i2(r-, r'). (B.ll) 

These two terms are given by 

n2323(r-,r') = j 0*(r)(/>(r')|0(r2)|2|</.(r3)| V(r-23)' dr2 dr3 
= 0*(r)</.(r')(12|12), 

/(B 12) 
(t)*{r)4>{r')\<i){r2)\^f{r - r2)/(r' - r2) dr2 

~ j \4,{s)\^f{r-s)f{r' -s) ds. 

The normalized OBDM n{r,r') = An(r, r')/($|$) finally becomes 

n{r,r')^{N -2)<t>*{r)(t>{r) + j ds NsX*s{r)X8{r'), (B.13) 



where we have defined 



Xs{r) = f{r-s), Ns = 2 ^ . (B.14) 

J n(rj^ dr 



Here n(r) = |(/)(r)p is the mean-field density. Notice that eq. (B.13) has 
trace equal to N. The set {(l),xs\s S M^} is orthonormal in the limit con- 



sidered. Since the expansion of the OBDM in eigenfunctions, eq. (6.3), is 
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unique, we can compare eq. (6.3) with eq. (B.13) and read off eigenfunctions 
and corresponding eigenvalues directly. The condensate fraction is reduced 
to 1 — 2/N with (j) still being the corresponding single-particle state. The 
remaining 2 "particles" are distributed over infinitely many states Xs{f)i 
which are peaked around s. The occupation numbers dsNs for Xsif) are 
proportional to the square-density, n(s)^, of the mean-field. 
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Appendix C 



TF Boundary Conditions 
and the Virial Equation 



In this appendix we derive the free parameter A of the Thomas-Fermi so- 
lution eq. (8.12) in chapter ^ This is done explicitly by enforcing the 



virial equation. We show that this is equivalent to the boundary condition 
eq. 



.14) 



Even though eq. (8.12) is a solution to the modified GP equation eq. (8.9) 



for all A, it does not necessarily minimize the energy functional as discussed 
in the section 8.3 This can also be seen from the virial equation (with 



neglected kinetic energy), 



2Ev + 3Ei + 5Ei2 = 0, 



(C.l) 



which holds for all extremal points of the energy functional. Equation (C.l ) 



is derived from the energy functional using scaling arguments as in |PS03 

As an example consider the A = solution in eq. (8.12). This solution 
has a chemical potential shifted by 3g2 compared to the (72 = TF result. 
But the density is unchanged and so is Ey and Ej. Hence the usual virial 
equation —2Ev + 3Ej = for 5(2 = also holds for (72 7^ 0. Since Ej2 = 





8(72/2 7^ the virial equation eq. (C.l ) is not fulfilled, and hence the A 



solution is not extremal. Below we use the virial equation to calculate the 
value of A that minimizes the energy functional and the corresponding R 
and /i. We will also prove that this condition is in fact equivalent to assume 
p(xo) = Vxp{xo) = at the boundary. 

The general results for A, R, and 11 can be derived using the normaliza- 
tion and surface conditions eq. (8.13), and the virial equation eq. ( C.l| ). For 
convenience we introduce the variables fl = p/{2>g2) + 1, A = A/(2>g2), and 
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5 /2 

c = Nalg^ . The different energy contributions are 

Ey = -is r dxx\fl - — + A — ), 

Jo 6 X 

f'-^o ^ sill X 

Ei = 9s dxx'^ifL +A f, 

Jo 6 X 



El2 



^° , n, -i-sinx,,, ^sinx, 
dxx^ili - — + A )(1 + A ; 

D X X 



(C.2) 



wliere s = §2 /(2a). Direct integration of eq. (C.2), insertion of fi from the 



normahzation eq. (8.17), and some algebra gives the virial equation 



= -2Ev + 3Ei + bEi2 



Xo 



{4 



3A(xo cos Xo — sin xq))^ 



(C.3) 



We immediately see that this is in fact equivalent to eq. (8.15). There- 



fore, the solution we have explicitly found above minimizes the energy func- 
tional with boundary conditions p[xo) = VxP{xq) = 0. More generally, 
when we solved the MGP without considering the boundary terms in sec- 



tion (8.3), we found a one-parameter family of solutions (parametrized by 
A) . The virial theorem is merely a constraint on A for obtaining a minimum 
of E. 



120 



Bibliography 



[Adh02] S. K. Adhikari, Mean-field description of collapsing and explod- 
ing Bose-Einstein condensates, Phys. Rev. A, 66 (2002) 013611. 

[ADL04] E. Altman, E. Dcmler, and M. D. Lukin, Probing many-body 
states of ultracold atoms via noise correlations, Phys. Rev. A, 
70 (2004) 013603. 

[AEM+95] M. H. Anderson, J. R. Ensher, M. R. Matthews, C. E. Wieman, 
and E. A. Corneh, Observation of Bose-Einstein condensation 
in a dilute atomic vapor. Science, 269 (1995) 198. 

[AFT97] A. E. A. Amorim, T. Frederico, and L. Tomio, Universal aspects 
of Efimov states and light halo nuclei, Phys. Rev. C, 56 (1997) 
2378(R). 

[AG73] R. D. Amado and F. C. Greenwood, There is no Efimov effect 
for four or more particles, Phys. Rev. D, 7 (1973) 2517. 

[AN71] R. D. Amado and J. V. Noble, On Efimov's effect: A new 
pathology of three-particle systems, Phys. Lett. B, 35 (1971) 25. 

[AN72] R. D. Amado and J. V. Noble, Efim,ov's effect: A new pathology 
of three-particle systems. II, Phys. Rev. D, 5 (1972) 1992. 

[AS95] M. Abramowitz and I. A. Stegun (eds.), Handbook of Mathe- 
matical functions, Dover, New York, 9th ed. (1995). 

[ASRVKOl] J. R. Abo-Shaeer, C. Raman, J. M. Vogels, and W. Ketterle, 
Observation of vortex lattices in Bose-Einstein condensates. 
Science, 292 (2001) 476. 

[ATM+97] M. R. Andrews, C. G. Townsend, H.-J. Miesner, D. S. Dur- 
fee, D. M. Kurn, and W. Ketterle, Observation of interference 
between two Bose condensates. Science, 275 (1997) 637. 

[BayOl] G. Baym, Topics in the microscopic theory of Bose-Einstein 
condensates, J. Phys. B, 34 (2001) 4541. 

[BDZ08] I. Bloch, J. Dalibard, and W. Zwerger, Many-body physics with 
ultracold gases. Rev. Mod. Phys., 80 (2008) 885. 



121 



BIBLIOGRAPHY 



[BGOl] D. Blume and C. H. Greene, Quantum corrections to the 
ground-state energy of a trapped Bose-Einstein condensate: A 
diffusion Monte Carlo calculation, Phys. Rev. A, 63 (2001) 
063601. 

[BGH85] M. Brack, C. Guet, and H. B. Hakansson, Self consistent semi- 
classical description of average nuclear properties-a link be- 
tween microscopic and macroscopic models, Phys. Rep., 123 
(1985) 275. 

[BH03] E. Braaten and H.-W. Hammer, Universality in the three-body 
problem for ^Re atoms, Phys. Rev. A, 67 (2003) 042706. 

[BH06] E. Braaten and H.-W. Hammer, Universality in few-body sys- 
tems with large scattering length, Phys. Rep., 428 (2006) 259. 

[BHEOO] I. Bloch, T. W. Hansch, and T. Esshnger, Measurement of the 
spatial coherence of a trapped Bose gas at the phase transition. 
Nature, 403 (2000) 166. 

[BHK03] E. Braaten, H.-W. Hammer, and M. Kusunoki, Universal equa- 
tion for Efimov states, Phys. Rev. A, 67 (2003) 022505. 

[BHvKOO] P. F. Bedaque, H. W. Hammer, and U. van Kolck, Effective 
theory of the triton, Nucl. Phys. A, 676 (2000) 357. 

[BJ03] B. H. Bransden and C. J. Joachain, Physics of Atoms and 
Molecules, Pearson Education, 2nd ed. (2003). 

[BJK05] G. M. Bruun, A. D. Jackson, and E. E. Kolomeitsev, Multi- 
channel scattering and Feshbach resonances: Effective theory, 
phenomenology, and many-body effects, Phys. Rev. A, 71 (2005) 
052713. 

[BKK+05] R. Briihl, A. Kahnin, O. Kornilov, J. P. Toennies, G. C. 

Hegerfeldt, and M. StoU, Matter wave diffraction from an 
inclined transmission grating: Searching for the elusive ^He 
trimer Efimov state, Phys. Rev. Lett., 95 (2005) 063002. 

[BP96] G. Baym and C. .J. Pcthick, Ground-state properties of magnet- 
ically trapped Bose- condensed rubidium gas, Phys. Rev. Lett., 
76 (1996) 6. 

[BSH97] C. C. Bradley, C. A. Sackett, and R. G. Hulet, Bose-Einstein 
condensation of lithium: Observation of limited condensate 
number, Phys. Rev. Lett., 78 (1997) 985. 

[CGJT09] C. Chin, R. Grimm, P. Julienne, and E. Tiesinga, Feshbach res- 
onances in ultracold gases, arXiv:0812.1496v2 [cond-mat. other]. 



122 



BIBLIOGRAPHY 



[CHM+02] S. CowcU, H. Heiselberg, 1. E. Mazets, J. Morales, V. R. Pand- 
haripandc, and C. J. Pethick, Cold Bose gases with large scat- 
tering lengths, Phys. Rev. Lett., 88 (2002) 210403. 

[CKR+07] M. Combescure, A. Khare, A. K. Raina, J.-M. Richard, and 
C. Weydert, Level rearrangement in exotic atoms and quantum 
dots, Int. J. Mod. Phys. B, 21 (2007) 3765. 

[CKT+03] N. R. Clausscn, S. J. J. M. F. Kokkelmans, S. T. Thompson, 
E. A. Donley, E. Hodby, and C. E. Wieman, Very-high-precision 
bound-state spectroscopy near a ^^Rb Feshbach resonance, Phys. 
Rev. A, 67 (2003) 060701. 

[CMP07] A. Collin, P. Massignan, and C. J. Pethick, Energy- dependent 
effective interactions for dilute many-body systems, Phys. Rev. 
A, 75 (2007) 013615. 

[CPFG91] C. R. Chen, G. L. Payne, J. L. Friar, and B. F. Gibson, Nd 
zero-energy scattering, Phys. Rev. C, 44 (1991) 50. 

[Dan61] G. S. Danilov, On the three-body problem in the case of short- 
range forces, Zh. Eksp. Teor. Fiz., 40 (1961) 498. 

[DCC+01] E. A. Donley, N. R. Claussen, S. L. Cornish, J. L. Roberts, 
E. A. Cornell, and C. E. Wieman, Dynamics of collapsing and 
exploding Bose-Einstein condensates. Nature, 412 (2001) 295. 

[DE05] J. P. D'Incao and B. D. Esry, Scattering length scaling laws 
for ultracold three-body collisions, Phys. Rev. Lett., 94 (2005) 
213201. 



[DE06] J. P. D'Incao and B. D. Esry, Enhancing the observability of 
the Efimov effect in ultracold atomic gas mixtures, Phys. Rev. 
A, 73 (2006) 030703. 

[DGOl] J. L. DuBois and H. R. Clyde, Bose-Einstein condensation in 
trapped bosons: A variational Monte Carlo analysis, Phys. Rev. 
A, 63 (2001) 023602. 

[DGPS99] F. Dalfovo, S. Giorgini, L. P. Pitaevskii, and S. Stringari, The- 
ory of Bose-Einstein condensation in trapped gases. Rev. Mod. 
Phys., 71 (1999) 463. 

[DPS96] F. Dalfovo, L. Pitaevskii, and S. Stringari, Order parameter at 
the boundary of a trapped Bose gas, Phys. Rev. A, 54 (1996) 
4213. 



123 



BIBLIOGRAPHY 



[DRN07] A. M. Dudarev, M. G. Raizen, and Q. Niu, Quantum many-body 
culling: Production of a definite number of ground-state atoms 
in a Bose-Einstein condensate, Phys. Rev. Lett., 98 (2007) 
063001. 

[DvSG09] J. P. D'Incao, J. von Stecher, and G. H. Greene, Universal four- 
boson states in ultracold molecular gases: Resonant effects in 
dimer-dimer collisions, Phys. Rev. Lett., 103 (2009) 033004. 

[DZF+07] C. D'Errico, M. Zaccanti, M. Fattori, G. Roati, M. Inguscio, 
G. Modugno, and A. Simoni, Feshbach resonances in ultracold 
39k, New J. Phys., 9 (2007) 223. 

[Efi70] V. Efimov, Energy levels arising from resonant two-body forces 

in a three-body system., Phys. Lett. B, 33 (1970) 563. 

[Efi71] V. N. Efimov, Weakly-bound states of three resonantly- 
interacting particles, Sov. J. NucL Phys., 12 (1971) 589. 

[Efi73] V. Efimov, Energy levels of three resonantly interacting parti- 
cles, Nucl. Phys. A, 210 (1973) 157. 

[Efi90] V. Efimov, Is a qualitative approach to the three-body problem 
useful?, Gomm. Nucl. Part. Phys., 19 (1990) 271. 

[Efi91] V. Efimov, Force-range correction in the three-body problem: 
Application to three-nucleon systems, Phys. Rev. G, 44 (1991) 
2303. 

[EG06] B. D. Esry and G. H. Greene, Quantum physics: A menage a 
trois laid bare, Nature, 440 (2006) 289. 

[EGB99] B. D. Esry, G. H. Greene, and J. P. Burke, Recombination of 
three atoms in the ultracold limit, Phys. Rev. Lett., 83 (1999) 
1751. 

[Ein24] A. Einstein, Quantum theory of monatomic ideal gases, 
Sitzungsbcrichtc dcr Preussichcn Akadcmic der Wissenschaften 
Physikahsch-Mathematische Klassc (1924) 261. 

[ELG96] B. D. Esry, C. D. Lin, and C. H. Greene, Adiabatic hyperspher- 
ical study of the helium trimer, Phys. Rev. A, 54 (1996) 394. 

[FJ93] D. V. Fcdorov and A. S. Jensen, Efimov effect in coordinate 

space Faddeev equations, Phys. Rev. Lett., 71 (1993) 4103. 

[FJOla] D. V. Fedorov and A. S. Jensen, Correlation-induced collapse 
of many-body systems with zero-range potentials, Phys. Rev. A, 
63 (2001) 063608. 



124 



BIBLIOGRAPHY 



[FJOlb] D. V. Fedorov and A. S. Jensen, Regularization of a three-body 
problem with zero-range potentials, J. Phys. A, 34 (2001) 6003. 

[FJ02] D. V. Fedorov and A. S. Jensen, Regularized zero-range model 

and an application to the triton and the hypertriton, Nucl. Phys. 
A, 697 (2002) 783. 

[FJR94] D. V. Fedorov, A. S. Jensen, and K. Riisager, Efimov states in 
halo nuclei, Phys. Rev. Lett., 73 (1994) 2817. 

[FJT08] D. V. Fedorov, A. S. Jensen, and M. Th0gersen, Bose-Einstein 
condensates and Efimov states in trapped many-boson systems, 
Few-Body Sys., 43 (2008) 69. 

[FJT+09] D. V. Fedorov, A. S. Jensen, M. Th0gersen, E. Garrido, and 

R. de Diego, Calculating few-body resonances using an oscilla- 
tor trap. Few- Body Sys., 45 (2009) 191. 

[FKB+09] F. Ferlaino, S. Knoop, M. Berninger, W. Harm, J. P. D'Incao, 
H.-C. Nagerl, and R. Grimm, Evidence for universal four-body 
states tied to an Efimov trimer, Phys. Rev. Lett., 102 (2009) 
140401. 

[FRS79] A. C. Fonseca, E. F. Redish, and P. E. Shanley, Efimov effect in 
an analytically solvable model, Nucl. Phys. A, 320 (1979) 273. 

[FRS96] P. O. Fedichev, M. W. Reynolds, and G. V. Shlyapnikov, Three- 
body recombination of ultracold atoms to a weakly bound s level, 
Phys. Rev. Lett., 77 (1996) 2921. 

[FTDA99] T. Frederico, L. Tomio, A. Dclfino, and A. E. A. Amorim, Scal- 
ing limit of weakly bound triatomic states, Phys. Rev. A, 60 
(1999) 9(R). 

[FWG03] H. Fu, Y. Wang, and B. Gao, Beyond the Fermi pseudopoten- 
tial: A modified Gross-Pitaevskii equation, Phys. Rev. A, 67 
(2003) 053612. 

[Gaj06] M. Gajda, Criterion for Bose-Einstein condensation in a har- 
monic trap in the case with attractive interactions, Phys. Rev. 
A, 73 (2006) 023603. 

[GBC99] S. Giorgini, J. Boronat, and J. Casulleras, Cround state of a 
homogeneous Bose gas: A diffusion Monte Carlo calculation, 
Phys. Rev. A, 60 (1999) 5129. 

[GFTOl] A. Gammal, T. Frederico, and L. Tomio, Critical number of 
atoms for attractive Bose-Einstein condensates with cylindri- 
cally symmetrical traps, Phys. Rev. A, 64 (2001) 055602. 



125 



BIBLIOGRAPHY 



[GME08] A. O. Gogolin, C. Mora, and R. Egger, Analytical solution of 
the bosonic three-body problem, Phys. Rev. Lett., 100 (2008) 
140404. 

[GPS08] S. Giorgini, L. P. Pitaevskii, and S. Stringari, Theory of ultra- 
cold atomic Fermi gases, Rev. Mod. Phys., 80 (2008) 1215. 

[IGO+04] S. Inouye, J. Goldwin, M. L. Olsen, C. Ticknor, J. L. Bohn, and 
D. S. Jin, Observation of heteronuclear Feshbach resonances in 
a mixture of bosons and fermions, Phys. Rev. Lett., 93 (2004) 
183201. 

[JF03] A. S. Jensen and D. V. Fedorov, Efimov states in asymmetric 

systems, Europhys. Lett., 62 (2003) 336. 

[JHP02] S. Jonsell, H. Heiselberg, and C. J. Pethick, Universal behavior 
of the energy of trapped few-boson systems with large scattering 
length, Phys. Rev. Lett., 89 (2002) 250401. 

[JKTF07] A. S. Jensen, T. Kjaargaard, M. Th0gersen, and D. V. Fedorov, 
Eigenvalues of the one-body density matrix for correlated con- 
densates, Nucl. Phys. A, 790 (2007) 723c. 

[JNW06] L. M. Jensen, H. M. Nilsen, and G. Watanabe, BCS-BEC 
crossover in atomic Fermi gases with a narrow resonance, Phys. 
Rev. A, 74 (2006) 043608. 

[Jon04] S. Jonsell, Universality of the three-boson system close to a 
Feshbach resonance, J. Phys. B, 37 (2004) S245. 

[JRFG04] A. S. Jensen, K. Riisager, D. V. Fedorov, and E. Garrido, 
Structure and reactions of quantum halos. Rev. Mod. Phys., 
76 (2004) 215. 

[KFM+09] S. Knoop, F. Ferlaino, M. Mark, M. Berninger, H. Schobel, H.- 

C. Nagerl, and R. Grimm, Observation of an Efimov-like trimer 
resonance in ultracold atomdimer scattering. Nature, 5 (2009) 
227. 

[KGJ06] T. Kohler, K. Goral, and P. S. Julienne, Production of cold 
molecules via magnetically tunable Feshbach resonances. Rev. 
Mod. Phys., 78 (2006) 1311. 

[KH09a] S. Kreuzer and H.-W. Hammer, Efimov physics in a finite vol- 
ume, arXiv:0811.0159vl [nucl-th]. 

[KH09b] S. Kreuzer and H.-W. Hammer, Efimov physics in a finite vol- 
ume, Phys. Lett. B, 673 (2009) 260. 



126 



BIBLIOGRAPHY 



[KMW+06] T. Kracmcr, M. Mark, P. Waldburger, J. G. Danzl, C. Chin, 
B. Engeser, A. D. Lange, K. Pilch, A. Jaakkola, H.-G. Naegerl, 
and R. Grimm, Evidence for Efimov quantum states in an ul- 
tracold gas of cesium atoms, Nature, 440 (2006) 315. 

[LDD77] T. K. Lim, S. K. Duffy, and W. G. Damer, Efimov state in the 
^He tnmer, Phys. Rev. Lett., 38 (1977) 341. 

[LKJ07] M. D. Lee, T. Kohler, and P. S. Julienne, Excited Thomas- 
Efimov levels in ultracold gases, Phys. Rev. A, 76 (2007) 012720. 

[LMF+08] T. Lahaye, J. Metz, B. Frohlich, T. Koch, M. Meister, A. Gries- 
maier, T. Pfau, H. Saito, Y. Kawaguchi, and M. Ueda, d-Wave 
collapse and explosion of a dipolar Bose-Einstein condensate, 
Phys. Rev. Lett., 101 (2008) 080401. 

[LPS97] E. Lundh, G. J. Pethick, and H. Smith, Zero-temperature prop- 
erties of a trapped Bose-condensed gas: Beyond the Thomas- 
Fermi approximation, Phys. Rev. A, 55 (1997) 2126. 

[MFK+00] S. Moszkowski, S. Fleck, A. Krikeb, L. Theui31, J.-M. Richard, 
and K. Varga, Binding three or four bosons without bound sub- 
systems, Phys. Rev. A, 62 (2000) 032504. 

[MHUB06] E. J. Mueller, T.-L. Ho, M. Ueda, and G. Baym, Fragmentation 
of Bose-Einstein condensates, Phys. Rev. A, 74 (2006) 033612. 

[MMT07] D. Mozyrsky, I. Martin, and E. Timmermans, Coherent meso- 
scopic quantum tunneling in boson-fermion mixtures, Phys. 
Rev. A, 76 (2007) 051601. 

[MRD+09] N. P. Mehta, S. T. Rittenhouse, J. P. D'Incao, J. von Stecher, 

and C. H. Greene, A general theoretical description of N -body 
recombination, arXiv:0903.4145v3 [physics. atom-ph]. 

[NFJ97] E. Nielsen, D. V. Fedorov, and A. S. Jensen, Three-body halos 
in two dimensions, Phys. Rev. A, 56 (1997) 3287. 

[NFJ98] E. Nielsen, D. V. Fedorov, and A. S. Jensen, The structure of 
the atomic helium trimers: Halos and Efimov states, J. Phys. 
B, 31 (1998) 4085. 

[NFJ99] E. Nielsen, D. V. Fedorov, and A. S. Jensen, Efimov states in 
external fields, Phys. Rev. Lett., 82 (1999) 2844. 

[NFJGOl] E. Nielsen, D. V. Fedorov, A. S. Jensen, and E. Garrido, The 
three-body problem with short-range interactions, Phys. Rep., 
347 (2001) 373. 



127 



BIBLIOGRAPHY 



[NG99] M. Naraschewski and R. J. Glauber, Spatial coherence and den- 
sity correlations of trapped Bose gases, Phys. Rev. A, 59 (1999) 
4595. 

[NM99] E. Nielsen and J. H. Macek, Low-energy recombination of iden- 
tical bosons by three-body collisions, Phys. Rev. Lett., 83 (1999) 
1566. 

[PC97] D. R. Phillips and T. D. Cohen, How shoH is too shoH? Con- 
straining zero-range interactions in nucleon-nucleon scattering, 
Phys. Lett. B, 390 (1997) 7. 

[PDJ+09] S. E. Pollack, D. Dries, M. Junker, Y. P. Chen, T. A. Corcovilos, 

and R. G. Hulet, Extreme tunability of interactions in a 7Li 
Bose-Einstein condensate, Phys. Rev. Lett., 102 (2009) 090402. 

[PJP09] L. Platter, C. Ji, and D. R. Phillips, Range corrections to three- 
body observables near a Feshbach resonance, Phys. Rev. A, 79 
(2009) 022702. 

[Pla09] L. Platter, Universality and leading corrections in few-body sys- 
tems, Few-Body Sys., 45 (2009) 211. 

[P056] O. Penrose and L. Onsagcr, Bose-Einstein condensation and 
liquid helium, Phys. Rev., 104 (1956) 576. 

[PPOO] C. J. Pethick and L. P. Pitaevskii, Criterion for Bose-Einstein 
condensation for particles in traps, Phys. Rev. A, 62 (2000) 
033609. 

[PS02] C. J. Pethick and H. Smith, Bose-Einstein Condensation in 
Dilute Gases, Cambridge University Press (2002). 

[PS03] L. Pitaevskii and S. Stringari, Bose-Einstein Condensation 
(The International Series of Monographs on Physics), Oxford 
University Press, USA (2003). 

[RCC+01] J. L. Roberts, N. R. Claussen, S. L. Cornish, E. A. Donley, 

E. A. Cornell, and C. E. Wicman, Controlled collapse of a Bose- 
Einstein condensate, Phys. Rev. Lett., 86 (2001) 4211. 

[RF94] J.-M. Richard and S. Fleck, Limits on the domain of coupling 
constants for binding N-body system,s with no bound subsys- 
tems, Phys. Rev. Lett., 73 (1994) 1464. 

[RFOl] R. Roth and H. Feldmeier, Effective s- and p-wave contact in- 
teractions in trapped degenerate Fermi gases, Phys. Rev. A, 64 
(2001) 043603. 



128 



BIBLIOGRAPHY 



[RHBE95] P. A. Ruprccht, M. J. Holland, K. Burnett, and M. Edwards, 
Time-dependent solution of the nonlinear Schrodinger equation 
for Bose-condensed trapped neutral atoms, Pliys. Rev. A, 51 
(1995) 4704. 

[RZCD+07] G. Roati, M. Zaccanti, J. C. C. D'Errico, M. Modugno, A. Si- 
moni, M. Inguscio, and G. Modugno, ^^K Bose-Einstein con- 
densate with tunable interactions, Phys. Rev. Lett., 99 (2007) 
010403. 

[SEGB02] H. Suno, B. D. Esry, C. H. Greene, and J. P. Burke, Three-body 
recombination of cold helium atoms, Phys. Rev. A, 65 (2002) 
042725. 

[SFJ02] O. Sorensen, D. V. Fedorov, and A. S. Jensen, Correlated 
trapped bosons and the many-body Efimov effect, Phys. Rev. 
Lett., 89 (2002) 173002. 

[SFJ03] O. Sorensen, D. V. Fedorov, and A. S. Jensen, Correlated N- 
boson system,s for arbitrary scattering length, Phys. Rev. A, 68 
(2003) 063618. 

[SFJ05] H. H. Sorensen, D. V. Fedorov, and A. S. Jensen, Correlated 
gaussian method for dilute bosonic systems, AIP Conf. Proc, 
777 (2005) 12. 

[SJ87] P. J. Siemens and A. S. Jensen, Elements of Nuclei, Addison- 

Wesley (1987). 

[Sky56] T. H. R. Skyrme, The nuclear surface, Phil. Mag., 1 (1956) 
1043. 

[Sor05] H. H. Sorensen, Correlations in many-body systems with the 
Stochastic Variational Method, Master's thesis, University of 
Aarhus (2005), arXiv:0502126vl [cond-mat.soft]. 

[SSJF05a] T. Sogo, O. S0rensen, A. S. Jensen, and D. V. Fedorov, Semi- 
analytic solution to the N -boson problem with zero-range inter- 
actions, Europhys. Lett., 69 (2005) 732. 

[SSJF05b] T. Sogo, O. S0rensen, A. S. Jensen, and D. V. Fedorov, The 
zero-range approximation applied to the N-boson problem, J. 
Phys. B, 38 (2005) 1051. 

[Sto97] H. T. C. Stoof, Macroscopic quantum tunneling of a Bose con- 
densate, J. Stat. Phys., 87 (1997) 1353. 

[SV98] Y. Suzuki and K. Varga, Stochastic Variational Approach to 
Quantum-Mechanical Few-Body Problems, Springer (1998). 



129 



BIBLIOGRAPHY 



[SZS+04] C. A. Stan, M. W. Zwierlein, C. H. Schunck, S. M. F. Raupach, 
and W. Ketterle, Observation of Feshbach resonances between 
two different atomic species, Phys. Rev. Lett., 93 (2004) 143001. 

[TFJ07] M. Th0gersen, D. V. Fedorov, and A. S. Jensen, Trapped Bose 
gases with large positive scattering length, Europhys. Lett., 79 
(2007) 40002. 

[TFJOSa] M. Th0gersen, D. V. Fedorov, and A. S. Jensen, N-body Efimov 
states of trapped bosons, Europhys. Lett., 83 (2008) 30012. 

[TFJOSb] M. Th0gersen, D. V. Fedorov, and A. S. Jensen, Spatial corre- 
lations in Bose gases, AIP Conf. Proc, 998 (2008) 85. 

[TFJ08c] M. Th0gersen, D. V. Fedorov, and A. S. Jensen, Universal prop- 
erties of Efimov physics beyond the scattering length approxi- 
mation, Phys. Rev. A, 78 (2008) 020501 (R). 

[TFJ+09] M. Th0gerscn, D. V. Fedorov, A. S. Jensen, B. D. Esry, and 
Y. Wang, Conditions for Efimov physics for finite range poten- 
tials, Phys. Rev. A, 80 (2009) 013608. 

[Tho35] L. H. Thomas, The interaction between a neutron and a proton 
and the structure of H^, Phys. Rev., 47 (1935) 903. 

[TZJ09] M. Th0gersen, N. Zinner, and A. S. Jensen, Thomas-Fermi 
approximation for a condensate with higher-order interactions. 
Submitted to Phys. Rev. A, arXiv:0907.5349 [cond-mat. quant- 
gas]. 

[UL98] M. Ueda and A. J. Lcggctt, Macroscopic quantum tunneling of 
a Bose-Einstein condensate with attractive interaction, Phys. 
Rev. Lett., 80 (1998) 1576. 

[VHBL+98] L. Vestergaard Hau, B. D. Busch, C. Liu, Z. Dutton, M. M. 

Burns, and J. A. Golovchenko, Near-resonant spatial images of 
confined Bose-Einstein condensates in a 4-Dee magnetic bottle, 
Phys. Rev. A, 58 (1998) R54. 



[vSDG09] J. von Stecher, J. P. DTncao, and C. H. Greene, Four-body 
legacy of the Efimov effect. Nature Phys., 5 (2009) 417. 

[WGS98] N. K. Wilkin, J. M. F. Gunn, and R. A. Smith, Do attractive 
bosons condense?, Phys. Rev. Lett., 80 (1998) 2265. 

[Wig55] E. P. Wigner, Lower limit for the energy derivative of the scat- 
tering phase shift, Phys. Rev., 98 (1955) 145. 



130 



BIBLIOGRAPHY 



[WLO+09] A. N. Wenz, T. Lompe, T. B. Ottenstein, F. Serwane, G. Ziirn, 
and S. Jochim, A universal trimer in a three- component fermi 
gas, arXiv:0906.4378vl [cond-mat .quant-gas]. 

[Yan62] C. N. Yang, Concept of off-diagonal long-range order and the 
quantum phases of liquid He and of superconductors, Rev. Mod. 
Phys., 34 (1962) 694. 

[ZDD+09] M. Zaccanti, B. Deissler, G. D'Errico, M. Fattori, M. Jona- 
Lasinio, S. Miiller, G. Roati, M. Inguscio, and G. Modugno, 
Observation of an Efimov spectrum in an atomic system. Na- 
ture Physics. 

[ZT09a] N. Zinner and M. Th0gersen, Higher-order terms in an attrac- 
tively interacting EEC: Mean-field stability and macroscopic 
tunneling lifetimes, In preparation. 

[ZT09b] N. Zinner and M. Th0gersen, Stability of a Bose-Einstein con- 
densate with higher-order interactions near a Feshbach res- 
onance, Phys. Rev. A (in press), arXiv:0903.2261v2 [cond- 
mat. other]. 



131 



